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On the Quantum Theory of Reciprocal Fields and the 
Correspondence Principle 


By J.- RAY SKI 


Department of Physics, Nicolas Copernicus University, Torun, Poland 
Communicated by M. Born; MS. received 18th May 1951 


ABSTRACT. A method of quantization of non-local fields whose interaction is expressible 
by means of integral equations is given. ‘The S-matrix is constructed by Yang’s method. 
| The non-local field theory is free of the usual convergence difficulties. It is possible to 
% satisfy generally either the postulate of reciprocity or the postulate of correspondence, but it 
¥ seems impossible to satisfy both principles simultaneously. 


§1. GENERAL INTRODUCTION AND SUMMARY 

HE ideas of non-localizability and reciprocity (originated by Born (Born 

| and Peng 1944, Born and Green 1949) and developed further by Yukawa 

& (1950)) are very attractive for the following reasons: ‘he assumption that 
the wave function depends not only upon the operator of position but also upon 
the displacement operator constitutes a vast generalization. Some of the classical 
numbers have been replaced by operators, so that the evolution of quantum 
mechanics is seen as a natural continuation from classical theory. On the 
other hand, the postulate of reciprocity introduces a new symmetry into the 
fundamental laws of nature, similar to that introduced by restricted relativity: 
— as relativity takes account of the velocity of light and quantum mechanics of 
| Planck’s constant, so reciprocity introduces, in a very elegant way, a third 


/ fundamental constant, i.e. the universal length ». It is to be expected that the 


non-localizability will remove the usual convergence difficulties. 

The basic ideas seem very promising, but the theory is not yet sufficiently 
developed, and the decisive problems which need solution are: (i) how to interpret 
the non-local fields; (ii) how to take account of interaction between non-local 
fields; (i11) how the non-local field should be quantized ; (iv) whether it is possible 
to secure a reciprocity of interaction; (v) whether it is possible to secure corres- 
pondence with the local field theory. 

In order to avoid complications we investigate the above-mentioned problems 
on the simplest possible field model: a complex scalar field % with a scalar inter- 
action with a real scalar field ¢. In the frame of the traditional field theory these 
fields obey the equations 


(L—m*)p= —gdp, ...... (1) (Ll-)p= gb, ....-. (2) 
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where g is the coupling constant assumed small compared with unity (weak 
coupling). The equations (1) and (2) should be generalized for the case of a 
non-local field y. The real field ¢ is assumed as localizable for simplicity. 

We shall investigate first the non-local field % in vacuo. A correspondence 


cule wiil enable us to set up expressions which may be interpreted as densities. 
The quantization of the field im vacuo is performed by assuming suitable 


commutation relations in conformity with the postulate of reciprocity. It will 


be shown that, apart from a slight difference in densities, this theory of the non- | 
local and reciprocal field is fully equivalent to the usual quantum theory of a local — 


scalar field zm vacuo. 


In §3 we shall set up some integral field equations in order to take account of © 


a non-local interaction. In this way we avoid the infinitesimal point of view 


intimately connected with differential equations. The integral equations in — 


question constitute a straightforward generalization of the equations of Yang 
and Feldman (1950) adopted for the case of scalar fields, and may be written 
either in terms of ‘retarded’ or ‘advanced’ potentials. The form of the generalized 
A function appearing in the integral field equations need not be specified for 


the purposes of §3. The perturbed fields may be quantized by a method called — 
here a ‘direct quantization’, avoiding completely the Hamiltonian formalism and ~ 


the Schrédinger equation. The method of direct quantization consists simply in 
the assumption that the ingoing and outgoing waves obey the reciprocal com- 
mutation relations for the case of fields without interaction. The commutation 
relations for the perturbed waves then follow automatically from the field equations 
which may be solved by iteration. The S-matrix may be computed by means 
of the method of Yang and Feldman. General rules for an immediate trans- 
cription of the traditional S-matrix elements into the non-local form are inferred. 
The existence of an S-matrix proves the compatibility of the field equations and 
the consistency of the whole scheme. an 

In §§4 and 5 two alternative assumptions for the generalized A function are 
investigated. The first alternative provides us with a formalism satisfying 
generally the principle of correspondence: the formalism goes over into the usual 
one in the limit A—-0. However, the postulate of reciprocity is then violated 
(i.e. restricted to the case of a field im vacuo). On the other hand, the second 
alternative secures the reciprocity of interaction but violates the principle of 
correspondence in orders of approximation higher than g*. The fundamental 
discrepancy between the traditional field theory (or a theory standing in a full 
correspondence to the local theory) and the completely reciprocal field theory 
manifests itself in the extremely relativistic region, but is without influence 
upon effects whose main contributions lie in the non-relativistic region (as, for 
example, the Lamb shift). The non-local field theory is free from the usual 
convergence difficulties. 


§2. THE RECIPROCAL SCALAR FIELD IN VACUO 


In conformity with Born’s (Born and Peng 1944, Born and Green 1949) 
and Yukawa’s (1950) ideas of non-localizability and reciprocity, we assume that 
the wave function % depends not only upon the operator, of position x, but also 
on the displacement operator d,,: 


b=y(x,d), where [x}-d;\=6.,5°) 9 ee mere (3) 


[L? 


ey 
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1. We postulate the following reciprocal commutation relations : 
[x3 [x,,, ws]] =\{[d,, id; ws], See ene (4) 
' Par xu 0" eo. (5), Pata, ime. ew. (6) 
The relation (6) is reciprocal, due to (4). Equations (4), (5) and (6) are equivalent 
} to the equations of Yukawa (1950), with the only difference that his A has been 
replaced by our mA?. 

Following Yukawa, we use a representation in which x, are represented by 


) diagonal matrices. In this representation y is a matrix <x’||x”) which may be 
written also as a function of two variables, 


= (NAT hee eles Mea aA &) rea aw a (7) 
In this representation our equations (4), (5) and (6) assume the form 
‘ 7 Pab(x, 7) =A] gh(x, 7 \pgtetiole re eS et nese (4’) 
rs WaT) = 05 weiss (Ot (le = er aT yO. eae. Ss (6’) 
7 


Let us apply Fourier analysis to the wave function ¢(x,7). Its Fourier 
transform must be of the form 


VP T)=N(PWP) (TP)? + Per) na eee (8) 
in order to satisfy (4’) and (5’). M(p) is a normalizing factor extracted from 
%(p) for convenience. Equation (8) could still contain an arbitrary, though 
regular, function of 7,, but we purposely drop this factor in order to avoid higher 
spin values. As has been shown by Bloch (1950), this may be consistently 
achieved by means of a supplementary condition. If («,7) also satisfies (6’) 
then the Fourier transform ¥(p) contains a third delta function 6(p,2+m?), and 
the function 6(7,?+p,7A*) in (8) may be replaced by 8(7,2—m?A‘). 

In order to establish a correspondence and an interpretation of the non-local 
_ formalism we postulated in a previous paper (Rayski 1950) that for any gauge 
invariant polynomial a(x, 7) of wave functions and their derivatives the observable 


pees a(n) = | GU Gat \a 2 ae ted 9) DAA ets (9) 
The same postulate has been proposed independently by Moller (see Bloch 1950). 
By putting 3 

N(p) = { | d*78(p,7,)8(r,2-+ pe)| SRS (10) 
we get from (8) Wp) = | Tori ee i eT (11) 


which is of the same form as (9). 
With the aid of (9) we form the scalar density corresponding to 


AX ae eae ee Te Es (12) 


Remembering that the wave functions are actually matrices, we have 


(x) = | L7UX,7) = | d4r(x' |p| x") 


=|{ d'r d4x" (xe! |b* |x") <x" |]x"), vee (13) 

which may be also written as follows: 
o(m)=| | dirdtsy*(x—4r, s\Ww—ds7). ae, (14) 
Since Weer MIENU(K Fr) a | lean. (15) 
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where C is an arbitrary constant, the density (14) may be written also in the form 
o(x)=[[ dirdisp*(y,sh(y,7) tees (16) 
with VSx=s(r-s). ** 7 9 = eee (17) 
Thus the rule for constructing an observable scalar density at the point x is: 
multiply %(x,7) by %*(x,s), translate the point x to y and then integrate over 
d‘r d*s. "The same rule may be easily generalized for other densities, for instance 
the current density. The above considerations enable us to write immediately 
an equation corresponding to (2) in the case of a non-local field 4: we have only 
to replace the product appearing on the right-hand side of (2) by the expression 
(16). 
Quantization 
We quantize the non-local field in vacuo by assuming the following commuta- 
tion relations: 
[b(x.7) .*(x,s) | =2A(x—x,7, 5), | (18) 
[le 7), He’ r'V] = TEACH 9), PG sI=O. fo 


Ihe generalized A function is defined in the following way: 
1 . 
= 4 , 
IN CEO any | d*pexp (pA Dkr, 8) Wien ee (19) 


the Fourier transform being 

A(p,7,8)=A(p)p(P,1)P(Ps8), ne es (20) 
where A(p) is the Fourier transform of the usual A(x) function, well known from 
the theory of Schwinger (1949): 


Alp) = Zito) 2-72) Son Pye an meee (21) 
while p(p,7) is written for 
= p(P, r) =N(p)8(p,7,)8(7,7 AD) es (22) 
with N(p) given by (10). 
Since | arp pn=le Cy wae eee (23) 
we have I] dtr-d4sA(w,758\ = A(X) ae eee (24) 
We notice the following properties of the generalized Jordan—Pauli function: 
al A 7) S) = Nr A(x 7s) == Ags AAs) ee (255) 
0 0 ” 
Te am (755) = S, om N(& 0,5) = 0 rr ee (257) 
(S12 ACT Ss) =0, 2 eee (25%) 


which show that the A(x, 7,s) function is a reciprocal concept and that the 
commutation relations (18) are compatible with the field equations (4’), (5’) and 
(6’). We notice that A is an odd function of x,, 

AC 0,755) == A(x575 5) ee ree (26) 
as usual, but it should be stressed that the generalized A function is neither anti- 
symmetric in x) nor vanishes outside the light cone, in contradistinction to local 
field theory. ‘This new property is a characteristic feature of non-local fields 
and is closely connected with the finite dimensions of particles. Due to the 
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# existence of the delta functions 8(p,2+m?), 8(p,7,), 8(7,2+p,2A4), etc., in the 
) generalized A function the four-vectors 7, and s, are space-like and their compo- 
5) nents tend (separately) to zero for A->0. In this case the generalized A function 
© goes over into the traditional one, so that the correspondence with the local field 
| theory is fully taken into account. The non-commutability of wave functions 
j at points x, and x,’ lying outside each other’s light cones is appreciable only for 
1} distances comparable with the absolute length A. 

The wave functions satisfying the commutation relations (18) are of the form 
© (8), where (p) are now matrices identical with those used in the local theory of 
) a complex scalar field in vacuo. 

Applying the rule (9) we may introduce the charge and current density, 
} energy density, etc., quite analogous to the scalar density (16). It may be 
= easily shown that such densities obey the continuity equations, and it may be 
} verified that, although the densities differ from the traditional ones, the total 
| charge, energy and momentum operators (obtained by integrating the corres- 
ponding densities over three-dimensional space) are identical with the traditional 
4 ones. ‘Thus, except for a slight difference in densities, this theory of a non-local 
} and reciprocal field zm vacuo is fully equivalent to the usual local field theory and 
{ may be written in a canonical form. 


§3. THE RECIPROCAL FIELD, INCLUDING INTERACTION 

We have to determine the proper form of the field equations in the case of 
i interaction and to quantize the perturbed field %. The considerations of the 
) foregoing section led us already to the plausible assumption for the generalized 
} equation (2) describing the influence of the non-local field % upon the local field 4, 
(ae) A( os) U0) (27) 

) where v(x) is given by (16) and (17). The problem of a generalization of the 
» equation (1) describing the influence of the field ¢ upon the non-local % is more 
| difficult. Itis obvious, from the very outset, that we have to avoid the infinitesimal 
(| point of view to take full account of the non-local character of fields and their 
| interaction. The integral equations, avoiding completely the infinitesimal 
| approach, constitute the most natural form for the description of a non-local 
| interaction; therefore we look rather for integral than differential field equations. 
_ First of all we may replace (27) by the equivalent integral equation written either 
| by means of retarded or advanced potentials, 


Ha) = $a) + de [APA (xa WYA(y' 5), HO Vs + (28) 


where the product of the non-local wave functions has been written in a sym- 
_metrized form. ‘The symbol [dP’ stands for 


| dP’ = | | | d4x! d4y’ d4s! = | | | diy! dtr'dss', et (20) 

where x’ =y'+4(7' +5’), ae Est (17) 
ret =) 

while Neo A(X) Ea AN(G) ee eee Be ole ttre eevee (30) 


The A(x) functions involve the mass constant « of the real field g, but this cannot 
introduce any confusion with the delta functions A(x,7,s) connected with the 
non-local field and involving the mass constant m. 
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The ingoing and outgoing waves ¢'", 6°" obey the Schrédinger—Gordon 
equation im vacuo, in 
(T—12)4"*=0. nes 


Denoting half the sum of the ingoing and outgoing waves by the superscript 
(9), we have also 


$2) = 6x) + de [ PD ex UO SWOT). G2) 
‘Now we may postulate analogous integral equations accounting for the influence 
of the local field ¢(x) upon the non-local %(x*,7). We suppose that there exists 


a generalized A(x,7,s) function constituting a counterpart of the local A(x), so 
that we may introduce the ‘retarded’ and ‘advanced’ functions 


ret —_ 
AR4Y(2e, r, s) = ANG ifs 5) ae 4A(x, iP S). Soe (33) 


The form of the generalized A(x,r,s) need not be specified for the purpose of 
thissection. Atany rate, the application of the operator [_], — m* to the generalized 
A(x, 7,s) must not yield a (minus) Dirac delta function which was characteristic 
for the local theory. With the aid of the generalized A(x, 7, s) function we assume 
the following equation: 


Hy, 7) =, 7) +48 | APA(y—y', 7,8) Y's7), PH) ee (34) 


and similar equations written in terms of advanced or retarded potentials which 
need not be written down explicitly. ‘The ingoing and outgoing waves are 
assumed to obey the reciprocal equations in vacuo, i.e. the equations (4’), (5’) and 
(6’). The form of the equation (34), and similar equations for the retarded and 
advanced potentials, is suggested by analogy with (32). In particular, it should 
be noticed that we have had to use the translated variables y (instead of x) in order 
to obtain the corresponding density (16) at the point x. Since the application of 
the operator []—m? to the generalized A function is not supposed to yield a 
Dirac delta function, equation (34) is not reducible to a pure differential 
equation, in contradistinction to (32) which reduces to (27) due to the local 
character of the field ¢. 


Quantization 

The usual Hamiltonian formalism and the Schrédinger equation (which are 
typical for an infinitesimal point of view) are not adequate for the quantization of 
fields with a non-local interaction. On the other hand it is possible to quantize 
directly the fundamental integral equations written in terms of retarded or 
advanced potentials by assuming the interaction-free commutation relations for 
the ingoing and outgoing waves. The same relations hold for $® and 4. 
The commutation relations between the perturbed wave functions then follow 
automatically from the field equations and may be computed by iteration. 
In such a direct method of field quantization the role of the Heisenberg picture 
is predominant. It seems doubtful whether, in the case of non-local fields, it 
will be possible to go over to the Schrédinger picture or to the intermediate 
picture of 'Tomonaga and to restore a Schrédinger equation. 


By subtracting the equations with advanced potentials from those with 
retarded potentials we get 


L(yor)— oy, 7) =48 | dP'AQ—y's7,' 00,7) @)}, (35) 
G(x) — (a) =e | dP'ACe—v' YY, 8) WOT). eee (36) 
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« By applying the method described by Yang and Feldman we may construct the 
/ S-matrix for the non-local field theory. Since the ingoing and outgoing waves 
| satisfy the same commutation relations, there should exist a unitary matrix 
4) transforming simultaneously the ingoing into the outgoing waves 


=o, CoS OL Ns) MOR Lean (37) 
W with the same S, provided our equations are compatible with each other. From 
) (35), (36) and (37) we have 


W(y,7), S]= 4245, [dPAG-y.n)HO67) 0}, (38) 
($x), S]= —4{S, [ aPA@—e) YAO) YO Eee. (39) 


From (32), (34), (38) and (39) the S-matrix may be computed as a power series 
in the coupling constant g. But we may immediately guess the result: the 
reciprocal S-matrix elements may be obtained simply from the traditional ones 
by exactly the same substitutions which convert the corresponding field equations 
of the local theory into our (32), (34), (38) and (39). ‘Thus, in order to obtain 
the new S-matrix elements we have to replace in the traditional S-matrix elements 
every (x) and P(x) by PO(y,7) and ~*(y,s), every A(x’—x”) and 
A(x’ — x") involving the mass constant m by A(y’—y", 7’, s”) and A(y’—y", 7’, s”) 
while the local 4 (x) and its A functions remain unchanged. Moreover, the 
integrals { d+x'd*x” ... have to be replaced by the threefold integrals [dP’ dP”... 
given by (29).t ‘Thus we have proved the existence of the S-matrix by construc- 
tion. ‘The existence of an S-matrix constitutes, at the same time, a proof of the 
compatibility of our quantized equations and of the consistency of the whole 
scheme. It offers at the same time a possibility of confirming the formalism 
by experiment. 


§4. THE CASE OF FULL CORRESPONDENCE AND 
RESTRICTED RECIPROCITY 
Our next task is to determine the form of the non-local substitute for the 
traditional A function. We cannot take over literally the usual expression of 
Schwinger of the A function 
er Someta ACY 8h ere. Se (40) 
and replace in it the A(x) function by its non-local counterpart given by equations 
(19) to (22) since the non-local A(x,7,s) does not vanish outside the light cone, 
in which case (40) is not invariant. One of the possibilities of introducing a 
generalized A(x,7,s) may be accomplished by the following procedure: we may 
first replace sgn, by the invariant function (x), which is unity in the future 
light cone, minus unity in the past light cone and zero outside the light cone, 


Lifton? <0, x50 
n(x) = OM Ome Ue aes (41) 
—1 for x7 <0, x) <0, 
and write in the case of a local field 
ONG Ve Pe 7) CVA C2 er Pra (42) 


This procedure of constructing an S-matrix has been applied by the author in a previous 
paper (Rayski 1951). 
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which is completely equivalent to (40) owing to the fact that the local A(x) vanishes | 
outside the light cone. Now we may postulate a covariant non-local A function, 


A(«,17,8)= —4n(x)A(w,7,5). ka nee (43) 


From (33) and (43) we easily see that the function A™*(x,7,s) yields contri- 
butions not only from the past light cone, where it is equal to — A(x,7,s), but also 
from the domain outside the light cone, where it is equal to — $A(x,7,s). Thus 
the superscript ‘ret’ is now justified merely by the fact that no contributions come 
from the future light cone. The contributions coming from outside the light 
cone mean an action propagating with higher velocity than that of light, but this 
does not contradict relativity since such actions cannot be used for sending signals. 
Such effects are limited to distances comparable with the length A, and thus are 
only a consequence of the finite dimensions of elementary particles in the interior 
of which we cannot possibly trace causality. 

Let us consider an S-matrix element. The variables 7’,r’,...5',s”... 
appearing in it are connected either with the unperturbed 4 functions or with 
the generalized A(x,7,s) functions. Since their Fourier transforms '(p,7) 
or A(p,7,s) contain time-like p,: p,? = —m?, the four-vectors 7, and s,, are space- 
like and their components tend to zero for A tending to 0. Therefore any 
S-matrix element goes over into the local one in the limit A=0*. Thus there 
exists a complete correspondence between this formalism and the local field 
theory. On the other hand, it is the unperturbed field (x, 7) only that satisfies 
the reciprocal equations (4’) and (5’), but not the perturbed field (x, 7) (given 
by (34)). In this formulation the postulate of reciprocity is restricted to non- 
local fields am vacuo. 


So.) LEE CASE OR AVE UiteR Ee CLPROCiMay AND 
RESTRICTED CORRESPONDENCE 


In this section we show that it is possible to guarantee the reciprocity of the 
perturbed wave function ¢, but, as we shall see later, in this case the postulate of 
correspondence cannot be satisfied in higher orders of approximation. The 
complete reciprocity may be achieved by means of an alternative definition of 
the generalized A function. Instead of introducing the generalized A function 
into (42) we may directly generalize the local A(«) function on the same lines as 
we have generalized the A function by means of equations (19) to (22). We 
assume that the generalized A is 


7 ee tee : 
ING IS) = ony | UPA(D, 15S) CXP\CDj,,), ise eee eee (44) 
where the Fourier transform may be written in the form } 
A(p,7,s)=A(P)p(P,7)—(P,5), nnn ae (45) 
A(p) being the Fourier transform of the usual A(x) function, well known from the 
theory of Schwinger, ~ 
y 6 A(p)=Pol(p2+m), aaa (46) 


* Born and Yang have informed the author that they have convinced themselves by a direct, though 
lengthy, calculation that the equations (38) and (39), which define S in the non-local theory, go 
over, for A+0, into the corresponding equations of the local theory. 
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) while p(p,7) is given by (22). Owing to the delta functions 6(7,2+p,?A4), 
6(p,7,,)).+- in the Fourier transforms (45), the new A(x, 7, s) satisfies the equations 


ee 2 sr, s)=G 2K, 5)50, ..--, (47) 
es pe 
a ace TiS) 25,5 Oke Aataes = 024 oe ncaa (48) 


| which show that the generalized A function is a reciprocal concept. It should be 
} noted that the application of the operator [],—m? to the generalized A does not 
* yield a Dirac delta function, in contradistinction to the local field theory, whence 
| the equation (34) is not reducible to a pure differential equation. 

| The unperturbed function $(y,7) appearing in (34) is assumed to obey 
} (4’), (5’) and (6’), while the new A function obeys (4’) and (5’). Therefore 
| the application of the operators 7,?—A‘L], or 7,(0/dy,) on the right-hand side of 
p (34) yields identically zero, whence the left-hand side yields zero also. In this 
} way it is proved that not bale (x, 7) but even the perturbed wave function 
) %(x,7) satisfies the reciprocal equations (4’) and (5’). By applying to (34) the 
| operator [],—m?, using (47), and changing the variables, we get an integro- 


differential equation (EL Smiinasilary, Ole > aa (49) 
) where f(x,7) is the following function of % and ¢: 


fle,r)=4g | dP(),+ 7,24 2m) — 2,7, 8)Ua',r), dla! + Hr’ +5')D. 
eee (50) 


| Equation (49) is reciprocal since ¢b(x,r) obeys (4’) and (5’) and f(x, 7) is self- 
) reciprocal. ‘Thus the perturbed wave function #(x,7r) satisfies the reciprocal 
) equations (4’), (5’) and (49). In this way the non-local field theory, including 
| interaction, is put into a completely reciprocal form. 

| We mentioned above that this completely reciprocal formulation violates, to 
 acertain extent, the principle of correspondence. This is due to the circumstance 
| that the Fourier components of the A(x, 7, s) function are not restricted to p,? = — m?, 
| but p,? may assume any value, even positive. For space-like p, the variables 
-r and s (which are connected with p, through the delta functions 6(p,7,,) and 
) 8(7,2+p,7A*), etc.) may be time-like and their components do not vanish in 
| the limit A 0, in contradistinction to the case investigated in § 4. 

Let us discuss the results of the above displayed formalism. It may be easily 
/ seen that the new S-matrix elements yield non-vanishing probabilities only for 
/ processes in which the total energy and momentum are conserved. But the 
probability amplitudes differ from the traditional ones by a set of factors of 


the type 
oe J(p, 9) = I EEN CASAS Os Ne Bec (51) 


/ where p, and q, are energy and momentum four-vectors for real or virtual states 
of the non-local particles. By introducing the notation 


Nisa) = | *r EXP (32947 (DM) +P A) avons (52) 
the factors (51) may be written, from (22), as 
J(b, 9) =N(d)N (2, 9). eres (03) 
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The computation of the integral (52) yields the following result: N(p,q) is a 
function of the argument, invariant and symmetrical in p and q, 


G2IND.G,) =P de be eee (54) 
of the form 
ee plea a 0 ; 
N(P, 9) = aye anya for azZ0, p, ae ee (55’) 
1 u 
N(p, 9) = Ze 4/ =aexp (4/—a) Vfor 720.0908 ae ee (539) 


If we compute directly the normalizing factor N(p) from (10) then we get simply 
<0 


N(p)= 41/27 for p2 a SY OL ee (56) 

Hence the factor (53) becomes 
J(p,q) = Sain Va for. 9.24.20,” seer (57’) 
J(p, 9) = = 7 oP (=4/=a), for -a <0 7 PRES ee (57") 


Thus J(p,q) is a convergence factor strongly diminishing the interaction for 
|a|->co. In consequence of the existence of the factors J(p,q) the formalisms 
discussed here are free of the usual convergence difficulties. 

The case a>0 occurs always when at least one of the four-vectors p,, q, is 
not space-like. (This is always the case for the formulation given in §4, where 
p, is always time-like.) In this case the convergence factor tends to unity for 
A—0, as it should do in a theory obeying the principle of correspondence. The 
case a<( is rather exceptional and may occur only for virtual processes in higher 
orders of approximation than g’, where both p, and g, may be space-like simul- 
taneously. In such cases we get a strange result violating the principle of 
correspondence: J(p,g)->00 for A+0. For A finite, but such p, and gq, that 
a->—0), the convergence factor J(p, g) tends also to infinity, so that we encounter 
an unexpected case of resonance. ‘The question whether this singularity does, 
or does not, introduce new divergences needs a special investigation. 

This strange feature of the reciprocal interaction may be partly removed by 
the following device: The normalizing factor N(p) given by (10) is identical with 
N(p, p) as defined by (52) and may be regarded as a limit 


Np) = lum NOP, 9). eee (58) 
Meera 


This limit transition is, however, ambiguous. If we perform the transition to: 
the limit from the domain a >0 we get, from (55’), 


| er ie “e 
NP ee NP. Daze > okies Meseeaes (59’) 
a—>+0 
but if we let gp so that a<0, then we get, from (55”), 


N(p)= lim N(p,q)=0. Betas (59") 
I=>D 


a—>-—0O 


—— 


SS eee 
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In view of the above ambiguity we may define N(p) appearing in (53) according 
to whether p and g appearing in (53) yield a>0 or a<0. Now the convergence 


factor becomes sin «/a 
for a>0, 


SCRE 2 ai alates Ai eRe. (60) 
0 fOr a-<0; 


Consequently, any virtual transition between two states characterized by space- 
like four-vectors of energy and momentum of virtual particles is made impossible 
and the formalism is free from any resonance for a->0 and from the usual 
divergences. 

The higher order corrections to scattering processes will be sensitive to the 
irregular behaviour of the S-matrix elements for a<0 which may constitute a 
possibility of experimental check of the reciprocal formalism. At any rate, all 
the effects restricting the principle of correspondence lie in the relativistic region 
(since small perturbations cannot change the time-like character of the four- 
vector of energy and momentum of the (virtual) particle). ‘Thus the correspondence 
principle is still valid in the non-relativistic limit. 


§6. CONCLUSION 

We have seen that it is possible to formulate a quantum theory of non-local 
scalar fields in full agreement either with the postulate of correspondence with the 
local quantum field theory or with Born’s postulate of reciprocity, but it seems 
impossible to satisfy both postulates simultaneously. If the reciprocity is quite 
generally taken into account then the correspondence is restricted, and vice versa. 

Finally we have to examine whether it is possible to apply similar ideas to more 
complicated fields and interactions, above all to quantum electrodynamics. It is 
difficult to decide if, in case of a spinor field zn vacuo, we should retain equations 
(4) and (5) and replace (6) by the Dirac equation, or whether we should look for 
a suitable equation of the first order to replace (4) and to secure full reciprocity 
with the Dirac equation. The first alternative spoiling some of the reciprocal 
aspects is simpler, and it seems questionable whether a general reciprocity should 
be maintained also in the case of spinning particles, as the reciprocity establishes 
a symmetry between the operators of position and displacement but not between 
position and the rotational degrees of freedom. 

The extension of the above-described formalism (without the replacement of 
(4) by an equation similar to the Dirac equation) lies at hand. However, some of 
the traditional features of quantum electrodynamics will be lost in consequence 
of the non-localizability. ‘The continuity equation for the (perturbed) current 
formed with the aid of the correspondence principle (9) cannot be maintained 
and the gauge invariance of the second kind will also be violated. In the opinion 
of the author the lack of gauge invariance and of the continuity equations for the 
perturbed current do not disqualify such a formalism of quantum electrodynamics. 
We should remember that the introduction of new universal constants is necessarily 
accompanied by some renunciations. The inclusion of the finite velocity of light 
necessitated the sacrifice of the notion of absolute time, the introduction of Planck’s 
constant, the sacrifice of the traditional concept of state and measurability. It 
seems quite natural that the incorporation of a fundamental length will necessitate 
also a renunciation of some of the notions which we were accustomed to and liked 
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so much. It seems that this time we have to abandon the details of kinematics 
of densities resulting in the continuity equations, as it is obviously impossible to 
trace the currents inside the charged particles. 

In spite of the lack of a continuity equation, the integral charge conservation 
law is guaranteed simply by the gauge invariance of the first kind. On the other 
hand, the gauge invariance of the second kind is by no means a necessary condition. 
for the vanishing of the rest mass of the photon. The value zero for the photon 
rest mass may be secured by other devices, e.g. by a realistic compensation by 
means of charged bosons or by a formal mass renormalization which, in the frame 
of a convergent formalism, becomes a mathematically correct procedure. In this 
respect the new formalism does not constitute any drawback in comparison with 
the traditional quantum electrodynamics which was gauge invariant only formally, 
so that we were also obliged to renormalize the (infinite!) self mass of the photon. 
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ABSTRACT. A classical treatment of the domain energy terms of a homogeneous 
ferromagnetic solid leads to a formula for the internal field contributions from these terms. 
With this result, modifications in the resonance condition of ferromagnetic resonance 
arising from self energy, exchange energy, magnetocrystalline anisotropy, and applied or 
intrinsic stress are obtained and are applied to various crystalline anisotropy and stress 
conditions of interest in ferromagnetic resonance experiments. Finally, the bearing of 
the results on the anomalous g-values obtained in resonance experiments is considered. 


Sa INTRO DIT CiOin 
HIS paper is chiefly concerned with the forces which influence the direction 
of the magnetization vector I of a ferromagnetic substance in an external 
magnetic field and with the effects of these forces upon some of the 
phenomena of ferromagnetic resonance. A macroscopic viewpoint is adopted 


* Now at Armour Research Foundation, Chicago, Illinois, U.S.A. 


Ferromagnetic Resonance and the Internal Field 969 


throughout, and all microscopic interactions are replaced by formulae for the 
| magnetic potential energy density of the material which involve the macroscopic 
» magnetization vector. 

The usual case considered in calculating theoretical magnetization curves 
) is that of an ellipsoidal ferromagnetic specimen in a static magnetic field (Stoner 
and Wohlfarth 1948); in this case the magnetization is uniform throughout the 
specimen and is time-independent. In ferromagnetic resonance experiments, 
_ however, a very small high-frequency oscillating magnetic field is applied to the 
+ ellipsoidal specimen, as well as a large static field oriented at right angles to the 
» oscillating field. ‘The magnetization of the body may then be neither space nor 
» time independent, although, because of the relative magnitudes of the static 
# and oscillating applied fields, variation of the magnetization vector with time and 
© position willbesmall. One of the principal problems of ferromagnetic resonance, 
however, is the determination of this small space-and-time dependence of the 
direction and magnitude of the magnetization vector. In the purely static case 
the direction of the magnetization vector is determined from the general condition 
that the total magnetic potential energy of the body be a minimum. Because 
the magnetization is here uniform throughout the specimen, minimization of 
» the expression for the magnetic energy density with respect to angle variables 
} allows the direction of the magnetization vector to be uniquely specified. 

In the present work the assumption will be made that the same general 
condition of minimum total magnetic potential energy can be applied to the 
oscillating field case when the magnetization of the body is homogeneous to 
yield the direction of the magnetization vector to a good order of approximation. 
When this minimization procedure is valid, the direction of the magnetization 
vector which is found is its instantaneous equilibrium direction I,, defined as 
the direction in which the magnetization vector would remain or reach quickly 
if the time variation of the oscillating field were suddenly halted at the instant 
under consideration (and if retarded fields were of negligible importance). 

If the oscillating field is of very high frequency, the actual instantaneous 
direction of the magnetization vector I may not coincide with the instantaneous 
equilibrium direction. The actual instantaneous direction and magnitude of the 
magnetization at a given point may be found, however, from the solution of the 
differential equation of motion satisfied by the magnetization vector; this 
equation involves the instantaneous internal, or local, field H' at the point in the 
body considered. ‘The instantaneous equilibrium magnetization vector is, by 
definition, parallel to the instantaneous internal field, since minimum magnetic 
energy of the specimen results when the magnetization vector lies along the 
internal field at every point in the body. 

In §2 the variational principle is applied to the total magnetic potential energy 
of a homogeneous non-conducting ellipsoidal ferromagnetic body to yield a 
general method for calculating the internal-field vector when an external magnetic 
field is applied, and it is shown how the instantaneous direction of the magneti- 
zation vector may be determined. ‘The contributions to the internal magnetic 
field from the various magnetic potential energy terms present in a ferromagnetic 
specimen are evaluated in §3. In §§4-6 the influence on the ferromagnetic 
resonance condition of the various energy terms is discussed in detail for macro- 
scopic single crystals and for polycrystalline aggregates. Finally, in §7 the 
bearing of the results of the paper upon the anomalous Landé g-factors found 
in resonance experiments is considered. 
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§2. CALCULATION OF THE INTERNAL FIELD 


Before deriving a general expression for the internal field H' in terms of the 
magnetic potential energy density, it is desirable to define further quantities and 
to consider the dependence of these quantities upon space and time. Since 
the externally applied field H is the sum of a static component Hy and a high- 
frequency oscillating component h, the internal field H' will be a corresponding 
sum of H,! and hi, where h/H, and hi/H,! are very much smaller than unity, 
H, and h will be independent of position within the body providing it is a non- 
conducting ellipsoid; if the ellipsoid is electrically conducting, h will decay 
with penetration depth because of the skin effect. Similarly, the magnetization 
vector I will be a sum of a static component I, and an oscillating part i. Because 
of the smallness of the oscillating field h compared with the static field, the 
oscillating component of the resulting magnetization will likewise be far smaller 
than the static component even at resonance, when the oscillating component is. 
a maximum. 

Finally, it is easy to show that the magnitude of the total magnetization 
vector is very nearly independent of time. The vector differential equation of 
motion of I, at a given point in the body, is 


7) ; 

; i oo eee. (1) 
where y=ge/2mc =gup/h is the magnetomechanical ratio, g the Landé splitting 
factor, and zz the Bohr magneton. Note that the g-factor is assumed, subject 
to experimental verification, to be a constant independent of field orientation, 
body shape, etc. Equation (1) is just an expression for the rate of change of 
angular momentum per unit volume of the system. It shows that I.dI/dt=0 
at all points in the magnetic body, and that therefore J is time-independent. 
However, a small additional term should be added to equation (1) in order to 
account for damping caused by interactions between spins in the magnetic 
material and between the spin system and the crystalline lattice (Kittel 1948, 
Bloembergen 1950). Such damping may destroy the time-independence of J, 
but since damping is usually found experimentally to be small, J still remains 
nearly time-independent; therefore the effect of damping will be neglected in 
the present treatment. It will now be shown how an expression for the internal 
field may be derived by extending a method used by Landau and Lifshitz (1935). 

_ There are several contributions to the total magnetic potential energy 
density which may influence the direction of the magnetization vector in a 
ferromagnetic body. The following will be considered here: (a) the magnetic 
potential energy density of the body in the external magnetic field, (5) the potential 
energy density arising from the magnetic self-energy of the body associated with 
its shape, (c) the Heisenberg exchange, or Weiss, energy density, (d) the magneto- 
crystalline anisotropy energy density, which is thought mainly to arise from 
spin-orbit coupling between neighbouring atoms, and (e) the strain magnetic 
energy density connected with magnetostriction and dependent upon the state | 
of strain of the body. | 

It should be noted that domain wall energy need not be explicitly considered | 
in the present treatment because the static external magnetic field will always 
be taken sufficiently large for the specimen to be considered a single domain. 
The domain wall energies and thicknesses for zero applied fields have recently | 
been considered by Lilley (1950) and by Kittel (1949 b). | 
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f The magnetic potential energy density of a uniformly magnetized magnetic 
_ body (non-conducting ellipsoid) may usefully be considered as a function of the 
) following variables*: E=E(I,X;,01/0X;,), where X; represents the three rectan- 
) gular coordinates. Hence the total potential energy of the body will be 


al 
{ E (LX, ax) Fhe 


{ Applying the condition that the total magnetic potential energy of the body 
} must be stationary with respect to small variations of the magnetization vector I, 
‘ one obtains 


{28 (t X, gy) a= | are Zax qanax | oe) 


Now the magnitude of I has been anne to be a constant independent of 
} time. Therefore the variation 6I is not completely arbitrary but must lie in the 
{ plane perpendicular to I,,. This is because the variation consists of small 
i changes in the direction of I from the instantaneous equilibrium direction I,,. 
¢ Since the variation is otherwise arbitrary in magnitude and direction, however, 
| the integrand of the above integral must itself be zero. One thus obtains 


(BERCLO20Ph me oA Dra he asi. (2) 


i) where ® is the differential vector operator 


(-a*?ax Lzaraxs]} 


Equation (2) can only be satisfied if the vector WE is parallel or antiparallel 
} to I,, or is zero, since SI, is perpendicular to I,,. However, ®E cannot, in 
* general, be zero as long as an external field is applied. Further, by definition, 
) I,, is always parallel to the internal field H'. Therefore (@Z) must be equal 
to Hi apart from a proportionality factor. This factor can be shown to be unity 
' by applying the operator ® to an energy density term for which the corresponding 
contribution to the internal field is already known. Hence 


TS ay ee i, eee ae eee (3) 


; and the instantaneous equilibrium direction of the magnetization vector at any 
| point in the body may be determined from the condition that 
| bine a 3 tall) ee cent a) 
) Since Hi itself will usually be a function of I and its individual components, 
' equation (4) must be used to determine the actual equilibrium direction of I. 
» A help in solving this vector equation is the auxiliary equation I, x Hy'=0, which 
» expresses the fact that the static component of the magnetization lies in the 
direction of the static component of the internal field. Equation (3) will be used 
| later to compute some of the contributions of the different energy densities to 
| the internal field. 
! When the applied frequency is very high, the instantaneous direction of the 
| magnetization vector will not coincide with its instantaneous equilibrium direction. 
| The actual instantaneous direction may be found using the equation of motion 
(equation (1)), which may be rewritten as 

lol 10: 


olan =i et tes =(Iy +i) x (Hy'+h') =I, x hit+ix Hy ...... (5) 


* The index here and in all further equations takes on values from 1 to 3. 
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since ix hi is a very small second-order term and may be neglected. Becauses 
ys lope 10? cycles/sec/oersted, the left side of equation (5) is very small unless 
the applied frequency is high. When the frequency is low Ix H'=0, and the 
instantaneous equilibrium and actual directions of the magnetization vector 
are always effectively the same. But when the frequency is high this is no longex: 
the case, and the direction of i must be found by solving equation (5). Sinces 
I, is effectively the saturation value of the magnetization and its direction is fixed} 
by the relation I, x H,'=0, the instantaneous direction of the total magnetization: 
vector I=I,+i may then be found. 


§3. INTERNAL FIELD CONTRIBUTIONS 

In this section the individual contributions to the total internal field from the) 
five potential energy densities mentioned in § 2 will be calculated using equation (3), 
The five energy density terms may be written as functions of the magnetization: 
vector as follows: 

(a) External-field energy density: F..,= —H.I. 

(b) Self-energy density (for a non-conducting ellipsoid only): E,=}=2 N*[?. 

(c) Exchange energy density (Stoner and Wohlfarth 1948)*: 


Foxe = ~tNgl I+ FE S(VE£). (V I,)] 


ol ol 
ay 
4N,1.1+3 cu | xx | 


(d) Crystalline anisotropy energy for single crystals (Becker and Doring 1939): 
Hexagonal crystals: E,= Ky,'(1—«,?)+ K,’(1 —«,?). ; 
Cubic crystals: EB, = Ky[a4?05? + a22097 + a47047] + Ko[la2%2057] 

= $Ky[1 — (a y+ a4 + «54)] + Keloy7a9"a9"]. 
(e) Strain energy density for cubic single crystals (Becker and Doring 1939)f: 
Eig = — ¥ Ajo BaP Py — 3g 11 {61%2 Pye + 9% Pog + O31 P 3}. 


In the above, «,, %) and « are the direction cosines of the magnetization vector 
with respect to the X,, X,, X; axes, taken along the cubic axes of the cubic crystal 
and taken with the Xj axis in the direction of the principal crystalline axis of the 
hexagonal crystal, so that «,=J,/J=(1.i,)/I, where i, is a unit vector along the 
X, axis, etc. The P,,’s are the components of the stress tensor referred to the 
same coordinate system, Ky’, K,’ and K,, K, are crystalline anisotropy constants 
which may be determined ee ere Recent values of these constants 
are quoted by Kittel (1949b). ‘The N%,’s are the shape demagnetization constants 
applying when the coordinate system is taken along principal axes of the ellipsoidal 
specimen. N, is the Weiss molecular field constant, and C is a constant 
approximately equal to $a°N,,, where a is the grating spacing of the material. 
The results of more accurate calculations of the quantity C are given by 
Lifshitz (1944) and by Stoner and Wohlfarth (1948). Finally, Aypp and A,,, are 
the saturation values of the magnetostriction constants of the ferromagnetic 
crystal in the [100] and [111] directions respectively. Numerical values of 
these constants for iron and nickel are given by USES and Doring (1939, 
pp. 277-280). 


* This formula has been derived only for substances having cubic structure. 
+ This expression is valid only for the region of elastic strain and for stress tensor elements 
independent of position (homogeneous stress). 
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Since the total energy density is to a very good approximation the sum of the 
( individual energy density terms, the contribution to the internal field of each 
| energy density term may be obtained by applying the operator ® to each term in 
turn. One then finds the following results: 


phd CRs OS ei er en ere 1 (6) 

(H',); SS Dy N* I; Caceoach (7) 

Hi... =N,I+ OSA od eee glee renee. ry (8) 

’ , 2 

Hi, = [= + = {1 _ (3) \ (3) ie tlexagondl-4 © ane ee (9) 

(A= “ I - a Ll21? (j,k, 2=1, 2,3 and permutations) 
Cibich oe) nes (10) 

; 3 
(H')j= RP = [(Aroo—Asa)Oje+ Aral Pyle we (11) 


¢ Some of these equations are written directly in terms of vectors, others in terms 
_of the jth component of the corresponding vector, in the interest of simplicity. 
From the above expressions the direction and magnitude of each contribution 
: to the total internal field may be obtained. Previously, Bozorth (1949) obtained 
| the approximate maximum magnitude of the internal field contributions arising 
» from several of the energy density terms considered here by equating the terms 
> separately to H'.I. Further, Kittel (1948) computed the quantity Ix H', for 
; the case of crystalline anisotropy for several special directions and planes in 
) single crystals. Finally, Brown (1940) has used a different variational approach 
to obtain general equations equivalent to (I) x H,');=0, expressing the balance 
of torques on the magnetization vector when no oscillating field is applied, and 
i including the static effect of all the potential energy terms considered here but 
' with these terms expressed somewhat differently. ‘The present results, however, 
| give both the magnitude and the general directional dependence of the various 
static and oscillating internal field contributions in convenient form. It should 
) be noted, however, that the internal field contributions due to crystalline anisotropy 
| and strain in cubic crystals are not unique, because the relation «7 +a)? +a3?=1 
may be used to transform the expressions for the energy densities to different 
| forms from those used here. Nevertheless, such indeterminacy does not affect 
» the term I x H' which determines the direction of I; therefore both the direction 
| of I and the effect of the above energy terms upon ferromagnetic resonance 
| phenomena remain unique. 


§4. FERROMAGNETIC RESONANCE EXPERIMENTS 


| The foregoing internal field contributions may be used to compute the 
| dependence of the magnetization curves of ferromagnetic single crystals along 
| different crystal directions upon ellipsoid shape, crystalline anisotropy and applied 
stress. Part of such a programme has been carried out by Stoner and Wohlfarth 
| (1948), who, however, used the method of direct minimization of the total magnetic 
| energy density to obtain the dependence of the magnetization components upon 
an applied static magnetic field. In solving such a problem, the actual internal 
| field contributions derived here are unnecessary, although they shed additional 
_ light on the physical situation considered. In considering, however, the influence 
| PROC. PHYS. SOC. LXIV, II—A 63 
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of ellipsoid shape, exchange anisotropy, crystalline anisotropy and strain upon 
ferromagnetic resonance phenomena, explicit expressions for the internal field 
contributions are useful. 

In a ferromagnetic resonance experiment two external magnetic fields are 
applied to an ellipsoidal ferromagnetic body. The specimen may be a single 
crystal or polycrystalline, non-conducting or conducting. The fields are usually 
taken along two mutually perpendicular principal axes of the ellipsoid. One of 
the fields is large and static but can be varied in magnitude over a wide range; 
the other is produced by incident microwave electromagnetic radiation and is 
very small in magnitude. In making ferromagnetic resonance experiments, the 
frequency of the incident radiation is commonly held constant and the magnitude 
of the static field varied. It is found that a maximum of the absorption of micro- 
wave power in the specimen occurs for a unique static field strength. This 
resonant field strength, H,", is that field strength which makes the effective Larmor 
frequency of the spins of the ferromagnetic material, (y/27)H, equal to the 
frequency of the applied radiation v. Since the effective field H* involves all 
the internal field contributions, theoretical determination of the resonant field 
strength requires information concerning such contributions. 

It is usually assumed in ferromagnetic resonance calculations that the externally 
applied static field is always of sufficient strength to cause the resulting static 
component of the total magnetization vector to be substantially equal to the 
saturation magnetization of the material and to lie in the direction of the static 
field. Crystalline anisotropy and applied stress may tend to pull the static 
magnetization vector away from the field direction, but if the field is sufficiently 
strong (as is usually the case experimentally) such effects may be neglected. 
Then any influence of shape, crystalline anisotropy, or strain on the resonance 
phenomena arises from the interaction between their internal field contributions 
and the very weak oscillating component of the magnetization vector produced 
by the oscillating magnetic field of the incident radiation. 

There are two particularly important quantities which may be obtained from 
the results of a ferromagnetic resonance experiment. ‘These are the Landé 
g-factor determined from the resonant magnetic field strength, and the damping 
factor which may be calculated from the shape of the microwave absorption 
resonance curve. ‘lhe importance of these two factors in giving insight into 
the coupling between ferromagnetic electrons and between these electrons and 
the crystal field has been recently discussed by Kittel (1949a) and Van Vleck 
(1950). Inthe present treatment we are especially concerned with the calculation 
of the Landé g-factor for different ferromagnetic specimen shapes, orientation of 
applied fields, crystalline nature and applied stress, although other quantities, 
such as crystalline-anisotropy and magnetostriction constants, may also be obtained 
from the analysis of resonance experiments. It is not expected that the g-factor 
of a given ferromagnetic material should depend appreciably upon any of the 
above conditions; however, its calculation requires knowledge of the effects of 
these factors upon the resonant magnetic field. Such knowledge may be 
embodied in the explicit form of H', and, given an experimental H,' for a specific 
microwave frequency, the corresponding g-value may be computed from the 
resonance condition v=(ge/4mmc)H™. The presence in an experimental 
material of any anisotropy which is not recognized by an appropriate modification 
of H*" can obviously cause an incorrect value of g to be obtained from the experi- 
ment. 


Ferromagnetic Resonance and the Internal Field 975 


§5. THE GENERAL RESONANCE CONDITION 
: The ferromagnetic resonance condition may be obtained from the solution of 
} the equation of motion (equation (5)) when the explicit form of the internal field 
components is known. The equation of motion may be solved in any convenient 
| coordinate system; the most convenient is often the rectangular system 
( (X,, X,, X,), defined by coincident ellipsoid and single-crystal principal axes. 
For the present work it is useful to introduce another rectangular coordinate 
} system (X,’, X,’, X,') rotated with respect to the (X,, X,, X;) axis with the static 


| 
| 
| 
| 


| the X,’ axis. By solving the equation of motion in this primed system, general 
} resonance conditions applicable to arbitrary orientations of the applied fields 
} with respect to ellipsoid principal axes, single crystal axes, or an applied or 
intrinsic stress system may be obtained immediately. 

In order to solve the equation of motion in the primed coordinate system, 
the internal field vector must be resolved in the primed system. The internal 
field contributions given in equations (7), (9), (10), and (11) involve magnetization 
» vector components referred to the unprimed system; the transformation of 
] the internal field expressions and the magnetization vector components to the 
| primed system is straightforward and is accomplished by introducing the nine 
| direction cosines y,; connecting an X;’ axis to an X; axis. One finds that the 
/ internal field components in the primed system due to shape effects, crystalline 
} anisotropy, or strain may each be written in the general form 


(H’,),/ x, y Ne ily» Sad Dene (12) 


) where the index ‘a’ denotes anisotropy due to shape, single-crystal structure, or 
i strain. The N*,, matrix elements are effective demagnetization constants due 
to anisotropy referred to the primed system; they will, in general, involve the 
2 © 

The equation of motion for arbitrary N*,, values may now be solved using 
equation (12). The components of the over-all internal field, including applied 


fields, are of the form CEE Pe emg BON alee ie (13) 
k 


mhete 1, =h, H,’=0, H,’=H,. Since the static field-, is taken as far greater 
/ than any static anisotropy field, the static magnetization vector will lie very nearly 
along the direction of the static applied field. In this case, the components of 
mene Magnetization vector are: J,’=7,' I,’=i,', I,'=1,. The solution of the 
} equation of motion with the above form of H” was first carried out by Kittel 
| (1947, 1948). In Kittel’s analysis the [.NV*,,] matrix was diagonal and the 
; following expression for H*" was obtained: 


H°8 = {[ Ho" + (N22 — N*53)Io]l Ho" + (N*11— N*s5)Zo]}?, «++ (14) 


+ where H,,' is the resonant value of the static field. 

| It turns out that the same result is an excellent approximation when the 
| [N*,,| matrix is not diagonal, provided that Hy is large enough to make J, lie along 
| its direction. In the general case of arbitrary orientation between applied 
fields and crystal or ellipsoid axes, the off-diagonal elements of [N*,,] are not, 
| in fact, zero. Nevertheless, the same resonance equation applies very closely. 
| In order to compute the resonance shifts due to the various different sources of 
| anisotropy, it is thus only necessary to calculate the diagonal elements of the 
' corresponding [NV*,,] matrices. When two or more sources of anisotropy are 
| 63-2 
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present in a specimen to be used in a resonance experiment, the matrix elements 
occurring in the resonance condition are the sums of the matrix elements arising 
from the individual anisotropies, as can readily be seen from the ne of 
equation (5) in H' and of equation (12) i in N*, 

In the Appendix general expressions for the (N? 11— N%33) and (N? oN oat 
factors occurring in the resonance condition are given for shape, crystalline 
anisotropy and strain anisotropy contributions. Special cases of these expressions 
for situations which have arisen or might be important in resonance experiments 
are considered in the following section. 


§6. ANISOTROPY CONTRIBUTIONS TO THE RESONANCE 
CONDITION 
(i) Shape 

Ferromagnetic resonance experiments have thus far been carried out with 
spherical magnetic specimens of high resistivity (Hewitt 1948, Bickford 1950, 
Yager et al. 1950) and with either high resistivity ferritic or conducting metallic 
disc-shaped specimens (Griffiths 1951). In these cases the applied magnetic 
fields are always taken along principal ellipsoid axes (considering a thin disc as 
an oblate spheroid) and the general treatment in the Appendix is unnecessary} 
the results given in equations (23) are only required when field directions and 
principal axes do not coincide. ‘The shape effects derived below apply both 
to polycrystalline and to single crystal materials. 

The effective fields for the above specimen shapes are easily obtained by 
substituting the shape demagnetization constants N*, directly for the N%, 
elements occurring in equation (14). For a non-conducting isotropic sphere 
N*, = N*,= N*,=47/3 and H,¢ reduces to Hj". ‘There is no shape anisotropy 
here and hence no resonance shift. For a non-conducting thin disc the demagneti- 
zation constants for fields lying in the plane of the disc are approximately 7?/m, 
while that for a field normal to the disc is (47 —27?/m), where m is the ratio of 
disc diameter to thickness (Osborn 1945, Stoner 1945). There are two field 
orientations of experimental interest, both applied fields lying in the plane of 
the disc, and static field perpendicular to the disc. ‘The demagnetization constants 
in these two cases are, respectively, N,=N,=7?/m, N,=47—27?/m, and. 
N,=N,=7?/m, N3=42—2n?/m. ‘The corresponding effective fields are 


e Tt 37? Tr he e€ 3 
He = | (2 SF 0) He sek | Hi (47 =) | 


(non-conducting discs). ...... (15) 


When the resistivity of a specimen is not so high that it may be considered | 
as virtually non-conducting, the skin effects begin to become important. Because: 
of this effect, the magnetic fields in the specimen are inhomogeneous since there: 
is a decay of the amplitudes of the oscillating components of the field as they? 
penetrate into the disc. In this case the expression for the magnetic potential! 
energy should contain terms arising from and depending upon the distribution of | 
V.I, which would cause the simple expression for E, given in §2 to be incorrect. 
A good approximation to the correct value of the effective field in this case may| 
be obtained, however, by remembering that at the high frequencies employed’ 
the skin depth is always considerably less than the diameter or thickness of the 
disc. ‘Then the demagnetization constants in the direction of oscillating field 
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{components have approximately the values they would have if the disc were of 
| infinite diameter. Therefore, for parallel fields, N,=0, N,=47, N3=7?/m, 
Hwhile for perpendicular static fields, N,=0, N,=0, N,=417—27?/m. The 
j effective field expressions are 


* 7 2 1/2 . , Qn 
[20 0 Se): [in (w- 2) 


(conducting discs). ...... (16) 


«Experiments designed to compare Sphere With gais. (parallel field orientation) 
for a polycrystalline ferrite of high resistivity were performed by Hewitt (1948), 
} who obtained reasonably good agreement between the g-values, indicating that 
{ the resonance formulae for these cases are at least good approximations. Further, 
} an experiment comparing g, and g, for a supermalloy disc specimen was carried 
J out by Kittel, Yager and Merritt (1949) who, found a difference between g, and 
}g, of 25%. Using the above expression for H,, which is slightly more accurate 
than that used by Kittel et al., one obtains a difference of approximately 1:1%; 
this can probably be ascribed to inexact knowledge of Jp. 


(11) Exchange 

In order to evaluate the effect of exchange anisotropy on the resonance condition 
it is necessary to add the exchange internal field contribution given in equation (8) 
to the equation of motion and solve it with this added term. If the magnetic 
material is non-conducting, however, the contribution will be zero since V°I 
7 will then be zero and Ix N,I is identically zero. Because of the skin effect 
present when the material is conducting, V*I is not zero and becomes V%1. 
Exact solution of the equation of motion with such a term taken in conjunction 
' with Maxwell’s equations is a lengthy complicated process since the material 
becomes triply refracting to microwaves. ‘The results of such a solution by the 
} author (Macdonald 1950) concur with those of a perturbation treatment by 
Kittel and Herring (1950) in indicating no appreciable resonance shift from this 
cause for metals at room temperature. 


(ii) Magnetocrystalline Anisotropy 


Thus far no resonance experiments have been reported on hexagonal-close- 
| packed single crystals. Experiments to determine the anisotropy constants K,' 
| and K,’ might be carried out with the static magnetic field always in a plane 
| perpendicular to the principal crystal axis or, alternatively, with the static field 
lying in a plane containing the crystal axis. If @ is the angle between the static 
| field and the crystal axis, the direction cosines in these two cases are, respectively, 
bY¥a=1, yu=Ya=0 and y,=sin#, ys,;=cos, yox3=0. The corresponding 
, contributions from equations (24) (see Appendix) to the effective field are 


2 
(N°, — N%s3)lo =0; (N° — N°s5)L = = ig (K,’ +2K,’) eoeevee (17) 
0 


and (N°, —N'%s3)Lo= ae cos 26 + = sin? 0(1 + 2cos 26) ; 
0 0 


2K,’ Koa: 
(N25 -- N33) = ie cos” 6 + A SLL 1S. = 9) la ca (18) 
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When 6 =0 the last two equations reduce to a result first obtained by Kittel (1948), 
who considered terms in K,’ only. For cobalt the factors 2K,'/I,) and K,’/Ip 
are approximately 6,000 and 700 gauss respectively, showing that crystalline 
anisotropy can affect the resonance condition greatly for this material. 

Two interesting experiments have been carried out on cubic single crystals. 
Kip and Arnold (1949) investigated the dependence of the resonant field for a 
silicon-iron single crystal disc upon the angle @ between the static field and 
the [001] crystal direction when the static field was constrained to lie in the 
(010) plane. In this case the non-zero direction cosines are y1,=y23=C¢089, 
13 = —Y31=8in8@ and y.2=1. The contributions to the effective field from 
equations (25) are 


KS 
(Nog aN slg = cos 40; (NS Nes) a= aI, > (3+.cos-48) + 5 ae 
0 


These results were first obtained by Kittel (1948, 1949a). Kip and Arnold 
found the expected dependence of H," upon 6 at constant applied frequency and 
obtained a good value of K, for silicon—iron as well as a g-value independent of @. 

Recently Yager, Galt, Merritt and Wood (1950) investigated the effect of 
rotating the static field applied to a single-crystal ferrite sphere in a (110) crystal 
plane. If @ is the angle between the applied field and the [001] direction, the 
non-zero direction cosines applying to this case are 


Yu=Vw2=V(1/2)ys3 = V (1/2) c08 8, 31 =¥92 = — V(1/2)yi3 = V/(1/2) sin @. 
The general ee (25) then reduces to 


CNS SIN Te sin? @ — 3 sin? 20) + igi sin’ 6(1 — 60 cos? 6+ 65 cos? @) ; 


(No — IN a3) lo= 2K, (1 —2 sin? 6— Sie 28) - 
TR 8 


1S 3 
Tél, (7 —5 cos? 0). 

The first-order terms in these equations have been reported earlier by Bickford 
(1950). Yager et al. found the 0-dependence indicated by the above equations 


and an excellent value of K, for the ferrite investigated: the K, terms were 
apparently negligible for this ferrite. 


(iv) Strain 

From the general equations (26) given in the Appendix it is relatively easy 
to compute the strain contribution to the resonance condition when the applied 
fields lie in a (010) or (110) crystal plane and none of the stress tensor elements is 
zero. Here, however, we shall consider some even simpler situations of more 
experimental interest. 

In the case of uniform pressure P applied to a ferromagnetic ellipsoid, the 
stress tensor elements are given by P,,= — Poi;, and there is no contribution 
from the general equations (26) because there is no strain anisotropy. This is 
not the case, however, for plane stress even in the simplest and experimentally 
most interesting case of a purely radial stress T, (Ty is positive for tension, negative 
for compression) applied at the circumference of a circular disc. Here the only 
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| non-zero stress tensor elements are (for a disc with its normal in the [010] direction) 
' P,,;=P33=T. On substituting these results into the relations (26), and simpli- 
/ fying, one obtains 


| Fr ; , aT: 
(MN oon aN ay Om aay (y12" Ts 32°) } (N oe Es N*53)1y = Oe sa ° 


(Y22" — — 39”). 


In an experiment to measure the effect of applied stress of this character, produced, 
for example, by means of a draw-strap around the disc, the applied magnetic 
fields would be taken either in the (010) crystal plane (the plane of the disc) or 
with the static field perpendicular to the disc. For these two conditions the above 
equations reduce to 


Parallelihields? (N“@,— Ns.) 1p =05) (Ns. — Nhs) Lp = 3Ar00 Lol Los 


Perpendicular static field: (N*,, — N**3)Ip =(N**o. — N53) 1p = =HMrooTo : 


In neither case is the contribution to the effective field zero, although rotation 
of the fields around the normal to the disc has no effect (crystalline anisotropy 
neglected) since there is no strain anisotropy in the plane of the disc. These 
non-zero results are somewhat unexpected and are discussed later in connection 
with actual experiments. 

Another simple stress system of experimental interest consists of a unidirec- 
tional stress T applied in the [f;, 8,, 83] crystallographic direction. Here the 
stress elements are P,,;=8,6,T. In the case where the applied fields lie in the 
plane of a single-crystal disc cut so that this plane coincides with the (010) 
crystallographic plane, and where the unidirectional stress is also applied in this 
plane, both the magnetostriction constants Ajo) and ,,,; may be obtained by 
_ measuring the change of H,' as the fields are rotated with respect to the crystal 
axes, keeping the stress direction fixed with respect to these axes or, alterna- 
tively, keeping the stress direction fixed and shifting the stress direction with 
respect to the crystal axes. No experiments of this type have been reported as yet. 

The above strain contributions have all applied to cubic ferromagnetic single 
crystals. Ferromagnetic resonance experiments are often made, however, on 
polycrystalline materials in which the individual crystallites have approximately 
random orientation. ‘The extension of the strain formulae to isotropic poly- 
crystalline material may be made in two different ways. In the first method 
the complete resonance condition, including self-energy, crystalline anisotropy, 
and strain contributions, must be averaged over all possible directions of the 
crystallite axes, preferably taking into account absorption line broadening due 
to relaxation processes, etc. Since this would be extremely complicated to carry 
out, even if crystalline shape and size distributions were known, it is simpler to 
make the approximation of isotropic magnetostriction usual in dealing with 
polycrystalline aggregates of ferromagnetic crystallites. This approximation 
consists of taking Ayo) =Ayy, =A. The effective strain demagnetization constants 
in equation (26) then simplify to give 


ri 3A , y RI 3 3X ’ (} 
(N**,, — N*53)I9 = Tt. (Ps3 ane ys (N** 49 — N**55)Io = 7 (Pss = Poo ), 
0 sees 0 
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where the primed stress elements are referred now to the primed (field) coordinate 
system. 

When isotropic stress Ty is applied in the plane of a polycrystalline disc, the 
effective field contributions are given by equations (22) with A,o9 replaced by A. 
The sign of these factors depends both upon the sign of A and upon whether 
the applied stress is tensile or compressive. ‘The change in the resonance field 
caused by these factors may be quite appreciable; for example, Kittel (1949 b) 
gives for 3AT,/I, the values 600, 600, 4,000 gauss for iron, cobalt and nickel, 
respectively, when TJ) is taken as the breaking stress of the material. 

No resonance experiments deliberately undertaken to observe stress effects 
on single crystals or polycrystalline materials have been published as yet. 
However, the resonant field dependence on applied plane stress has been 
verified by Macdonald (1950, 1951 .a,b)*, who investigated resonance shifts 
for thin nickel films evaporated on mica. Here, applied isotropic plane stress 
(independently determined from magnetometer measurements) was found 
to be produced by differential contraction between mica and nickel on cooling 
after evaporation. Such stress was dependent on film thickness and produced 
an apparent dependence of the g-factor on film thickness. No dependence of 
the g-factor on applied stress was observed, however, when the stress-corrected 
resonance condition was used. In the course of making resonance measurements 
from room temperature to the Curie point on a nickel disc silver-soldered to a 
copper block, Standley (unpublished) observed a decrease in the apparent 
g-value with increasing temperature. Calculations by the author (Macdonald 
1950) have shown that such decrease was only apparent, being caused by the 
temperature-dependent plane tension produced in the nickel disc by the greater 
expansion on heating of copper as compared with nickel. It may be postulated 
that a similar small decrease in g with increasing temperature for nickel found by 
Bloembergen (1950) and reported to be within experimental error was caused by 
the same phenomenon. Further, the 4°% decrease in g with decreasing temper- 
ature over the range from room temperature to — 153° c. found by Bickford (1950) 
for a Fe,O, single-crystal disc may be susceptible to the same general explanation, 
although the magnetostriction constant for magnetite may depend significantly 
upon temperature in this region. ‘The dependence of the resonant field upon 
unidirectional stress applied to a polycrystalline specimen has not been verified 
in detail experimentally so far, although Macdonald (1950) observed some 
directional dependence of H," upon the angle between the applied static field in 
the plane of an evaporated nickel disc and the coplanar direction of a magnetic 
field applied during annealing of the specimen. 


§7. CONCLUSION 


In all the resonance experiments discussed in §6 the g-values obtained were 
substantially greater than 2-00, the value to be expected for free electron spins. 
In fact the great majority of all g-values measured so far have significantly 
exceeded the free electron value, usually lying in the range between 2:1 and 2:3. 
One of the objects of the present work was to investigate fully the modifications 
in the resonance condition caused by internal fields arising from various types of 
anisotropy, with the hope that the results might possibly explain some or all of 
the anomalous g-values that have been reported. 


* See also Griffiths (1951). 
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The bearing of the present results on the g-value problem may be summarized 
| as follows: 

} (i) The use of the effective field expressions given in (16) for a conducting 
| disc decreases the experimental discrepancy between g, and g, found by Kittel, 
) Yager and Merritt to about 1°%. Since this is within the limits of experimental 
i error, g is independent of field orientation. Nevertheless, the value of g found for 
| supermalloy still remains about 2-19. 

(ii) The work of Kittel and Herring (1950) and Macdonald (1950) shows that 
| the high g-values found for metals cannot be ascribed to exchange force effects ; 
# such an explanation would not, in any event, apply to non-conducting materials. 

(iii) The results of experiments on anisotropic single crystals are in excellent 
agreement with theory. There appears to be no g-dependence upon special 
crystalline directions, but the experimentally measured g-values are anomalously 
# large. 

(iv) The g-factor has been found independent of applied stress, and the use 
© of the correction terms to the effective field found when isotropic plane stress is 
@ present in a disc-like ferromagnetic specimen has satisfactorily explained an 
} apparent dependence of g upon evaporated film thickness and upon temperature 
) for nickel discs. Stress effects may also be pertinent to the complete interpretation 
) of the g-value temperature dependence found by Bickford (1950) for magnetite 
| at low temperatures. 

| (v) The internal field results of the present paper, taken in conjunction with 
the experimental results of Macdonald on stressed and stress-free nickel films, 
would seem to render inapplicable Birks’ (1948) suggestion that high g-values are 
due to internal anisotropy or strain fields. 

There has been much speculation as to the reason for the anomalously large 
g-values. A critical survey of the reasons suggested for the effect has been given 
by Kittel (1949 a), see also Van Vleck (1950). None of the explanations seems able 
' to account for the phenomenon in detail, although it is likely that the effect is 
connected with spin-orbit coupling as postulated by Polder (1949) and by Kittel 
(1949a). There is as yet no detailed theory of the phenomenon, and it remains 
. essentially unexplained. 


APPENDIX 


EFFECTIVE DEMAGNETIZATION CONSTANTS DUE TO ANISOTROPY 

The problem considered here is the transformation of equations (7), (9), 
(10) and (11) for anisotropy internal-field contributions (referred to the fixed 
coordinate system (X,, X2, X;)) to the rotated coordinate system (X,’, X.', X3') 
in order to obtain expressions of the form of equation (12) for the internal field 
contributions in the primed system. ‘The N*,, demagnetization elements due 
to anisotropy defined by the transformation may then be used in equation (14) 
to give the general resonance condition for any orientation between primed and 
unprimed axes. 

(i) Shape (cf. equation (7)) 

Since equation (7) is linear in J,, the N*; shape demagnetization constants 
may be considered the diagonal elements of the demagnetization matrix referred 
to the unprimed principal-axis coordinate system. A coordinate tranformation 
of this matrix to the primed system then gives the desired result : 


N ix = a ViY nN) 
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Therefore the diagonal-element differences appearing in the resonance condition 
are 


(Min = N*33) hy bats = V3 1N*, } (N* 29 a N*33) ac [yer ss Yar IN*,. serene (23) 


(ii) Magnetocrystalline Anisotropy: Hexagonal Crystal (cf. equation (9)) 


The above procedure for obtaining N*,, elements is not applicable to this 
case because Hi, is here not a linear function of the J; components. Instead, 
the H', vector must be resolved directly in the primed system and the unprimed 
I, components occurring in the resulting (H’.),’ field components expressed in 
terms of the I,’ components. A linear dependence of the (H'.),’/ components 
upon J,', as in equation (12), is then obtained by dropping terms containing 
powers of 7,’ greater than unity. Such neglect is justified as long as J) >2,'; 
this condition holds in practice. One obtains 


P 2 ; ' } 2 
(N°,—-N 33) aT T,? [(Ky +2K, Mysx" a Yu") +2K, (3y11° a ¥31°)¥31° 15 | 


(24) 
. c 2 i} , 2 2 / ¢ 
(N%2 — N%33) = I, [(Ky! + 2K9')(y31" — yar) + 2K 9’ (3y21 — ¥s1")ya1° I] 
where I has been replaced by J. 


(iii) Magnetocrystalline Anisotropy: Cubic Crystal (cf. equation (10)) 


To simplify the succeeding formulae, the following functions of the direction 
cosines are defined : 


if SS Say 2 2. h —_ 
LL 7 Yjl Vet > l,m, n— YY m2Y n3> 
a m,n lacy nl Ay, m,n Be 2h, nt “te aes lc 


The desired demagnetization element differences obtained as in (ii) are then 
c Cc 2 
(N°, — N33) = Ie [Ka(fss — 3/13) — Ko{gs, 3,3— (81,3,3 + 83,3,1+83,1,3)}) | 


25) 
Cc iC 2 ( 
(N°x2 — N33) = 73 Aisa — 3f25) — Kote, 2,3 — Bs, 3,3+83,3,2+83,2,3)}I- 


The terms in K, are in agreement with results obtained by Van Vleck (1950) 
from a quantum-mechanical treatment. General expressions for the second- 
order terms have not been given previously. As Van Vleck points out, it is 
noteworthy that the classical and quantum-mechanical treatments give exactly 


the same results, at least to first order. 
(iv) Strain (cf. equation (11)) 
Define the following quantity in terms of the stress tensor elements referred 
to the unprimed crystal-axis coordinate system: 
3 
Six = P [(Axoo — A111) 95x + Ayu] Pipe 

Then the expression for the internal field contribution may be written 

(3), ie 2 Siied is 


‘Transforming the matrix [.S,,] to the primed system, one obtains 
sti 2* 
N ik oe - x Viv km Som 
m 
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Therefore the general expressions for the strain demagnetization element 
differences are 
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ABSTRACT. The theories of the mean free path of electrons in various types of materials, 
under various conditions, are discussed. The conductivity, thermo-electric power, Hall 
coefficient, fractional change of conductivity in a magnetic field, and the Nernst, Ettinghausen 
and Righi-Leduc coefficients are all calculated upon the Lorentz—Sommerfeld theory for the 
various mean free path theories. The variations are discussed with particular respect to the 
changes of the reduced chemical potential 7*. The magnitude and sign of the coefficients. 
are found to depend upon the mean free path theory assumed, that is, upon the type of 
material. In particular it is shown that the Ettinghausen and Nernst coefficients should be 
negligible in the majority of semiconductors with ionic lattices. Some of the assumptions 
usually made in semiconductor theory are shown to lead to erroneous conclusions. Experi-- 
mental evidence is discussed and is found in general to support the theory given, both in 
impurity metals and semiconductors, except perhaps in the case of the thermo-clectric 
power. No experimental results in non-metals are known for the thermo- and galvano- 
magnetic effects apart from the Hall coefficient, and it appears likely that such measurements 
will produce interesting information. 


§1. INTRODUCTION 


ROM the earliest days of the free-electron theory of conduction in solids the 
electrical properties and their temperature variations have been expressed 
in terms of two functions, namely the distribution function of the electrons 

in the solid, and the mean free path of the electrons. It is now accepted that the 
appropriate distribution function is the Fermi—Dirac, which under some conditions. 
reduces to the simpler Maxwell—Boltzmann law. For the mean free path no. 
universally applicable function has been obtained hitherto. On the contrary, 
the mean free path function differs according to the character of the solid (whether 
it possesses an atomic or an ionic lattice, for example), and according to the 
temperature and electron density. In the past the appropriate mean free path 
function has been by no means invariably used in discussions of experimental 
results, and assumptions of the mathematically simplest function have not 
infrequently resulted in an apparent failure to explain the experiments. More- 
over this same avoidance of the correct procedure has led to theoretical results 
which, in the main, are approximations for the limiting cases of metals, or of 
materials containing so few conduction electrons that the Maxwell—Boltzmann 
statistics apply, leaving a wide gap into which many materials necessarily fall. 
The objects of this paper are to lay down a more general procedure employing the 
most appropriate mean free path function in all cases, and enabling intermediate 
substances to be considered. As a secondary objective, fuller consideration is. 
given to the less familiar galvanomagnetic and thermomagnetic effects, from 
which valuable information is obtainable, but which have been much more 
neglected by the experimentalists. Therefore, in addition to conductivity, 
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thermo-electric and Hall effects as normally considered, we shall consider also 
-the Nernst, Ettinghausen and Righi-Leduc effects, and finally the effect of change 
of resistance in a magnetic field. Since the theory of metals is already well 
understood, non-metals only will be considered here. Finally a comparison 
is made between the results predicted by the theory and those observed expe-i- 
mentally. In order to do this, simple approximate formulae are obtained for 
the variation of the electrical properties with temperature, which are applicable 
over limited temperature ranges. ‘The experimental results for many materials 
are found to be in good accord with the approximate formulae, and it seems 
likely that many more would be explained if the approximations were no longer 
used, and the full theoretical expressions compared with experiment. 


§2. THE DISTRIBUTION FUNCTION AND DEGENERACY 


_ It is necessary first to explain the manner in which the Fermi—Dirac distri- 
bution function will be employed. According to this function the concentration 
of free electrons 7; is given by 


ng=4ah-8(2m* RT)? Fyo(n*) ns ee (1) 


where TJ is the absolute temperature, A and R are Planck’s and Boltzmann’s 
constants, m* is the effective mass of the electron and Fy). (7*) belongs to a group 


of functions defined by 6g Fs 
F ( *) = (GS ee (2) 
m( CER RESTO egal 5 oes lee 

In this definition 7 = E/RT the ‘reduced’ energy of an electron, the value of E 
being the energy of the electron relative to the bottom of the conduction band. 
n* is the reduced chemical potential, and its value is calculated by integrating 
the number of electrons in each energy state, and equating to the total number 
of electrons. ‘The equation (1) may be regarded as expressing 7* in terms of 
the concentration m;. In several important papers the reduced chemical potential 
7* has already been used very effectively as a parameter indicative of the degree 
of degeneracy, and it will be so used here. It will be observed that * is inde- 
pendent of the mean free path but is dependent on the distribution of energy 
levels, and on the distribution of electrons in the energy bands. When the 
variation of 7* (and thus the degeneracy) with temperature is considered, non- 
metallic solids divide naturally into two groups: semiconductors in which the 
number of free electrons is zero at absolute temperature but increases as temper- 
ature rises, and impurity metals where the number of free electrons is a constant 
independent of temperature. In the first case it appears that 7* is large and 
negative at low temperatures, rising to a maximum as the temperature rises, 
and then decreasing again at still higher temperatures. In the second case 7* 
is large and positive at low temperatures, and decreases as the temperature rises. 
These two groups consist usually of substances in which the electrons are due to 
a non-stoichiometric excess of one component of the lattice, or to impurities. 

In more detail, for impurity metals, since m, is constant, and equal to the 
concentration of impurity centres, equation (1) expresses the variation of 7* with 
temperature. At low temperatures, or more precisely when m, h® >47 (2m*RT)3 


and 7* >0 we have  (3m\8 : 
BaSee eT aenistit)) 
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showing that 7* is inversely proportional to the temperature. At high tempera- 


ture, i.e. when 7, h?< 47 (2m*kT)3 


ie at n,h® 
n *—In | soa | G5 0-010 (4) 


showing that * varies but slowly with temperature, becoming large and negative 
as the temperature rises. At intermediate temperatures, however, the full 
expression (1) must generally be used. ; 

In a semiconductor, however, the problem is more elaborate, as the number of 
free electrons varies with the temperature. Formula (1) still holds for the relation. 
between m, and 7*. 

The equation for * is then : 


ny, =4rh-3(2m*kT)?? [1 + exp(n’+7*)] Fye(n*)  —.. . (5) 
where 7’ =E’/RT, and m, is the concentration of impurity centres with energy 


levels a distance E’ below the conduction band. 
For low temperatures (E’ > RT), then the following approximation will hold 


good: Ae 
Si 2 APU geil ey | ee ea ae 
Lees 3 hag In | Soar | Sn ee (6) 


which is such that 7* is inversely proportional to the temperature, and is large and 
negative. ‘This is the case of the classical semiconductor. 

As the temperature rises and RT becomes of the order of £” and greater, then 
the above approximation breaks down. ‘The value of 7* increases from its large 
negative value inversely proportional to T whilst the approximation (equation (6)) 
holds good, then as it breaks down the value of »* fails to increase so rapidly, 
reaches a maximum and begins to decrease again. Eventually, at sufficiently 
high temperatures, the situation will resemble that of the impurity metal at high 
temperatures since all the available free electrons will be free, and the value of 7* 
will decrease logarithmically with temperature as shown in equation (4). 


§3. REVIEW OF MEAN FREE PATH THEORIES 
In metals and substances with atomic lattices it has been shown (Wilson 1936), 
that if the energy FE of an electron above the bottom of the conduction band is 
small, i.e. H<RT, then the mean free path is, to a first approximation, inversely 
proportional to the temperature, and independent of the energy E: 


ITS) ES) (7) 


where /, depends upon the constants of the lattice. Shifrin (1944) has shown that 
the presence of degeneracy which was assumed absent in Wilson’s derivation, does 
not affect formula (7) to the first approximation. When the energy E£ of anelectron 
is high, then its mean free path is given by (Wilson 1936) 

L= bh Ta® | >| OS VE eee (8) 
where 7=E/kRT, and is the reduced energy. This is the case which is applied 
to metals, and has of course been dealt with elsewhere in detail. "The assumption 
is made that formula (8) may be applied to all the ‘ conduction’ electrons with 
reasonable approximation. ‘The variation as predicted by formula (8) will not be 
discussed in this paper. 

In an ionic lattice the mean free path is dependent upon whether the electron 
energy E above the bottom of the conduction band is greater or less than hv, the 
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| energy of the characteristic vibrations of the lattice. For energies less than hv, 


4 Mott and Frohlich (1939) and Davydov and Shmushkevitch (1940) have shown 


| eet Delia ees MCS Dayne eo ee eee (9) 

! This formula will be applicable to the vast majority of electrons present in an ionic 
} lattice if the temperature and concentration of free electrons are low, so that few 
electrons have high energies. Shifrin (1944) has suggested that this formula is 
@ only true when the concentration of free electrons is so small that classical statistics 
| may be applied, and that the presence of degeneracy modifies the formula to the 


following: 
ollowing FATR Tut eh Fe R(t fy oes (10) 


where f=[1+exp (7 —7*)]"1 is the Fermi—Dirac function; and 7* = E*/RT is the 
reduced chemical potential. We shall see that this modification has very little 
effect upon the electrical properties of the material so long as 7*<—4. This is 
usually the case in ‘ classical’ semiconductors with ionic lattices, and it is to this 
group of materials to which equation (9) may be considered applicable. Shifrin’s 
modified, and more complicated, formula (10) whilst applicable to these materials 
need only be used when the degeneracy becomes appreciable, or 7*> —4. The 
modifications introduced by formula (10) are particularly important in connection 
_ with ionic latticed ‘impurity metals’. ‘These are materials which possess a 
constant concentration of free electrons in an ‘ impurity band’ of permitted energy 
levels, due to the presence of foreign impurities or an excess of one constituent 
above stoichiometric proportion. At low temperatures the degeneracy in such 
_ materials is high, as in metals, and thus the factor (1 —/f) has an appreciable effect. 
The application of formula (10) to ionic latticed impurity metals is restricted to 
- cases where the concentration of free electrons is such that few have energies 
greater than hv. 5 

For electron energies E greater than hv the mean free path has been shown to 
be (Mott and Frohlich 1939): 


PEE Ree Uhh oe SOW at ethic (11) 


where /, is a constant for temperatures greater than the Debye temperature 
(Wright 1951). ‘This equation may be applied to all ionic latticed semiconductors 
and impurity metals at high temperatures. 

In certain cases equation (11) may be applied to alltemperatures. Ifthereisa 
high concentration of free electrons in the impurity band of an impurity metal, 
then at low temperatures the Fermi energy E* will be large, viz. E*>hv. Under 
these conditions all the effective ‘ conduction’ electrons, i.e. those near the Fermi 
level, will have energies E> hy, and thus formula (11) may be used. Hence the 
use of equation (11) for all temperatures, for all the electrons, will result in a fair 
approximation. 

Finally it may be pointed out that in the original theories of semiconductors the 
mean free path was assumed to be a constant, independent of temperature and 


ey: ee (12) 


Many conclusions derived from this assumption are extremely erroneous, as will 
be seen in the following discussion. 

Summing up, we see that all the above expressions for the mean free path may 
be expressed in the form /=/, f(y) where J, is independent of the energy F (i.e. 7) 
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but is, in general, dependent upon the temperature, and f(7) is a function of 7 only, | 
The above results are collected in the Table. 

In the following discussion the general form /= 1, f(y) will be used, both J, ate 
f(y) depending upon the actual material and conditions under consideration. | 


Limits of 


Material SRC i Equation lo f() 
Classical case — (12) lls 1 
Atomic latticed materials Non-metals (7) ee 1 
Ionic latticed materials Very high T (11) ly n 
Ionic latticed semiconductors Low T, low r¢ (9) 1,(RT)1/2che/kT yi 


Ionic latticed impurity metals Low T, low n¢ (10) I(RT)2ehrkiP p12 eae 
Ionic latticed semiconductors Low T, high mg (10) I(RT)!eh/RL 12/1 —f) 
Ionic latticed impurity metals High 7¢ (11) Un n 


$4) “THE BLECTRRICAL BRORPERTIE'S 
We are now ina position to write down expressions for the electrical properties 
calculated from the Lorentz—Sommerfeld theory in a generalized form suited to a 
discussion of the differences made by the different mean free path functions. 
These expressions are : 


Electrical conductivity =A loamte2R7 1)/sh*> (en eee (13) 
Thermo-electric power O=-BRie- — « © | i) peewee (14) 
Hall coefficient R==Cin-e2 9. "yaa eee ee (15) 
Fractional charge of o ina 

magnetic field H dsjo= DA e722 kl (16) 
Nernst coefficient O=E Ri (2mtkiI 2 i eee (17) 
Ettinghausen coefficient P=3F/2nk “(7.0 Seer (18) 
Righi-Leduc coefficient S==Geli2QmtrT 2a) Oe eee (19) 
iWihere the? letters’ A,B... : G represent functions of the Fermi—Dirac 


functions F’,({y*). These functions differ according to which mean free path 
theory is used, and the various values are given in Appendix I. It must be 
remembered that J) is in general dependent upon the temperature. 

The Nernst coefficient Q is defined by the equation F,, =QH, dT/dx where F, 
is the electric field in the y direction, and the rest of the notation is as usual. The 
Ettinghausen coefficient P is defined similarly by dT/dy = PH,I,, where I,, is the 
current density in the x direction. ‘The Righi-Leduc coefficient S is defined by 
dT/dy=SH,dT/dx. 

The majority of the expressions for the constants A, B and C have been derived 
before, but the expressions for D, E, F and G given in Appendix J are original, since 
these effects seem only to have been considered in the limiting cases of metals and 
the ‘classical’ case. ‘The variations of the expressions A to G with y* are shown in 
Figures 1 to 7. In each Figure there are four graphs: graph I shows the variation 
with 7* of the appropriate constant if the mean free path is assumed to vary as 
shown by equation (12) or equation (7); graph II if the mean free path is given by 
equation (9), graph III if by equation (10) and graph IV if by equation (11). The 
variations are shown with respect to the variations of 7*, because 7* is independent 
of the lattice structure and therefore of the mean free path. The manner of the 


we 
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Figure 4. 


Ettinghausen Coefficient 
P= 3F/2n,k 
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Nernst Coefficient 
Q=EKl/(2mkT )* 
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' Figures 1-7. Variations with * of the functions 4—G of equations (13)-(19), as given by Appendix I. 
PROC. PHYS. SOC. LXIV, II——A 64 
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variation of 7* with temperature depends entirely upon whether the substance is a 
semiconductor or an impurity metal. It should be remembered, however, that 
when the variation of the properties with temperature is considered, then the vari- 
ation of J, must be taken into consideration, particularly in the case of ionic latticed 
materials at low temperature, when the variation of J, with temperature is 
exponential. 

In the next section approximate formulae are developed to give the variation of 
the properties with temperature. Where these approximations are not valid, 
which may be the case in practice, then the accurate value predicted by theory 
must be deduced from equations (13) to (19), and the graphs for the functions 
A-G. 

The following important points may be noted from the graphs : 

(a) The variation of the thermo-electric power with y* is completely different 
in ionic latticed materials depending upon whether the mean free path follows 
equation (9) or (10), i.e. whether the mean free path is dependent upon the degree 
of degeneracy or not. If7* is large and negative then the two equations give rise 
to identical results, but as 7* becomes large and positive, then equation (9) predicts 
that the thermo-electric power approaches asymptotically a minimum value of 
© =0 (graph II, Figure 2) whilst equation (10) predicts that the minimum value 
which will be approached is © = — k/e (graph III, Figure 2). If equation (10) is 
correct then the lowest value that the thermo-electric power of an ionic latticed 
material at low temperature can have is 86 microvolts. In the majority of 
semiconductors, however, 7*< —2 and thus equation (9) and equation (10) give 
rise to identical results. The difference will be important, however, in the case 
of impurity metals at low temperature. 

(b) In ionic latticed materials with low degeneracy, and at low temperatures, 
both the Ettinghausen and Nernst coefficients are zero. Any presence of such 
effects in these materials will indicate that Shifrin’s modified formula (10) should be 
used, and also that the degeneracy is appreciable, since even equation (10) predicts 
negligible values if 7*< —2 (graph III, Figures 5 and 6). This 1s of practical 
importance in the measurement of the Hall constant in which precautions must be 
taken to eliminate any E.M.F.s due to the Ettinghausen transverse temperature 
gradient producing a thermo-electric E. M. F. in the Hall potential probes. In the 
majority of ionic latticed semiconductors 7* is usually large and negative, and the 
graphs show that effects are extremely small. No actual measurements of these 
effects in such materials is known, though precautions to eliminate them are often 
mentioned. Experiments are in progress to measure these effects in an impurity 
metal, CdO, where they should be appreciable. 

(c) The actual directions of the Ettinghausen coefficient P and the Nernst 
coefficient Q are dependent upon the mean free path theory adopted. It should be 
noticed that both P and QO are independent of the sign of charge carried by the 
current carrier, yet in atomic lattices (graph I, Figures 5 and 6) the effects are 
positive, whilst in ionic lattices (graphs II, III, IV of Figures 5 and 6) coefficients, 
if appreciable, are negative. Furthermore, it will be seen that the assumption 
that /) is a constant, which is often made in the simple classical theory of semi- 
conductors, leads to completely erroneous conclusions in respect of both the sign 
and the magnitude of the coefficients in the case of an ionic latticed material. 

(d) It will be noted that in materials which are subject to equation (11), such as 
CdO, the change of resistance in a magnetic field and the Righi-Leduc effect 
should both be comparatively large. 
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(e) The thermodynamical relationship that Q = K P/T (where K is the thermal 
conductivity of the electrons) is obeyed for all values of * and all the theories of 
the mean free path. The empirical relationship RO/PS = K/oT is correct when 
»* is large and the material is highly degenerate, as in a metal, where the relation- 
ship is found to be obeyed. 

(f) If the assumptions that are made in the classical theory, viz. 7* < 0 and 
J=I,, are taken, the above equations lead to the values for the electrical properties. 
»shown in Appendix II, which are the same as those which have been obtained 
} previously. There are, however, no substances to which such simplifying 
/assumptions can be considered applicable. It will be noted from Appendix II 
| that R = 37/8n,e, which is usually assumed for semiconductors, is deduced only if 
| the mean free path is assumed constant. In the majority of semiconductors with 
‘ionic lattices it will be seen that R = 1/n,e is more nearly correct. 


f §5. COMPARISON WITH EXPERIMENTAL RESULTS 
q . aye 
F (i) Conductivity 
In atomic latticed materials such as silicon and germanium several investigations 
have been made of the variation of conductivity with temperature. Consider such 
» a semiconductor with a small activation energy EZ’ and number of impurity centres 
_m. At all temperatures the mean free path is given by equation (7), and thus the 
Bexpression for o may be found as in Appendix I: 
| l6am*ke?l, 
G2 SSS 
3h 
where it should be remembered that /, is independent of the temperature. At low 
_temperatures 7* is large and negative, so that equation (6) may be used as an 
' approximation for 7*, giving 


Ee ex pit) ae? ORI Aone (20) 


oC 


_ 8m," e7), (ny enw 
1 Ul A ll aA 4 
| Thus In o is proportional to 1/T as is usually found in poorly conducting semi- 
‘ conductors. As the temperature increases, however, 7* increases beyond the 
| limits where equation (6) is a good approximation, and eventually reaches a 
| maximum, after which it decreases. Shifrin (1944) has shown by calculation that 
_ »* decreases in such a way that In(1 + exp 7*) decreases linearly with temperature. 
The conductivity should therefore follow equation (21) at low temperatures, 
reach a maximum, and then decrease proportionally as T increases. Precisely 
| this variation has been observed by Putley (1949) for Si and Get. At higher 
» temperatures still, intrinsic conduction occurs and masks the impurity conduction. 
Pearson and Bardeen (1949) obtained results for silicon very similar to those of 
' Putley’s. ‘They also showed that as the concentration of impurity increased, the 
material became an impurity metal, as was shown by the Hall coefficient (and thus 
the concentration of free electrons ,) becoming independent of the temperature. 
The conductivity in atomic latticed impurity metals is given by equation (20), 
but »* varies with the temperature in the manner described in §2. At low 
temperatures, when 7* >2, we may use equation (3) so that 


2 2/8 
| wore el if (02) 


SE SSS 
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. + Hutner, Rittner and du Pré (1950) have shown that the agreernent is good except below 
150° k. when the additional scattering due to the impurity centres (not considered here) must be 
considered to give good agreement. 
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At high temperatures, when 7 * < — 2, we may use equation (4) for 7* so that 
— 4h, 1 
°~ 3(Qam*ky2 T32 
Although some of Pearson and Bardeen’s specimens became impurity metals 
when the concentration of impurities was sufficiently high, neither of the above 
variations with temperature was compared with those found experimentally, 
though such a type of variation seems to be indicated by the graphs given in their 
paper. 
In ionic latticed semiconductors, the theory indicates that at low temperatures 
In o should be proportional to 1/7, as in the usual semiconductor theory : 


o =(2nkT)%4m*4n1? e],h-82 exp(2hv—E)/2RT —..... (24) 


Such a variation has been observed in the majority of poorly conducting semi- 
conductors, and was, at one time, considered to be a test of whether a material was 
a semiconductor. It should be noted that in equation (24) the variation of the 
mean free path with temperature has been allowed for, and the slope of the graph 
of Ino against 1/T should be (2hyv — E’)/2k, and not, as is often assumed — E’/2R. 
Usually, however, 2hv < E’. When 7* increases and approaches zero, then the 
approximation given by equation (6) breaks down, and equation (24) is no longer 
applicable. Such a change has been noticed by Busch and Labhardt (1946) in 
SiC. The graph of Ino against 1/7 is no longer astraight line. ‘They found that 
use of the accurate expression for the conductivity (see Appendix I) was in good 
agreement with the theory. When 7* increases and reaches its maximum, then at 
higher temperatures 7*, and thus Fy/.(7*) and m,, will no longer vary rapidly with 
temperature. Shifrin (1944) has shown that for many specimens m, should vary 
roughly proportionally to 7. In this case the variation of the conductivity will be 
mainly controlled by the exponential term in the mean free path, and so the 
conductivity should decrease exponentially as the temperature rises. No simple 
analytical expression can be obtained. Such a variation of conductivity with 
temperature has been observed in several materials, but is capable of a different 
explanation, and such results will be discussed later. 

In an ionic latticed impurity metal, we may consider first the case in which 
there are sufficient free electrons, such that equation (11) is valid at all temperatures. 
In this case Appendix I tells us that 
_ 167m*e?l,R 

Leoecrpsne 

At low temperatures, when 7 * >2, this approximates to 

318 e2hn,,48 Up 
oat “4m* kris Gp oe EP PPEOR I ie «Fe oho oherenn: 
The variation of conductivity with T predicted by this equation, has been found | 
to hold very well in the case of specimens of CdO (Wright 1951). It should be: 
noted, however, that at temperatures less than the Debye temperatures, /, is no } 
longer independent of the temperature, and the conductivity should approach a 
uniform maximum value as 7 decreases to zero. At high temperatures when| 
»* < —2, then equation (25) approximates to 

8e?1,ny 1 / 
3(2nm*k)1? T12 6 ORS, FO aw Le | ee a aterel ene) 


TEias) See (25) 
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The second type of ionic latticed impurity metal is that in which the concentra- 


j tion of free electrons is sufficiently small that the vast majority of free electrons 
¥ has energies less than Ay, and thus equation (10), for the mean free path, is valid. 
} It is not quite certain whether this type exists since the concentration of impurity 
* centres, and thus free electrons, must be large for the material to become 


an impurity metal. The conductivity of such materials, if they exist, is given at 
low temperatures when 7* >2 by 


Bit ig POs loa teittigion mw, antlirali ier. (28) 


aa D322 


The conductivity should therefore increase logarithmically as 1/7. At high 


temperatures the conductivity should vary as equation (27), as should all ionic 
latticed materials at very high temperatures. 

Interesting experimental results have been obtained using specimens of PbS 
by Devajatkova, Maslakovec and Sominskij (1941), with specimens of PbSe 
by Bauer (1940), with specimens of PbTe by Putley (private communication), 


» and with specimens of V,O3 by Foéx (1949), which have not been fully explained 


theoretically. For the curves given for the resistivity of PbS against temperature, 
there appears a definite indication of a variation in agreement with equation (26) 
asa highly conductiveimpurity metal. Shifrin has suggested that the conductivity 
is proportional to exp(hv/RT), and thus indicates either a highly conductive 


} semiconductor or an impurity metal as in equation (28). It is difficult from the 


graphs given in the paper to decide between the two variations. Unfortunately, 
no measurements of the Hall constant were made, which would have instantly 
distinguished between a semiconductor and an impurity metal. In PbSe, also, the 
conductivity appears to increase in proportion to exp (hv/RT) and is therefore 
following equation (28) as an impurity metal, or is acting as a highly conductive 
semiconductor. Similarly in V,O, the conductivity increases exponentially at 


low temperatures as expected of semiconductors, reaches a maximum, and then 


decreases such that a graph of Ino against 1/7 isa straight line. In PbTe the 
conductivity is also found to be proportional to exp (hv/RT), in this case for both 
excess and defect conduction, in the region of temperature where the Hall constant 
is practically independent of temperature. ‘The conductivity is therefore repre- 
sented by equation (28). In both V,O, and PbTe it seems likely that the material 
is a semiconductor, and the conductivity follows equation (28) wher. the majority 
of the impurity centres are ionized. 


(ii) Thermo-electric Power 
In general the results for the thermo-electric power have agreed less with theory 
than the other electrical properties investigated. ‘This may be due to the fact 
that the theory of the thermo-electric power contains many over-simplifying 


assumptions. 
In atomic latticed semiconductors the theory predicts that 
k[2F,(7*) | 
e=-5| z —n* ° eccece 29 
eL Fa) 7 aa 


When the temperature is low, then equation (6) is applicable, and F, = Fy, 


Let hae me) 
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This is the type of variation which was considered typical of a semiconductor, 
but there are but few atomic latticed semiconductors, and none seem to obey this 
type of equation. As the temperature rises so 7* increases and eventually reaches 
its maximum, and then starts to decrease again ; so the thermo-electric power should 
decrease in proportion to 1/T at low temperatures, and at higher temperatures 
reach a minimum, and finally increase in proportion to In J. Putley’s measure- 
ments (1948) on silicon and germanium found that the thermo-electric power 
increased rapidly with temperature, over part of its range roughly proportional to 
temperature, thus contradicting the above theory. 

Atomic latticed impurity metals should obey equation (29) and at low tempera- 
tures we may approximate for 7* by equation (3) when »* >2 


hen (8mr\ #8 2k | 
ae 


The thermo-electric power should therefore be proportional to the temperature. 
‘This is more in accord with results found for specimens of silicon and germanium 
by Putley, but these specimens, as far as conductivity is concerned, acted as semi- 
conductors. At high temperatures when »*<—2 the variation with tempera- 
ture should be the same as that predicted for semiconductors, i.e. proportional 
to In T. 

In ionic latticed semiconductors at low temperatures, the thermo-electric. 
power is the same as equation (30), except that the constant 2 in the parentheses 
should be changed to 5/2. ‘The thermo-electric power should therefore be 
inversely proportional to temperature. ‘This type of relation is found in many 
good insulating semiconductors. In NiO (Wright and Andrews 1949) the 
numerical values which may be deduced from the graph of © against 1/T were 
found to be in good relation with those deduced from both conductivity and Hall 
coefficient measurements. As the temperature rises, so 7* increases to its 
maximum, and then decreases; but for values of * in the region of zero then the 
only accurate formula is the full one 


R (5 Fs)(y*) 
O22 5- (3 pelt NY i Pty *) lh Geta ees 32 
€ \3 Fyp(y*) ‘ oS 


The thermo-electric power should therefore decrease to a minimum and then 
increase, and eventually, at sufficiently high temperatures, become proportional to 
Ingi: 

Impurity metals, with an ionic lattice and a high concentration of free electrons 
so that equation (11) applies to all temperatures, have at low temperatures when 
et ke Qn? (ayes 
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Such a variation has been found definitely to exist in specimens of CdO which 
acted as impurity metals. 

A similar type of variation has been observed by Busch, Lundt and Spodlin 
(1941) in SiC, but whether the material acted as an impurity metal in other respects 
isnotknown. At high temperatures, 7* < —2, then 


= R nh? 
o--* E —In (sone - 7m) | Ss eae (34) 
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so that as the temperature increases, the rate of increase of ® with temperature 


a) = 


decreases until it is proportional toln7. If the concentration of the free electrons 


_ in the impurity metal is low, we have at low temperatures 


_ Rk 87 \23 mk oe 
o-- [1+ (=) 1G +5 )r]. Sr alerasece (35) 


This formula indicates only a slow variation with temperature. At high tempera- 


a _tures such a material would follow equation (34). Considering the case of the 


experimental values obtained by Devajatokowa, Maslakovitz and Sominiskii 
(1941) on PbS, they obtained values which seem at low temperatures to be 
proportional to the temperature, but whether according to equation (35) or 


- equation (33), depends upon the specimen. Both types of variations do seem to 


be present. At high temperatures there are indications that © is proportional 
to InJ. These observations would indicate that the specimens of PbS used 
were impurity metals, and this agrees with some of the deductions based upon the 
variation of the conductivity with temperature. 


(iii) The Hall Coefficient 


The function C in equation (15) is seen in Figure 3 to vary between unity and 
1-33 for all types of materials and alltemperatures. Variations of C with tempera- 
ture is therefore very slow, and C may usually be taken as a constant, so that the 
variation of the Hall coefficient with temperature is purely dependent upon the 
variation of the concentration of free electrons. In semiconductors at low 

temperatures, then, InR should be proportional to 1/7, and this has been observed 


' in many semiconductors. As the temperature increases this simple variation 


should disappear and only the accurate formula suffices (Busch and Labhardt (1946) 
found this to be the case in SiC). At higher temperatures, R should decrease 
linearly with temperature, and finally approach a minimum constant value. In 
impurity metals the concentration of free electrons is independent of temperature; 
thus so is the Hall coefficient. Several substances have recently been found 
to illustrate this type of material, such as some specimens of Si, PbSe, PbTe and 
CdO. ‘The invariance of the Hall coefficient with temperature is a sure indication 
of an impurity metal. 


(iv) Change of Resistance in a Magnetic Field 


Very little experimental work has been done on this phenomenon. It will be 
seen from Figure 4 that except for material for which 7*>0 the variation of D 
with 7* is of little consequence. Thus in the case of semiconductors D may be 
considered a constant, and the variation with temperature will depend entirely 
upon the remainder of equation (16). _It will be seen that investigations of this 
effect will therefore give direct information upon the variation of the mean free 
path with temperature. Putley (private communication) has shown by observa- 
tions at one temperature on PbTe that do/o is proportional to H?, and that the 
fractional change of Hall coefficient with magnetic field dR/R is also proportional 
to H?. Recently Estermann and Foner (1950) have measured do/o for specimens 
of germanium, and their results are in reasonable accord with theory. ‘The only 
other results are those of Kapitza (1929) on much more impure specimens of 
germanium, which gave much smaller values for do/o. ‘This is to be expected, 
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since the degeneracy is much higher, and thus, as will be seen from graph I, 
Figure 4, the value of D should be much less. 

For impurity metals 7* may be greater than zero and the variations of D must bd 
taken into account. For atomic latticed materials at low temperatures 


do Sar \ 2/8 e712 HT? 
paar - (=) he ; oo (36) 
whilst at high temperatures 
QGae Te in 
me = 4 Ink T3 orp PP US RM: STI. ° Peis tecerese (37) 
In ionic latticed impurity metals with a high concentration of free electrons, at 
low temperatures 7 * >2 
do = (3n,\? (el, RH\2 1 
—_ - (%) (sr) a: tes (38) 
whilst at high temperatures n* < —2 
do. 225 7elfH 1 
aa — 256 Ink ip Oe | 8  Beyeliate: wile (39) 


If the material has but comparatively few free electrons, and is still an impurity 
metal, then at low temperatures 


—$ —«= —~—————_ @ 5 
o 9m 
whilst at high temperatures all ionic materials will obey equation (39). No 


experimental results have been obtained to compare with the predictions made 
above. 


(v) The Nernst, Ettinghausen and Righi-Leduc Coefficients 


No measurements of these coefficients have been made upon materials other 
than metals. ‘The assumptions made in the theory appear to be as far-reaching 
as those made in the theory of the thermo-electric power, where experiment does 
not agree very well with theory. It would therefore appear that experimental 
investigation may lead to a further understanding of the deficiences of the theory, 
and thus give rise to a more adequate theory. In this section we shall discuss the 
manner in which the present theory predicts that the effects should vary with 
temperature. 

In semiconductors with an ionic lattice, it has already been pointed out that 
both the Ettinghausen and Nernst coefficients should be negligible unless the 
specimens are very degenerate. In semiconductors with an atomic lattice, 

eel! 1 3a 
Ca mk TUR tte (41) and P= nk? (42) 


so long as 7* is negative, which is the usual case. 

In an impurity metal with an atomic lattice, at low temperatures with »*>2, 
the Nernst coefficient Q is independent of temperature, and inversely proportional 
to the concentration of impurity electrons, whilst the Ettinghausen coefficient is 


given by 
8r\28  m 
d 5a (=) in, Tipe | Se eae (43) 
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_ At high temperatures the variations are those predicted by equations (41) and (42). 
| Ifthe impurity metal has an ionic lattice with a large number of free electrons, then 
+ at low temperatures the Nernst coefficient is again independent of temperature, and 
as the temperature rises and »* becomes less than —2 


o= ~15klan' 
32(2m*RT) 2" 


_At low temperatures, the Ettinghausen coefficient P is given by equation (43) but 

with a negative sign. “As the temperature rises P approaches a constant negative 

value. If, however, there are but few impurity electrons in the impurity metal, 
then at low temperatures 

Rls ho[kT 

O=- ane 

3(2m*)12 ; 

and as the temperature rises the value of Q will gradually approach to zero. The — 

Ettinghausen effect is constant at low temperatures when 7*>2, but reduces to 

zero when the temperature has risen such that »*< —2. 
The Righi-Leduc coefficient S in atomic semiconductors is 


— Tate: I 


S = T6(am*RyP TH owiestes (40) 
whilst in ionic latticed semiconductors 
[ecg Gan 


iv? (2m*)1? 
_ In an atomic latticed impurity metal at high temperatures equation (46) would 
hold, but at low temperatures, with 7*>2, 


1/3 
s=—38(=7) = eens} 
In an ionic latticed material with a high concentration of free electrons, at low 
temperatures 7* >2 e oa sel es Salil Mice a 
8a 2m*kT 
whilst at high temperatures 7*< —2 and 
ee ees AN Te (50) 
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_ Ifthe material has a low concentration of free electrons, then at high temperatures 


n*<-—2 and equation (47) will hold, whilst at low temperatures S will have 
one-third the value given by equation (47). 


(vi) Conclusions 


In this section we have attempted to compare the theoretically predicted varia- 
tions with temperature of the phenomena with those found experimentally. It 
has only been possible to do this by deriving approximate formulae for the electrical 
properties which only apply over specific temperature ranges. ‘The accurate 
formulae, given in Appendix I, must be directly compared with experimental 
results, as Putley (1949) did with Si and Ge, before it can be decided that the 
formulae do not explain the observed variations in any particular case. 
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Expressions for the Electrical Phenomena when I= lyn’ 

In monatomic lattices with low concentration of free electrons equation (1) is 
applicable, and J,=1,/T and r=0. 

In ‘classical’ semiconductors with ionic lattices (‘classical’ semiconductors 
being taken as those in which 7* <0), J) =1,(RT)!2e""*? and r=}. 

In ionic latticed materials at very high temperatures, and in impurity metals 
with a high concentration of free electrons at all temperatures, J) =/, andr =1. 

In the following F,,, is written to indicate F,,,(7*). 


_ 167m*eRTI, qin 8 (2st) eee 
er Teeel © a ae @=-5| G5 ia -9*], 
hi aly 3 (27+) FipF ep do feeesks 2[ (27 +3) Fore |? 
ne 2 (r+1)? WE cas o 2m*kT (r+1)F, . 
of Rl, (7+2)(27 +3) Fai Porte — (27+ $)(7 + 1) F Persie 
Ca (+ TPF? ; 


2. aol Fy jg Ek +1)(2r+ 3) FP, Fersiie— (7 +2)(2r+ teat) 
2n-k (@F+1)F, (+2) Ea — + 3)e4 1) FP saF, 
le 
S= — Om* RT * 
2(r + 1)(r+2)(27 + 3) FF 41 Per sip — (7+ 1)(2r + 3) FP For ssie 
—(7+2)°2r+3)F 41° Fora |. 
CFD FAG 42) Fra 04048) Fe Fag} ] 

These relations are all deduced upon the assumption of relatively small magnetic 
fields. 

In the case of ionic latticed materials at low temperatures and low concentrations 
of free electrons equation (4) may be applied. ‘The electrical properties may be 
obtained from the above equations by setting J,=/,(RT)!2 e*/*? and making 
the following substitutions 


(7+ 1)F,=#( Piet 3 Fi) (27 +4) Fo, =3 Fie + $F 1 —iel op 


(7+2)F,=F(Pop + 3 Fo) (27-3) Persie = el ais + PF yp +75F ip 
(7+3)F, =F 504+ $F ie) (27+ 8) For taie = $F s/o + PP a0 + Fi 


AY PYPSEEN Dale Xora 


The various values taken by the functions 4 — G when n* is large and negative 
are shown in the Table. 


1 A B GE D E F G 
ly en* 2—7* 37/8 7/4 ql? /4 7/16 + iqir 
Lon?) 
a pf 3771/2en* 5 __»* 1-00 1-00 0 0 +1-00 
) : re 
AA ; 457 2257 [Sage Sar a's 
I 2en a Sy sae Hae iene gee hh Pe 
07 128 256 300 128 64 


All the formulae when /=/, agree with those predicted by the classical theory of 
semiconductors. 
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ABSTRACT. 'The nuclear magnetic moment is a two-valued function of the spin, and 
Schmidt’s argument is used in reverse to induce the orbital and spin g-factors. g,1s shown 
to be unity in odd-proton and zero in odd-neutron nuclei, as theory has assumed, but 
gg varies from nucleus to nucleus with a mean value about half that found in the free proton 
or neutron. The magnetic spin ratio G, defined as g,/g(p) in odd-proton and as g,/g(n) in 
odd-neutron nuclei, measures the approach of an empirical moment to a Schmidt line; 
it is shown to obey much the same rules in odd-proton and odd-neutron nuclei. 


§1. ORBITAL AND SPIN g-FACTORS 
A PLOT of the magnetic moment against the spin in either odd-proton or 


odd-neutron nuclei shows that all moments lie near to one of two lines, 

which are roughly parallel (Latham 1947). In odd-proton nuclei the 
magnetic moment increases with the spin (Figure 1), while in odd-neutron nuclei 
it is nearly constant (Figure 2) (nuclei lighter than A =20 are excluded). The 
two lines are explained if both spin and magnetic moment are due to a single 
odd nucleon, which can only exist in a 7-state with 1 and s parallel, or in an «-state 
with 1 and s anti-parallel. In developing this thesis Schmidt (1937) assumed 
that the g-factors of orbit and spin are those found in the free proton and neutron, 
but his lines are too far apart, and with many accurate magnetic moments 
available, the argument is better reversed. 
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When two vectors 1 and s combine to a resultant I, the magnetic splitting 
factor g; is obtained from g, and g, by the equation 


_ (1+1)4+04+1)—s(s+1) (+1) +s(s+1)-—d+1) (1) 
sr 21(I+ 1) é 21-1) Bal eae 
When | and s are parallel, write 7=(J—4) and the magnetic moment p(7), so that 

Quan) == (21 Wee ear” oe. cee eee (2) 
while when ! and s are anti-parallel, write /=(7+4) and the magnetic moment 
(a), and obtain: Ea (G\( Len) iene TeSys (3) 


Each equation is linear in the unknowns g,and g,, so that the latter can be obtained 
from each separately by inserting the mean empirical magnetic moments for a 
series of spin values. In the form here given the lines have the same slope g, and 
the same intercept g,, provided that 2,.(7) is plotted against (27—1) whenl ands 
are parallel, and —2(«) (I+1)/J against —(27+3), when they are anti-parallel 
(Figure 3). 
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Figure 1. "The magneticmoments of odd-proton Figure 2. The magnetic moments 
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The broken curves show the values pre- ured in nuclear magnetons. 
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Figure 3. The magnetic moment as a function of the spin in odd-proton and odd-neutron nuclei. 


The ordinate is 24(7) when 1 and s are parallel, and —2 p(«) awe when they are 
anti-parallel. The corresponding abscissa are 27—1) and —(2I+3) 
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The line obtained for odd-proton nuclei shows conclusively that the slope g; 
is unity, as theory (Bohr 1951) requires, while the intercept g, is not 5:58, but 
about half this value. In odd-neutron nuclei the points are less consistent, but 


| _ theory argues strongly that g, shall be zero and this gives as good a fit as any other 


value. The mean value of g, is about — 1-95 or half that found in the free neutron. 

This conclusion might have been reached by a simpler if less conclusive 
_argument, for the direction of the Schmidt lines is determined by g, and their 
distance apart by g,. As the empirical lines are parallel to those predicted, but 
nearer together, the value of g, postulated by Schmidt is correct, but his value 
of g, is too large. 


$2) THE MAGNETIC SPIN RATIO 


When the magnetic orbital factor is taken as unity in odd-proton and zero 
in odd-neutron nuclei, an empirical magnetic moment determines the magnetic 
spin factor g,, provided one knows whether I and s are parallel or anti-parallel. 
In odd-neutron nuclei the distinction is between negative in 7-nuclei and positive 
in a-nuclei. In odd-proton nuclei magnetic moments between (J—4) and 
I(I+14)/(1+1) are ambiguous, but only three stable nuclei have been found in 
this range, and of these !’I and 7°°Bi are so near one limit that the correct 
assignment is hardly open to doubt. The third *As is not well established. 
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Figure 4. The magnetic spin ratio of 7-nuclei. The circles and triangles show the mean values of 
the ratio for all known odd-proton and odd-neutron nuclei respectively. ‘The curves are 
intended only to show the general trend of the values. 


Figure 5. The magnetic spin ratio of «-nuclei. 


The more nearly the empirical g, approaches the values g(p)=5:58 and 
g(n)= —3-83 found in the free proton and neutron, the nearer the magnetic 
moment lies to one of the Schmidt lines, so introduce the magnetic spin ratio G, 
and define it as g,/g(p) in odd-proton and g,/g(n) in odd-neutron nuclei. G may 
be derived from the empirical moment using equations (2) and (3), or explicitly 


Zodd. z-state Gy(p)=p(P, 7)—(7—-4) 
a-state Gy(p)=J+14—(1+1/D)u(P, «) 4 
WW. odd. a-state, aGu(n)=(N,a) ©) 
a-state Gy(n)=—(1+1/Z)u(N, «). 

The magnetic spin ratio G is convenient, because it often obeys the same 
rules in odd-proton and odd-neutron nuclei and lends itself to a concise statement 
of certain empirical rules. 

(i) Among stable nuclei with A greater than 20 the range of G is 0-29 in 7] 


to 0:94 ir, “7Au in edd-proton nuclei, and the range is almost the same in 
odd-neutron nuclei, 0-31 in 43Cd to 0-95 in Pt, 
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(ii) The mean values of G for all nuclei of the same spin and type are plotted 
in Figures 4 and 5. In z-nuclei G is nearly independent of the spin, ranging 
round values of 0:52 in odd-proton and 0-42 in odd-neutron nuclei. In «-nuclei 
on the other hand G falls continuously from 0-9 when the spin is 4 to less than 0-5 
when the spin is 44. 

(iii) When an odd-proton nucleus with 7 protons is compared with an 
odd-neutron nucleus with 7 neutrons, the magnetic spin ratios are closely the 
same (Table). This rule is valid only when 7 <49 and when both nuclei have 


the same spin, but six pairs conform and no exception is known (Schawlow and 
Townes 1951). 


The Magnetic Spin Ratio of Nuclei having the same spin and either 
ry protons or 7 neutrons 


Yr Nuclei I G 
5 °Be UB 14 0-61, 0-61 
CIN 4 1-10, 1:02 
13 25g 27Al 24 0:45, 0:58 
17 “BISA 13 0:55, 0:58 
37 SUG AIR DS 24 0:66, 0:40 
49 SSS EO bra 4} 0:57, 0:53 


(iv) The addition of two neutrons without change of spin has been 
described as pushing the magnetic moment towards the Schmidt lines 
(Wangsness 1950, de Schalit 1950). In terms of the magnetic spin ratio this 
reads G(Z,N+2)>G(Z,N). All known stable isotopes obey this rule, the 
fifteen pairs established being H, Cl, K, Cu, Ga, Br, Mo, Ag, In, Sn(2), Te, Ba, 
Re, Tl. But in I and Cs the addition products are unstable, and G decreases. 

Combined with Bohr’s (1951) rule for quadrupole moments, this means 
that when the addition of two neutrons produces a stable nucleus of the same 
spin, the numerical value of the quadrupole moment is reduced, so that the 
nucleus is more nearly spherical. Five isotopic pairs (Cl, Cu, Ga, Br, Re) 
conform; Eu seems to be an exception, but is easily distinguished because the 
addition of the two neutrons drives the odd proton from a z- to an «-state, while 
in the five conforming pairs the state of the odd proton is unaltered. 
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X-Ray Diffraction by a Crystal Possessing Periodicities 
within the Unit Cell 


By J. W. JEFFERY 
Birkbeck College Research Laboratory, London University 


MS. received 16th January 1951, and in amended form 5th June 1951 


- ABSTRACT. The Fourier transform of a crystal is used to explain the appearance of the 
X-ray reflections from the true crystal lattice of tricalcium silicate as satellites of the strong 


reflections corresponding to a pseudo-lattice. Similar but more general diffraction effects 
| are also discussed. 


§1. INTRODUCTION 


N the course of an x-ray crystallographic investigation of tricalcium silicate 
(Jeffery 1950, 1951) a rather peculiar diffraction effect was encountered. The 
crystal symmetry is strongly pseudo-rhombohedral, with a marked pseudo- 

halving of the hexagonal a-axis, but the real symmetry is monoclinic, one of the 

| pseudo-hexagonal a-axes being the monoclinic b-axis. The monoclinic face-centred 
cell (a=33-1a, b=7'1a., c=18-64., 8B =94°) is twelve times the volume of the 
primitive rhombohedral cell (Hex. a=7'04.,c=25°0a.). The ‘ extra’ reflections 
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Figure 1. Equivalent portions of two zero layer Weissenberg photographs. 
(a) About that hexagonal a-axis which is equivalent to the monoclinic b-axis ; (b) about one of the 
other a-axes. 


corresponding to the monoclinic lattice are all fairly weak and occur as satellites 
of the strong rhombohedral reflections. In between these rhombohedral 
reflections and their satellites are large areas in which the monoclinic reflections 
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are undetectable, so that even the establishment of the dimensions of the monoclinic 
cell was not an easy task. This is shown in Figure 1, which consists of the corres- 
ponding parts of two zero-layer Weissenberg photographs. One, about the 
monoclinic b-axis, shows satellites grouped round each strong rhombohedral 
reflection, the second, about one of the other pseudo-hexagonal a-axes, shows 
only the rhombohedral reflections; both photographs are heavily exposed and the 
satellite spots show up much more strongly than would be the case if the photo- 
graphic densities were proportional to the x-ray intensities. 

A similar effect has been described by Vand (1951) in the more general case of 
long chain compounds. Here, within a large unit cell, regions of order occur 
corresponding to a smaller unit cell, but there is no simple relation between the 
two lattices. In this case strong reflections corresponding to the effectively 
infinite lattice with the large unit cell only appear around the points where 
reflections would be produced by an infinite lattice based on the smaller cell. 

The purpose of this note is to show that these effects can be simply explained on 
the basis of the Fourier transform conception. The results derived by Ewald 
(1940) will be used, and are summarized below. 


§2. PROPERTIES OF THE FOURIER TRANSFORM 

(i) The Fourier transform of a given electron density distribution describes the 
scattering of x-rays by that distribution. 

(ii) In the case of the contents of the unit cell of a crystal the Fourier transform 
is, in general, a continuous function. 

(iii) The Fourier transform of an infinite crystal lattice with points of equal 
weight is the reciprocal lattice, which also has points of equal weight. 

(iv) The Fourier transform of a real crystal, assumed infinite, is obtained by 
multiplying the Fourier transform of the crystal lattice by that of the unit cell 
contents, at each point in Fourier space. It is therefore a reciprocal lattice 
consisting of weighted points. 

(v) The Fourier transform of a small crystal is obtained from that of the 
infinite crystal by associating a function, called the shape transform, with each 
point of the weighted reciprocal lattice. 

The shape transform is the same round each reciprocal lattice point, but it is 
multiplied by the weight of the point in each case. ‘The Fourier transform of a 
small crystal therefore has a volume distribution of scattering power round each of 
the reciprocal lattice points corresponding to an infinite crystal. 

The shape transform can be calculated in simple cases, and in particular Laue 
(1936) has shown how to calculate it for any shape bounded by plane faces. Unless 
the crystal is very small the transform will not have an appreciable value except in 
the immediate neighbourhood of lattice points. 


§35 DHE SE XPLANAT LONGO E ft Heep TE ERAGON Ge Cans 
(i) Tricalctum Silicate 


We now apply the above considerations to the case of tricalcium silicate. Asa 
first approximation consider the monoclinic cell as a number of rhombohedral | 
cells, i.e. consider it as a very small crystal with a true rhombohedral lattice. The | 
Fourier transform of this will be a weighted rhombohedral reciprocal lattice with | 
the shape transform of the monoclinic cell associated with each lattice point. This 
will be an approximation to the Fourier transform of the monoclinic cell contents. 
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| To obtain the transform for the whole crystal, considered as infinite, we multiply 
| the transform of the cell contents by the monoclinic reciprocal lattice, with all 
i) points having equal weight. The transform of the cell contents has appreciable 
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©) Figure 2. Schematic diagram showing the production of satellite reciprocal lattice points in the 

case of tricalcium silicate. 

(a) Regions in which the Fourier transform of the monoclinic cell contents has an appreciable 
value. The shapes of the regions are the same but their relative extents are determined 
by the weight of the rhombohedral reciprocal lattice points. 

_ (6) Monoclinic reciprocal lattice points. 


_ values only near the rhombohedral lattice points. Therefore after multiplication 
only those monoclinic lattice points coincident with or near to heavily weighted 
_ rhombohedral points will have weights appreciably different from zero, and the 
? monoclinic lattice points will appear as satellites of the rhombohedral points 
_ (Figure 2). 

" (ii) The General Case 

In the more general case the contents of the large unit cell are made up of two 
} parts. ‘The first part is the region of order which can be considered as a very small 
/ crystal with a small unit cell. ‘The second part is the rest of the contents of the 
+ large cell. If the former predominates, its Fourier transform is a good approxi- 
1 mation to the transform of the large cell contents. ‘The argument then proceeds as 
) before, except that none of the reciprocal lattice points corresponding to the large 
cell are necessarily coincident with those derived from the smaller one. 


(111) Conclusions 
: These considerations give a descriptive explanation of the observed diffraction 
| effects which can form a foundation for a more detailed and quantitative analysis. 


§4. FURTHER CONSIDERATIONS 
Vand (1951) has shown that even in the more general case the simple structure 
of the small cell can usually be solved fairly easily and the phases of the reflections 
corresponding to this are then known. From a consideration of the Laue 
function giving the shape transform of the very small crystal contained in the true 
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cell, it is then possible to deduce the phases of the ‘satellite’ reflections and thus 
solve the complete structure. 

In the case of tricalcium silicate it can be shown (Jeffery 1951) that the increase 
in relative intensity of the satellite reflections, from zero at the origin to comparative 
equality for the high orders, corresponds to a periodic distortion of the rhombo- 
hedral lattice rather than a periodic variation of density from cell to cell within an 
undistorted lattice. 
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A Search for Sudden Large Changes of Energy of Beta-Rays 
in a Cloud Chamber 


By W. T. DAVIES anp D. F. SHAW 
The Clarendon Laboratory, Oxford 


MS. received 11th Fune 1951 


ABSTRACT. Measurements are presented of the changes of curvature of f-ray tracks in 
a cloud chamber containing air at approximately N.T.P. situated in a magnetic field of 
400 gauss. These are examined for evidence of processes involving the reversible 
production of electrons of rest mass a small multiple of the normal value. Indications from 
a study of earlier investigations, that such processes might exist, are not substantiated, 
and the results confirm and extend those of Bleuler. It is shown that this type of cloud 
chamber investigation is unsuitable for obtaining useful information on inelastic scattering 
of fast B-particles, and that multiple scattering and observational bias are adequate to explain 
the anomalies previously reported. 


Sa UNALR ODUGIE ON 

LOUD chamber experiments have yielded conflicting results on the 
inelastic scattering of fast f-particles. Several investigators* have 
concluded that scattering with an energy loss of more than about 20% 
occurs with a frequency 10 to 100 times that predicted by the theory of Bethe - 

and Heitler (1934) in which the loss of energy is due to emission of radiation. 

This has led to suggestions that the excessive frequency of sudden energy losses 
is due to the occurrence of an additional process having the appearance of a. 

radiative collision, but involving a third particle which must be non-ionizing. 
Measurements of the energy distribution of the x-radiation produced when . 
electrons of 1 to 3 Mev. are absorbed in different materials have confirmed the : 
accuracy of the Bethe-Heitler theory in this energy region. Calorimetric : 
experiments have shown that, in the absorption of fast electrons, any process } 
involving the emission of non-ionizing radiation, if it exists at all, must have al 
probability small compared with that for bremsstrahlung production. | 
The most recently reported search for inelastic scattering in the gas of a) 
magnetic cloud chamber is that of Bleuler (1942) who found no evidence for an} 
excessive number of cases of sudden energy loss in fluorine. While 19 cases: 


* A bibliography of earlier work is given under general headings at the end of the list of references. | 
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of single scattering were observed exhibiting apparent losses of over 10% of the 
incident kinetic energy, 25 cases were observed in which scattering was accom- 
| panied by an apparent gain of energy of more than 10°4. The expected number 
} of losses of this amount of energy or more, by radiation, was about three for these 
experimental conditions. Bleuler accounted for his observations on the assump- 
tion that the changes of curvature were not in fact due to changes of energy, but 
were simply the result of multiple scattering of the electron in the gas of the 
/ cloud chamber. Although this effect was encountered at an early stage in cloud 
» chamber development (Wilson 1923), it was left to Williams (1939) to present the 
first quantitative treatment. 

If the momentum of the electron is determined from the curvature of the 
track in a magnetic field, multiple scattering will introduce an error in the 
measured value of this radius of curvature (Bethe 1946). This error increases 
with increase of gas pressure and with reduction of the length of the track used 
in the curvature measurement (Bethe 1946, eqns. 12 and 13). Apart from the 
fact that the effect of multiple scattering on curvature is less likely to cancel out 
for a short track than for a long one, the curvature measurements themselves 
will be less precise as the track length is reduced. Thus, if a single event such 
as a nuclear deflection divides a track into two parts for the purpose of analysis, 
it may frequently happen that a change of curvature is recorded. 

Although one may explain the failure to observe cases of curvature decrease 
in other work than that of Bleuler as due to an observational bias which causes 
them to be rejected as of no significance, there is another feature of previous 
work which has not been satisfactorily explained, viz. the occurrence of cases 
of almost complete energy loss (see for instance Williams and Terroux 1930, 
Barber and Champion 1938, Davies 1939). In such cases, when the track 
observed after the process has been very short and heavily ionizing, the question 
of a multiple scattering error in the curvature hardly arises. 


§2. THE PRESENT INVESTIGATION 

It is conceivable that multiple scattering may not be the only process which 
could give rise to a decrease or an increase of curvature following a nuclear 
) collision. A process which might be imagined to take place is one in which 
| the rest mass of the electron is increased at the expense of its kinetic energy 
| relative to the nucleus. Such a process has in fact been previously suggested 
__ by Zwicky (1938)* and it is of interest to consider the implications in more detail. 
| Let E and p represent the total energy and momentum respectively of the 
i primary electron of rest mass my (measured relative to a colliding nucleus). Then 


; Doel iC ineCee 5 ow OO Meee (1) 
If the electron undergoes an elastic collision with the nucleus in which its 


' rest mass changes to nm, (n>1) and its momentum decreases by an amount Ap 
with the recoiling nucleus contributing to the momentum conservation: 


(B—Ap)*= BAe? = (niny)Pcte (2) 
Solution of equations (1) and (2) gives: 
Oe = alos} ia we nee (3) 


taking the negative square root since p>Ap>0 for n>1. 


* The possible existence of a spectrum of rest masses for the fundamental particles has also 
been suggested by Moller (1947) and Bhabha (1947). Also Janossy and McCusker (1949) have 
discussed the existence of cosmic-ray particles of rest mass 3m, and 10mp. 
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To comply with the results of the calorimetric experiments it would be 
necessary to assume that the inverse process can occur, i.e. the increase of 
translational kinetic energy with decrease in rest mass. For this case, taking as. 
the initial rest mass nm, and the final value m, the corresponding increase in 


momentum Ap’ is S Ap’ a {p2-+ (n? Ee 1)m,2c?}? —p “a. See (4) 


with the requirement Ap’>0 for n>1. 

In Figure 1 are plotted curves showing the relationship between Hp and 
A(Hp) on the assumption that the rest mass can increase from my to 2, 3 or 4m, 
or can change in the reverse sense. Included in this figure are the published 
results of Barber and Champion (1938) and Leprince-Ringuet (1935), chosen 
because they were in a form suitable for this analysis. It will be seen that, while 
this evidence cannot be said to corroborate the hypothesis, it certainly would 
not exclude the occurrence of such a process. 
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Figure 1. 


Besides the obvious discrete values of possible changes in momentum, the 
distribution shown in Figure 1 differs from that expected to arise from multiple 
scattering in two other respects which should be capable of experimental test: 
(a) for initial Hp values below 2,800 gauss cm. one would expect to observe only 
increases in momentum, whereas multiple scattering errors could contribute 
apparent increases and decreases; (b) for Hp values above 2,800 gauss cm. both 
increases and decreases could occur with the magnitude of the decreases 
significantly greater than that of the increases, whereas multiple scattering 
would give approximately equal positive and negative changes. 


§3. EXPERIMENTAL WORK 

In an effort to obtain further experimental evidence, particularly since 
Bleuler’s data were not given in sufficient detail for us to analyse in the manner 
indicated, we have made new measurements on the curvature of electron tracks 
in a cloud chamber situated in a uniform magnetic field of about 400 gauss. Air 
at approximately atmospheric pressure was used as the cloud chamber filling. 
The measured tracks of electrons from a **P source covered the range of Hp 
from 2,300 to 6,600 gauss cm. This is a region in which the effect of an integral 
change in rest mass as postulated would be most pronounced. 

The photographs were projected to original size on a screen on which were 
drawn reference circles of radii from 6cm. to 15cm. at lcm. intervals. The 
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| radius of curvature of each track was determined by fitting to one of these reference 
circles. Our criterion for circle fitting was that there should be no radial 
} deviations exceeding the width of the image of the track (4 to 1mm.) from a 
_ circular arc at least 5cm. long. This length was adopted, after trials by different 
observers on a selection of tracks, and by the same observer at different times, 
as the minimum which allowed a decision to be made between adjacent circles. * 
) The length and radius of curvature of each track, and single deflections greater 
/ than 3° were recorded. Examples of cases when a track did not fit the selected 
' circle along its whole length are illustrated in Figures 2, 3 and 4 (Plate+) and the 
, method of treatment detailed below. 
: (a) In Figure 2, track abc, a nuclear deflection of 4° occurs at b. It is 
necessary, therefore, to measure the radii of curvature of portions ab and bc 
separately. 

(5) In Figure 3, track defg, the portion ef is almost straight due to plural 
scattering. This results in a division of the track into two sections and again 
circles have to be fitted separately to the portions de and fg. Since the portion 
ef is less than 5cm. long, no attempt would be made to measure the radius of 
curvature of this section. 

(c) In Figure 4, track Ajk is typical of several found in which the curvature 
appeared to change continuously along its length. In such cases it was usually 
possible to represent the track very closely by two arcs of circles of different 
radii intersecting on the track. It often happened that the two selected circles 
best fitted the track when intersecting at a small angle (2° or 3°) so that in 
representing the track in this way a small deflection was introduced where none 
was at first apparent. It is not considered that this affects the significance of the 
results. 

The measurements covered 350 metres of track in the selected range of 
momenta; 252 cases were found in which a change of 1 cm. or more in the radius 


Table 1. Distribution of Track Length and Changes of Curvature 


! Total 
+ p (cm.) 6 7 8 9 10 11 12 13 14 15 Gals 
Hp range (gauss cm.) 2300-— 2700- 3100- 3600- 4000- 4400- 4900— 5300- 5700- 6200- 2300- 


2700 3100 3600 4000 4400 4900 5300 5700 6200 6600 6600 
|Actual track length (m.) Bail Sei GG Rey Gt) a) Se i eS Ses 
| Effective track length (m.) 52a Oe Sie 15 jeaen Lec Siee 14 Oe 112-6 ss 3°8 tO} 97025 90 


us 6 Tota 27) W300" 26 * Pt9" . 416 6 B 4 0 0 131 
— Dh oauma:?) 
Sy 495 Ae>4tem. 1 0 4 0 1 0 0 0 0 0 6 
+ aes Ap=4 cm. 1 D 0 0 0 0 0 0 0 p é 
tae 0 Ap=3 cm. 4 5 5 2 0) il 0 0) 0) 
Me st Ape tom) 11 as) a2 me) Re 000 56 
5“ TS Ap=icm. 10 8 5 feet 7 3 1 4 0 0 49 
¢ 
S- ie Ap=t'em, ©)'7 9 8 10 6 4 0 0 0 51 
fee o£ Ap=2cem.. 3 6 Sepliet Ss 5 Fe aytO 2 0 42 
Pee O53 Ap=3 cm. 0 1 2 4 5 2 2 1 0 1 18 
Efo¢gs 
BS EEE Ap=4em. 0 O26 1 0 0 1 1 9 1 4 
cya) 333 Ap>4cm. 0 1 0 0 0 2 1 1 0 1 6 
| a iotalee OMe i7ee e413 223° OF" 15: 12 3 2 3 121 
Number of tracks showing 
| no change (Agmen OSE 306) 9325-235. 175% 106 53. 33 She logs 


* Bothe (1948) has discussed various other criteria for circle fitting. We do not consider the 
difference between them to be significant when the radial interval between the circles to be fitted 
is as large as in our case. 

+ For Plates see end of issue. 
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of curvature of a track was observed. The distribution of track length and 
changes in curvature 1s shown in Table 1. Any sudden change in curvature 
occurring within 5cm. of either end of the visible portion of a track could not 
meet our requirements for acceptance. ‘Therefore, in computing the effective 
track length (row 4) the first and last 5cm. of each track have been excluded. 
In addition, the track length has been reduced from the experimental conditions 
(65cm. Hg pressure at 16°c.) to the equivalent at N.T.p. The effective track 
length for curvature changes was 90 metres. 


§4. DISCUSSION AND CONCLUSIONS 

It is evident from the results that cases of increase and decrease in curvature 
occur with equal frequency. ‘There was no apparent periodicity in the curvature 
changes such as would have accompanied changes in rest mass. In the lower 
momentum intervals the number of cases of curvature decrease significantly 
exceeds the number in which it increases. This may be attributed to the fact 
that the electrons producing these tracks entered the cloud chamber at an angle 
such that their path was convex to the chamber radius. If then, a change 
occurred to increase the curvature their path before striking the chamber wall 
would be less than if the curvature decreased. ‘The requisite minimum track 
length required for circle fitting would therefore be more readily achieved in the 
latter case. ‘This effect was less pronounced at the higher momenta where the 
electrons entered the chamber in a different direction. 

We verified that the observed changes were not produced by lack of 
homogeneity in the magnetic field or by distortions in the chamber, the cameras 
or the reprojecting system. ‘The spread of curvature changes is much greater 
than would be expected from our errors in circle fitting. We conclude therefore 
that most of the changes result from the effect of multiple scattering since this 
is the only other factor which might produce an almost symmetrical distribution 
in curvature about the mean value. The results of this experiment do not exclude 
the occurrence of a process, such as an increase in rest mass, with a frequency 
small compared with the random changes of curvature.* A more direct way of 
investigating this effect would be to observe the collisions between primary 
B-particles and atomic electrons. Ruark and Jones (1938) analysed the results 
of Champion (1932) and found no evidence for a change in rest mass, but the 
number of collisions observed was comparatively small. 

We analysed a selection of our tracks using the method of Groetzinger, 
Berger and Ribe (1950), by measuring the angles between successive chords 
1:5cm. long on 34 sections of track varying in length from 6cm. to 15cm. The 
root mean square deviation from the mean angle varied between 0-5° and 2:5°; 
the corresponding value obtained for measurements on true circles was rather 
less than 1°. In the case of an electron of initial Hp =5,000 which experiences 
a change of rest mass, the predicted values for.the root mean square angle of 
multiple scattering « in air at 65cm. Hg pressure are given in Table 2. These 
results show that this method cannot be used to detect a change of rest mass 
under our conditions. 

In the track length we examined, the Bethe—Heitler theory predicts about 
one case of radiation loss greater than 20°, whereas over 70 changes of curvature 


* The number of cases found of a change in radius of curvature corresponds roughly to the 
expected number of single scattering events in which the projected deflection is greater than 5°. 
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both positive and negative were observed corresponding to an energy change 
greater than thisamount. In order that a difference of 2°, between the number of 
apparent gains and losses may be statistically significant one would need sufficient 
track length to yield about 5,000 such changes. Under our conditions this 
would mean measuring nearly 15 kilometres of track, a tedious investigation 
unless assisted by such improvements in cloud chamber technique as those 
introduced by Champion and Roy (1948) and Hodson, Loria and Ryder (1950). 

The last point for discussion is the occurrence of cases of almost complete 
energy loss of which we observed three. If plotted on Figure 1, all three would 
lie close to the curve for n=2. We prefer, however, to regard them as due to the 
chance production of a 6-ray as the f-particle passes out of the illuminated region 
of the cloud chamber. A vertical section of the light beam gives a region, fairly 
accurately 1 cm. deep, of intense illumination bounded top and bottom by strips 
approximately 2mm. deep in which the illumination falls to a value insufficient 


Table 2. Predicted Root Mean Square Angles of Multiple Scattering for 
electrons of different rest masses in air at 65 cm. Hg pressure and 16° c. 


Initial conditions Final conditions 

Hp Rest mass a Hp Rest mass oa 
5,000 My de 4,000 2m ay" 
5,000 2m 25 6,000 Mo 1 Lots 
5,000 3m EY 6,500 Mo ae 


to photograph 6-rays. It is reasonable to assume that a 6-ray produced within 
the last centimetre of a track which is passing out of the light beam because of 
obliquity may be mistaken for a case of very large energy loss. We have counted 
in 27 metres of track 18 6-rays of projected length greater than 5mm.* The 
total number of tracks which passed out of the light in crossing the chamber 
was about 1,300. Hence one would expect roughly nine cases of this kind. 
Clearly it seems unjustified to consider any other process. 

We conclude that there is no reliable cloud chamber evidence that the 
frequency of occurrence of significant energy losses due to radiation is greater 
than predicted by theory. The results of most of the earlier cloud chamber 
experiments are invalidated by the limitations on the method imposed by 
multiple scattering. 
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ABSTRACT. ‘The theory of the first order Raman effect in crystals is presented in a form 
which takes into account the finiteness of the wavelength of the incident radiation. The 
theory is applied to calculate the relative intensities of the first order Raman line and of the 
Brillouin components of diamond, which are interpreted as the Stokes and anti-Stokes 
Raman lines produced by those elastic vibrations, whose wave vector Q fulfils the condition 
Q=Q’—-. (Q and Q’ are the wave vectors of the incident and scattered radiation.) 
The calculated intensities are in good agreement with experiments of Krishnan. <A formula 
derived by Leontowitsch and Mandelstamm is discussed and is shown to be inadequate to 
calculate the intensities of the Brillouin components. 
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§1. INTRODUCTION 
CATTERING phenomena in crystals can be divided into two groups: 
scattering with frequency changes and scattering without frequency 
shifts. A process of the latter kind occurs only with radiation whose wave- 
length is of the same order of magnitude as the lattice constants*. For longer 
waves, the scattered radiation which is of the same frequency cancels out by 
interference, if the crystal is perfect. In this case only those scattering processes 
_ yield a finite intensity, which are connected with vibrational transitions producing 
| frequency shifts between the scattered and the incident light. These latter 
processes are all of the nature of the Raman effect. They are called firsc order 
effects if only transitions v->+v +1 are involved. (Here v denotes the vibrational 
states of the crystal.) 

Max Born (Born and Bradburn 1947) has developed a general theory for the 
Raman effect in crystals, which yields explicit formulae for the intensity of the 
Raman spectrum of the first and second order. This theory has been found to be 
in good agreement with experiments of different kinds. The theory implies a 
simplifying assumption which can be completely justified for the Raman effect 
with optical waves; it is however not admissible when dealing with some of the 
phenomena which have been investigated experimentally during the last few 
years. Born’s approximation consists in assuming that the wavelength of the 
incident light is long compared with the dimensions of the scattering crystal, 
i.e. practically infinite. 

In this paper the theory of the Raman effect 1s presented in a form that takes 
into account the finite wavelength of the radiation. This allows us to investigate 
theoretically the Brillouin scattering of diamond, which appears to be the most 
important of the problems for which the finite wavelength has to be taken into 
account and which justify the increased complexity of the formulae. ‘The 
Brillouin components have been discussed previously by other authors; their 
treatment is however not quite satisfactory. 

The Brillouin components of diamond have been investigated by Krishnan 
(1947). Using the elastic constants of diamond, he calculated the frequency 
shifts and found good agreement with experiment. In order to derive the 
intensity values from the theory, however, he applied a simple formula of 
Leontowitsch and Mandelstamm, which is not adequate to give all the required 
details. Further, this formula does not allow calculation of the intensity ratio 
between the Brillouin components and the first order Raman line with frequency 
1332cm™!. These deficiences are removed in the present paper. 


§2. LATTICE DYNAMICS 
In this section we summarize some of the results obtained previously by Born 


(1942) and Born and Begbie (1947). 


The position vector of a lattice particle may be written 


r(,)=FO +e) p(t) 


where £0) ait lj Asclsaces) | PLN ale otos (2.2) 


* This is exactly true if the internal forces and the polarizability changes are linear functions of the 
displacements of the single atoms. The unavoidable non-linear effects are, however, negligibly 
small in most cases of light scattering in crystals. 
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here /,, /;, 1, are integers, a,, aj, a, are the three elementary vectors of the lattice 
cellandr(k)(k=1,2, ...p)are the position vectors of the p particles in the cell. 
l l 
(.) are x, (1 (Ca 
The second derivatives of the potential ® of the disturbed lattice will be denoted, 


J—I' PD 
P.p( hk! ) = aun ray ll Pp UE PS Rowse ve (2:3 
*(,)%(0) 
‘The equations of motion for small vibrations, of amplitude 
iB l l l 
me = a0 
e(,) 8 (,)/vme=X @) Xx, (.) eer (2.4) 
ey) J—I’ iM 
ae 3. () =-= Das( ae ) o, (.) har Mii tua oy se (2.5) 
l l 
where Dip (ay =(mmy)"? D4 Ge Pee a oe ea (2.6) 


The equations of motion (2.5) are satisfied by plane waves which may be 
written in the form 


The rectangular components of r 


by 


Ux (.) =) .(k) exp, (tut exp (qi.til) tn ual aL saan (2.7) 
or in the form Obs (.) = Wk) exp (—iwt) exp ii4 D3 (.) } St fueetid: eet (2.8) 
where WARDS ARiexp iid. t:(2) ei sieltialtl ait. aise (2.9) 
The equations of motion can be expressed in the compact form 
OA = Cg) We totes (2.10) and OV = DG) ere Cat) 
where Q =w? and W and V are one-column matrices : 
W={W,(1), W.(1), W3(1), W4(2)...,Wa(p)}. «eee (2212) 


The elements of the matrices C(q) and D(q) are 
l PAH 
(Nee , eee , = . 
Ca(q| RR ) D,(q| RR )exp U iq.r(kk )} *Dap (a exp {4-7 (sy) f 5) 
With!” tke )="(RV tk). 8 ge ee ene (2.13 


It is convenient to introduce normal coordinates, with the help of which the 
expressions for the kinetic and potential energy are transformed into sums of 


TT (0) vo) elon) 0 


The e(k|") are the components of the eigenvectors of the lattice, i.e. e{ 2 1\ig 
the displacement v(k) of the atom : when the lattice performs the normal 


vibration ( :) . They satisfy the orthogonality relations 


s(p)e(p) aan salon 


4 = 3(«B)3(RR’) 
and the equations of motion om 2 ese (2.15) 


Q(j)e(j) =D @ gh a nr cee (2.16) 
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For long elastic waves the amplitudes W(k)//m,, as defined by (2.8), (2.9) 
do not depend to the zero order in q on & and satisfy the ordinary equations for 
| waves in an elastic substance. ‘These can be written (p is the density of the 
os) Pe DiGi ee eee 7) 


where the elements of D'(q) are referred to the elastic constants c,,, by the matrix 
equation 


D'3(9) Ci 66 55 C65 Pes C16 qi 
D'22(9) Cog Con Cas C24 C46 C62 qo” 
D'33 (9) _| 5 4a &s3 Cag C35 C54 q3° 
D's3(q) Cos C24 Cag 3 (Cogt Cas) (Cas +56) 3 (Coa + C25) 29293 
D'31(9) C51 Cag C35 b(Cas tag) (C31 +655) 3( C56 + Cra) 24391 
D'12(9) Cie C2 C54 3 (Coat Cos) (C56 + C14) 3(Cr2 + Coo) 29192 
ee (2.18) 


To the first order in q the displacements W(k) of the different atoms in the 
elementary cell differ a little from each other. ‘This difference must be taken into 
account in calculating the intensity of the Brillouin components. The first order 
terms W” are given by the matrix equation 


he EG) Way ik oa ee aie hace (2.19) 
l 
(ON bp’) = 
where Cya(RR’) = 2D ag Co Fen cen (2.20) 
ety od l 
(1) Ye eel 
Cyn(q|RR’) = IBD og (e,) cS Gs geese (2.21) 


In a finite lattice the wave vector q is restricted by the cyclic lattice condition in 
the following way: 


We define the basic vectors b, of the reciprocal lattice by 


(lifa=p 
b= ae) Bee, Ys ANS) 
aoe | Vif aXB, a 
and write Q=21q,b,+27q.b,.+27gsb3. §= «20. (2.20) 
Then q.r(/) =2i(qyl,+ Yolo + sls) Meletatele (2:21) 
and 
h h h 
qs me a= Feb q3= * (Rate le 2, (t= 1))\o ae (2.22) 
where n? = N is the total number of atoms in a simple lattice. 
Equation (2.22) ensures that for any function of the lattice points 
SFC OHARA RY Ses OS ee (22 rer AB} ee Re. (2.23) 


§3. GENERAL THEORY OF THE RAMAN EFFECT IN CRYSTALS 

According to Placzek (1934) the scattered intensity corresponding to the 
transition between two vibrational states v, v’, is proportional to the matrix 
element [M]?,,. where M =cE is the electric moment induced by the electric field 
E of the incident wave and « is the polarizability tensor of the scattering crystal. 

M can, in first approximation, be considered as the sum of the contributions of 
the single atoms. In each atom the external field E acts with a different phase, 
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depending on its position. Further another relative retardation must be taken 
into account, because of the different distances of the scattering atoms from the 
observer. 


Let O28 Ol S287 [sale l= eee (3.1) 


be the wave vectors of the incident and scattered wave (neglecting the change of 
wavelength 4) and let 


Q=Q’—Q’ =27(s"—s')/A, Q=|Q|=(47/A)sindy, <..... (312) 
where y is the angle of deflection. Then an elliptically polarized harmonic wave 
E=R[Aexp{—i(wt—Q’.r)}] =4[A exp {—i(wt— Q’.r)}+ A*exp {7(at —Q’.r)}] 
Hee (333) 


produces the retarded moment 


M,= ee (|) exp (iQ. E () +u(,) |} Ao-Beede (pro =e 
(3.4) 


The polarizability of the atom _ depends on the displacements u, which is 


indicated by the symbol used. 
If s is a unit vector normal to s” the intensity becomes 


w= (0+ Oy IM Slip = (0+ yp)! BE anual Ae Ao sa eee (3.5) 


where lore im eee = [a oo putas lees eeceee (3.6) 


Since w,,,<o for visible or ultra-violet incident light, we may omit (w+ ,,,,)4, 
which has practically the same value for all Raman lines. The whole scattered 
intensity is now obtained by forming the thermal average over all the initial 
states v: 


[2 0, wr vv’ exp ~~ (e »/R ie) 


a) Pe ae exn( Ses ene €, =the energy of the state v. 


For natural incident light and scattered light observed without an analyser, we 
must average over the polarization of the incident and scattered radiation : 


YAl AL se sf : ; Cae jLaee 
0 i Se (po) oo of So N p,o= 12 Ce (oe et (3.8) 


fo= —B8e'so", (pxo), 5,8, =4(1—s,"2) | 


I, =4 E” is intensity of the incident light. 
§4. EXPANSION OF THE POLARIZABILITY IN TERMS OF 
NORMAL COORDINATES 


We expand the polarizability tensor in powers of the small displacements u (a : 


I l l 
tee(;,|4) =a.) (.) +a) (.) rs Anata 
I ’ = / / 
where a) (,) = 0k), Cone @ = Bes (es )u, ae See (4.2) 
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_ The zero-order terms do not depend on the cellindex/. ‘The first order terms 
| are Ses as produced by the contributions of the displacements u fe of any 
z point 3 to the change of the polarizability of a given atom z) the coefficients in 
the corresponding sum depend of course only on the difference (/—/’). 

The total polarizability of the perturbed crystal is 


=z [ae +002 (3) Jo (0-[eG) +46) 
Eile zzeeeCae ele) Jor el) eG) TI. 


AsQ.u () <1, put exp (iQ. u (.) } =1+7Q.u ) and express u in terms 


of normal coordinates, according to (2.14). Omitting the constant term «‘) and 


the product «{? (.) u (.) one obtains 


bug) Vie 


/ ai esate 
+2 GatD.())antaar fon o-()}e@) 


Introducing r(/')=r —r( _) one can immediately see that «{) is a 


k k 
5-function of (Q+q). With /—/'=/ one gets 


uate ONE ag epee oy 


where Cog QE =N ES vase Ore (i | =) exp {2Q. r(k)} 


go pea| aeo(F) 6 AG ; 1) exp fig. r(/)} 


wk 


l 


sete), (410 cnfa-+(‘)} a 46) 


Equation (4.6) shows that only those lattice vibrations are active in the Raman 
effect whose wave vector q equals — Q given by (3.2). ] 

For visible and ultra-violet incident light the scalar product Q.r (.) lath 
Zis restricted to the small values +1, +2, +3,.... belonging to the atoms which are 
not too distant from the atom k. This restriction is admissible as the «,,,_,, ian 
fall rapidly with 7. Hence one can expand the exponential and neglect terms of 
second or higher pee inQ: 


na) v3 | 0 
os. (ie 


l 


uU——_,-— 

A 
c 
A 
A 

“— 
iS 
NI 

— 


i 
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The physical meaning of (4.7) becomes more evident if the relative difference 


: Ral\, Ree cea je (*l 5) 
oy en(|8) | vm u,(«]2) 
Bes) 


of the displacements of the different atoms in the elementary cell is introduced. 
Equation (4.7) is transformed into 


(2) sata 
Brel) [iv () Js. (2) ]} 9) 


In the second term of this equation all terms not depending on / and k’ can be 
eliminated with the help of the general invariance relation 


i 
E 66 G5, dg! AN soo See (4.10) 


which follows from the invariance of «,, against rigid translations of the whole 
lattice. Hence, (4.9) becomes 

i) 

J 


X60 (e ) =NES {ini0(h)9, +E 664 ( fe) Ee .r(J) +A, (a 
2) | hu. (]2). Saat (4.11) 


+iQ.r(,)d, (x 


This expression behaves very differently for elastic vibrations (index e) and 
optical vibrations (index 0). For elastic vibrations it is convenient to express 


aN (x A as the sum of two terms: 
| =A (ix 


e 


A, (ax 


g ) + AD (i 
j 


e€ 


ae Recs (4.12) 


e 


Here A) ( kk’ | = represents the contribution of the zero order displacement 


Je 
W of equation (2.17) and A‘? (x8 =) the contribution of W of (2.19). Both 


Je 
terms are simple functions of the wave vector Q. 


As mentioned in § 2 after equation (2.16) the zero order amplitudes W(k)/4/m;, 
are independent of k; hence u (: °) =V(k)/\/m,=(W/+/m,) exp{ —7Q. r(k)} 
depend on k only through the exponential factor and one has 


=a) =iQ.1r(kk’)=iEO,r,(RR’). eee (4.13) 


A) (x 


e 


The second term of (4.12) has according to (2.19), (2.20), (2.21) the form 


we (x8 | : ) =t21 Oty (i iB ) ; uy, (: 2) ae ae (4.14) 
pee v7 e e 
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Substitute (4.13) and (4.14) in (4.11) and neglect all terms of second or higher 
order in Q, then 


| 
| 2 =i1N> (0) y i 
& Koc te =iN a Ooo), 1 Be Koo, u ee OF, A 


deals nie (i 
‘ 


i) [tC i) [boo 
=) , then 

| tea( ©) = Bote () (2), ee (4.16) 
ws) ors mo. 32L nl) ri 


+ 2ttee, ze a CQ | } Apr ree (4.17) 


This is the general nats for the polarizability changes connected with 
elastic waves. It consists of two parts: the second term (summation over », /, k’) 
describes the effective changes of the polarizability of the single atoms, the first 
the variations of the phase relations during a vibration. It will be seen in (6.47) 
that both terms have the same order of magnitude and that the phase term is 
indispensable to obtain agreement between theory and experiment. 

In the case of optical vibrations there does not exist a simple relation between 


the A, (ix Q ) is of the 


order 1, i.e. 10% times bigger than for those long elastic vibrations which produce 
the Raman effect with ultra-violet and visible light. One can, therefore, neglect 
in (4.11) all terms containing Q and obtain 

Q ) eet rar (4.18) 


Ore y 5 
bes PI eal ta hk’ A,,{ RR fe 


Summarizing, the stituation is this: With regard to the scattering power the 
essential difference between the optical and the elastic vibrations is represented by 


the different values of A,, (ie a) and A, (a fl ) . For optical vibrations the 


e 


ae Dea, (4.15) 


Arrange with respect to the w,, (i 


5 ) and the wave vector Q. One can only say that A, (i 


local distortion of the lattice, measured by A,,{ RR’ 


is large and the scattering 


power is strong. For elastic vibrations the local distortion is nearly zero and the 
scattering power at 0°K. is extremely weak. ‘This difference, however, is less 
pronounced with increasing temperature, which can be shown as follows: 

The thermal average (7,,.,,,)ay for a single transition vv + 1 of the vibration 


ta is, with the notation of Smith (1948), 


1 for Stokes lines 
One wav =2C (©) {He (*) oe ea 
ba ; : : exp (—hy,/kT) for anti-Stokes lines, 


where c (2) = eas, ts ee Aes pie, hires (4.20) 


IP OmNv ) ass exp] —tw() [ar 
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Hence, for the elastic branches, where hv «kT 


ail =kT [2anve(?). Pa (4.21) 


e 


According to (4.21) the intensity of the Raman effect with elastic waves is 
proportional to T and increases considerably if the temperature is raised. On the 
contrary, for the high-frequency optical waves, which are only weakly thermally 


excited at room temperature, C &) of (4.20) does not depend very much on T and 


there is no strong increase of scattering intensity with temperature. ‘This 
different temperature dependence favours the intensity of the Raman effect with 
elastic waves as compared with the optical Raman effect and compensates the 
fact that the local distortion is much bigger for the optical than for the elastic 
vibrations. 

Further, according to (4.21) the Raman effect for elastic waves should vanish. 
if T=O0°K. But as the condition hy<kT loses its validity even for the low fre- 
quency elastic vibrations if T+0, (4.21) is no more applicable andonly (4.20) gives a 


correct expression for the intensity. Here C( * } 40if T=O0°K. 


e 


§5. THE BRILLOUIN COMPONENTS OF THE SCATTERED LIGHT 


As mentioned in the introduction, there is no light scattering with unchanged 
frequency in perfect crystals if the wavelength of the light is greater than the lattice 
constants. Instead a set of six lines will appear in the scattering spectrum, which 
forms two groups symmetrical with respect to the frequency of the incident 
radiation. ‘The frequency shifts Av are usually extremely small and can be detected 
only by very subtle interferometric methods. Only under particularly favourable 
conditions does Av assume larger values which, however, do not exceed 10 cm~1 
for ultra-violet incident light. Such favourable conditions actually occur if the 
index of refraction and the eigen-frequencies of the scattering crystal are high. 
This happens in the case of diamond. ‘These lines had been predicted by 
Brillouin in 1922, six years before the discovery of the Raman effect, from an 
ingenious consideration, and they are called the Brillouin components of the 
scattered light. 

Brillouin interpreted the effect as reflection of the light waves on the elastic 
waves with Doppler shift. He showed that, in analogy to the Bragg condition, the 
wave vectors Q, Q’ and Q” of the elastic, the incident and the scattered wave must 
fulfil equation (3.2); further that the frequency shifts can be explained as Doppler 
effect, the reflecting waves moving with the velocity of sound. This explanation * 
implies that the thermal excitation of the eigen-vibrations is a necessary condition 
for the Brillouin scattering. 

The explanation suggested here regards the Brillouin components as Raman 
lines produced by the elastic waves. ‘They are the Stokes and anti-Stokes Raman 
lines of those vibrations of the three elastic branches of the crystal, whose wave 
vector Q fulfils the condition (3.2); the small number of Raman lines is due to the 
fact expressed by equation (4.6) that most of the vibrations of the crystal are 
ineffective, due to mutual interference of the scattered light. 

The treatment of the Brillouin components as Raman lines allows a rigorous 
calculation of their intensities. ‘The attempts to obtain the intensities from 
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) Brillouin’s own interpretation have so far not been quite successful. An intensity 
_ formula has been derived by Leontowitsch and Mandelstamm (1931), namely 


are 2 2 2 
T= Fy kTns { Ab tat) +2}, aoe. (5.1) 


¢ Tee Cy tly t+ 2lq, Cay 


_ where d is the wavelength of the incident light, m the index of refraction, T the 
temperature, c,, = elastic constants and p,,, Pockel’s elasto-optical constants. 

| This formula has been developed by an immediate extension of Einstein’s 
_ Scattering formula for liquids to crystals. It is based on classical continuum 
theory and on classical electrodynamics and statistics. The applicability of the 
theory is therefore confined to temperatures for which the classical approximation 
is still valid. According to (5.1) the intensity should vanish at 0°k., whereas the 
_ correct quantum mechanical formula predicts a finite intensity, which is in the 
_ case of diamond about 50 times smaller than it is at room temperature. 

A further simplification of Leontowitsch and Mandelstamm consists in 
_ neglecting the time dependency of the thermal fluctuations in the crystal. Hence 
' there is no frequency shift between the incident and the scattered light and the 
_ frequencies of the active lattice vibrations do not appear explicitly in the formula. 
_ Arather important detail is omitted by this approximation, as the intensities are in 
fact inversely proportional to the square of the lattice frequencies (see equation 
)=(4.21)). 

In Miiller’s (1938) derivation of (5.1), which is based on Brillouin’s theory for 
_ liquids, the wave character, i.e. the time dependence of the thermal fluctuations, 
is taken into account. Here the individual frequencies of the lattice vibrations 
_ vanish by averaging over all the wave vectors. Hence (5.1) represents the total 
_ scattering of a thermally excited crystal and not the light scattered by a single 
_ elastic wave for a single angle of deflection betwen the incident and observed light. 
_ (This was the case investigated experimentally by Krishnan.) Finally, (5.1) does 
_ not contain the phase term «{°)Q,, so that it neglects a decisive contribution to the 
intensity of the Brillouin components. 


§6. THE BRILLOUIN COMPONENTS OF DIAMOND 
(i) Lattice Dynamics of Diamond 
Frequency and eigen-vectors of the elastic vibrations of diamond can be 


calculated from equation (2.17), which for cubic crystals has the form(s unit 
vector with direction of Q) 


2 
Ase E SoS. Se 
2 = 

Sas As,?—-E cass Wi OF eae... (6.1) 

S25 SS As?—E 
where 

2 
pv? —c ~Cy—c w p m : 

ait _ py, 1 # _ 4, y= — =velocity of sound, p= — =3-5 for diamond. 
Cyo + Cg Cyn t Cag |Q| q 


This equation can easily be solved for some special cases. For the case 


| OD, =9,=9: 
Ey 2=s,(A- 1), E. =$;(A+2), 
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For the case O,, 0, =O,=0 
E,,.=0,E,;=s, A, 
W,7=(0)150);Wo=(G01)), We=(1,0,.0)-. j 
For the case O,=0, O,, Q,, 
Ey 3=tAls,?+ 57 + {A(s,7—s,7) +45,787}"7], Fo=0. «..0.. (6.5) 
The eigen-vectors are found by substituting (6.5) in (6.1). 
A general solution for all wave vectors can be obtained with the help of a 


perturbation method used in quantum mechanics, starting from the fact that the 
isotropic case can be solved rigorously. In this case one has ¢41— C44 =Cy9 + C44, 


A=1. 
We put A =1+yand introduce the matrices 


ign SWE Royal ls Sat seer O90 
FLOR SES re aS ee Sy Sen lige add nl POL anes Seal) 9 ot tear (6.6) 
ee Sie ONO nes. 
Then (6.1) reads (HOLE y HD). Wi EW Nie earn eraser (6.7) 
For y=0 we have Oe ET a ge eee (6.8) 
with the solutions He ¢ polar coordinates of s) 
Eo =0, W, =(sings—cos¢, 0), eee (6.9) 
E,'=0, W, =(cos@cos¢, cos?sind,—sin@);/ «ss (6.10) 
E, =1, W, =(sin@cos¢,sinOsing,cos@)=s, «>... +a (6.11) 


which express that in the isotropic medium there are one strictly longitudinal and 
two transverse waves, i.e. 
Wo s— Wo s=W. Wo =—0, We) s—l. 8 ees (6.12) 


Equations (6.9) and (6.10) are not the only possible solutions. Each linear 
superposition x, W, + x, W, can be chosen as well. Ify40 we expand 


FE BWW OW rece (6.13) 
and obtain the first order equation 
(H®—E,) WY =(E,Y-HPWO, aan, (6.14) 
The condition for the solubility of (6.14) is 
fora=1,2: (a@-Ey) x4 cx —0, cx, 1 (b— Ea, = 0 ee (6.15) 
for a=3: (d—E,™)x, =0, 
where (HW, .W,) =4 sin70'sin#2¢6 =a, 
(A W, . W,) = 4 sin?26 (cos 46 + 14) =6, 6.16 
(OW Wz) =4 1 sin’ cos@sin4dd Sc, aul tt Oe 
(H™s.s) = ¢s1n46 (cos 4¢ + 6) + cos40 =d. 
The solutions are 
BY, = Mat b+ Vg}, By =, 
2c (a—b+/g)'? 
(*1)1,2 oe (42)"4(a—b + 4/g)'? ) (X2)19 = (4g. ccccee (6.17) 
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Equations (6.17) together with (6.9) and (6.10) gives the eigen-vectors of the 
zeroth order in y. These are still normal to s and describe a pair of transverse 

/ waves. Inthe higher approximation of (6.14) the pure transversality of the waves 

| generally vanishes. 

; Up to now all formulae in § 6(i) referred to the displacements W of zero order 

in Q, which are independent of the index k. For the calculation of the intensity 


Q 


of the Brillouin components, however, the relative differences acn(ar | 


introduced in (4.14) must be known. They have been given by Smith (1948) by the: 
following matrix equation: 


OMe eos Ws 
o) OAW, § 0. SOUWAW acs, (6.18) 


: Q,/W, Q./W, 0 


where £/« = 0-66 is the ratio of the force constants of diamond. 
Comparison between (6-18) and (4.14) yields 


Wee ep 20) FORD foe sey. FO FO eee (6.19) 
Allthe otherT are zero. 


i lhe determination of the eigen-frequencies and eigen-vectors for the optical 
~ branches is more complicated than for the elastic ones. Using the notation of 
Smith (1948) equation (2.16) for diamond becomes 
F-A J K A B C 
J -G—hk of B A E 
i ET 670, E A 
Azo BA". Cote baa” ani. K 
1S eee ts Se J/&G—-K a 
CePA A* K Ws eH | 


soi(12 


DE ome: (6.20) 


The A, B,... LZ are known functions of the wave vector Q. Equation (6.20) can 
be written in matrix form 


My,e,+ M,.e,=0, M*e,+ M,,65=9. Mans, ereLe (6.20 a) 


€1, €, are the column matrices which represent the displacements e,, e, of the: 
two atoms in the elementary cell of diamond and 


pas UCN amen a AE AEG 
Pe Gehl lett | BA Bele cas (6.21) 
ori GEA 


M,, and M,, are both symmetrical, hence 


My=My, My. = Mp, M,,.§My.= My. M,,. v++44. (6.22) 
Multiply (6.20 a) by M,,, then 
Mere ce Mines =O gels, Vey Mateo 0. sss (6.23) 
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Substituting for M,,e¢. = — Myje, and My,e, = — Mq,e, from (6.20 a) 
(M,- My, M#)e,=0 ) (M,2- Mi My.)e,=0. SL O.0-o0 (6.24) 


We have thus reduced equation (6.19) of the sixth order to two equations 
(6.24) of the third ordert. Each of them gives the eigen-values. ‘The eigen- 
vectors e, and e, cannot be completely determined from the homogeneous 
equations (6.24): two complex factors A, and A, remain undetermined, which can 
be found by inserting A,e, and A,e, in (6.19) and by making use of the orthogonality 
relations (2.15). 

The special case used in the following is that for which 0,=Q,=Q, 
and Q.r(k’)<1. The eigen-vectors are found to be 


AO==HOe (= é a =): — (6.25) 


(11) The Symmetry of Diamond . 
Notations used in text and figures are those of Smith (1948). Figure 1 shows 
the elementary cell of diamond. 


mo eueecseases 


Figure 1. First neighbours of the two carbon atoms in the unit cell of diamond. 


The values of the rectangular components of r ( ue) for the first neighbours of 
Oand O’ are givenin Table 1. 


Table 1 
~ 1 2 3 eee me eee es We or! 
xy 3a —3}a —}a 2a xy — a 3a 3a —}a 
Xa ta —ta ta —ta Xe —ta ta —tha fa 
x fa da — ha — 4a % ha -la da da 
a=lattice constant =1-78 x 10-8cm. It follows from Table 1 that 
ZQ-r()=0, SQ. rd) =0; = (6.26) 
l 


+ A similar reduction has been given already by Born (1914). 
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The symmetry operations of the diamond lattice can be described by six 
transformations together with a re-labelling of the lattice points (see Born and 


i Begbie 1947). 


4 Uaiagd Ba (1) (2,3) (3,4) (4,2) 
_ = Or Oral i ir a ae ee (6.27) 
Lg (1) (2,3) (3,4) (4,2) 
= 0G a 2 o 
T= 0-1 0 Bind) Su C4572 aes RD tC 2) ea Rn (6.28) 
0 O-1 
Ora 0 Soy 
T3= ttn 0 (Ay -(2)) G4) ml) 2) (33.4)) tee oe (6.29) 
;: OO Osea 
Disks 2 (3) (1,2) (2,4) (4,1) 
1 O .0-=1 a re ee ces (6.30) 
r0* 0 (3) (1,2) (2,4) (41) 
DOS ote (2) (1,4) (3,1) (4,3) 
T;= 0 O-!I Wee way fet ee 8) a Fees (6.31) 
ce ae Oat (2) (1,4) (3,1) (4,3) 
0-1 0; (4) (1,3) @ 1) (3,2) 
T= O-SrOuired gli h Aa 1) eee ee cee re (6.32) 
9 Wi Ua (4) (1,3) (2,1) (3,2) 
By applying the transformation matrices 7, lattice tensors of different rank are 
transformed in the following way : 
Xd A OL =A (irstramk)\—— oa ass (6.33) 
ex” >>: T,,X,,0rX’ =TXT  “(secendrank) {-  ...... (6.34) 
a B’ 
Xe py = UT yy T pg T yy X ops! Caution g51) | not ere a (6.35) 
a By’ 


Application of (6.27) to (6.32), (6.35) and (4.9) gives the results of Table 2. 


Table 2. 
l oO me es Ber 3h” -4- l l 1 Die eS a: 
O06, ju OO’ : oa, reveu 
71,1 0 A —A —A A Oy 2.1 = Mo3 3 = M21 1 = A320, 3 E =f) E =f 
O92 Oe AA — A %93,2— A311 —%g2 2 = 13.1 EO) ti = Eo Bs, 
O33 3 DPA A—A—A 31,3 — O12 2= 13 3— Mel, 2 Bt Eee 
11,2 “aa, 2 OVE eo. BB O12, 3— a1 3 a9) DS IO VDye 18, 
22,3 %11, 3 USED oa as, O23 1—%32 1 —42) DD VDD 
CEE eer al Oe Baar * B 31,27 413 2 De a 
i 
G03 1s 0’0) 4% found with the help of the relation 
l if 
Koo, OO’ => — Koo, O’O OT ee ee ORCC ca) (6.36) 
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(iii) The Frequency Shifts of the Brillouin Components of Diamond 
Corresponding to the experimental work of Krishnan we shall now treat the 
case that Q” is normal to Q’ and that Q has the direction of the cubic diagonal: 
O-=0,=0, )Fromi(3.2)iwennd 


[OQ [eos ee ots ioe) ee (6.37) 


r Ce 


vy is the frequency of the incident light =1-185 x 101° sec~1, cy the velocity of 
lightin vacuum; »=2-61 =index of refraction of diamond for A =2536°5a., 


OF=07=0,—|Q|/V3=5-2 710" cist ere (6.38) 


Substituting (6.37) in (6.3) and using the values of the elastic constants given 
by Bhagavantam and Bhimasenachar (1946) and used by Smith (1948) 


€4,=9°5 x 10" dyne/om?, ¢,,=3-9 x 104dyne/cem*;  ¢7—5-0 x 10 dyneicms 


we find for the frequencies (the experimental values are given in brackets) 
Yo =*+9.cm 45-4em—);. rg=9em- (8S cigs): ge eee (6.40) 


(iv) The Intensity of the Brillouin Components of Diamond 
To calculate the intensities one begins with the summation over J, k, k’ of 


equation (4.17) which reduces the 135 elements «,,,, ce) of Table 2 to the 27 


terms « Next, the summation over p, v, 7 is performed according to 


oo,u\ > 
(4.16), (4.17), which again reduces the number of terms to the nine elements of the 
first order polarizability «,, g . Finally, the summation (3.5) and multiplication 
with the factor C Gs Jo (4.20) or (4.21) leads to the intensities. 

In order to Bra the summation (4.17) it has been assumed that the changes 


of polarizability « ) of the atom k depend only on the displacements of the 


if 
eo aN 
first neighbours and of the atom kitself. (/,/=0,1,2,3,4). Then using Tables 1 
and 2, equation (6.19) and the relations 


i ei l L 
= P\e al = — =— 
r(O) =0, r(O j=H(a, a, a), le) = Gear Koo, u (ee) % oa, u Gal 


equation (4.17) for elastic vibrations reduces to (a=lattice constant) 
OS ge l l 
Xoo, i ie =12N a” O.+ 20,4 X oa, Lt O @; Xoo, n kk’ Pt Ole 2 


For the optical vibrations with si <1 and. e(O)=—e(0’, | 


fay 
vn 


A, (7 


?) =2,(4.18) becomes 


hie (2) =16N 
~ \Jo 


I 
aed ( 0’) (Ga Ae ee (6.42) 
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With the notation of Table 2, the eigen-vectors (6.3) and (6.25) and equation 
(6.19) the summation over y, v, 7 in (4.16) and (4.17) gives the following «,, 
E for the elastic waves (M =2iNQ,,//m): 
, longitudinal : 
044 = Mq9 = Ag = M(3a' + 2aA + 4aB)/4/3 =), 
dye = My = Hyg =X gy = Nog = Xgo = 4 Ma(E + 0-66D)/1/3 =e, 


transverse : 
(1) &y= — a3 =2Ma(A—B)/V/2=c, %33=0, 
Hyp =%q1 =O, O13 = %g1 = — Hyg = — Xgq =2Ma(E + 0-66D)/+/2=d, 


(2) 41 = Ogg = — $0433 = 2Ma(A— B)//6 =f, 
— 049 = FQ] = — yg = — %gy = — Ogg = — Age =2Ma(E + 0-66D)/+/6 Fi 


For the summation (3.5) the average values of the polarization vectors (3.8) are 
| needed. A convenient pair of unit vector s” and s’, which fulfil the condition (3.2) 
| and are orthogonal to each other is 


Mes 1 4/31 BR ZIT Seba any 4 
'- (55 5 a), Ss -(-, OT — 77%) + (644) 


Put A, A* pat $,8,+A,A* Ae $,5,=5 oo,n» (fOr p=o, =v the second term has to 
be omitted), then 


MSs; 3:=0-033, Soo-05 = 0-322, Sup 93 =0-167, | 


Sie, 12 = — 0-020, Sie, we) = 0-244, Soi, 22 = — 0-172, 
Soi, 23 =0-113, Si, 12 = 0-005, Soo, 23 = — 0-424, 
S41,13 =9°008, Soo 93 =0-077, Sy 13 = — 0-031, | (6.45) 
Soo, 11 = 0-456, S33, o2 = 0345, Sai, 33 = 0-012, aed ra 
cE S42,31= — 9-034, Sor 12 = — 0-020, Soo 13 =0-114, | 
Py S13,52= — 0-038, Sor a3 = — 0-083, S31 29 = 0-024, | 
Soe, 33 =0°156, Sgg 11 = 0-489, Si. 03 =0°113. | 
All the other S,, ,,, are zero because they contain the vanishing factor s,’ =0. 


With these values one obtains according to (3.5) (the factors (w+w,,,)* are 
neglected): 


I¢p,) = (0-5126? — 0-670be + 1534e2)C aA | 
Tp, tx) = (0-375¢? + 1-023f? —0-498cd + 0-182fg+1-085d? wae. (6.46) 
440 92 Q 
+2-44992)C ts 2) . | 


The indices (B, 1) and (B, tr) denote the longitudinal and transverse Brillouin 
components. 

A, B, (E+0-66D) and « can be derived from the index of refraction and the 
elasto-optical constants of diamond, which have been determined byRamachandran 
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(1950). The method of derivation will be given in a following paper. The 
results found are as follows (R=n*/327Na, n=index of refraction, a=lattice 
constant, VN =number of unit cells per cm’) : 

A=Rpi,;—= —182 x 10-* cm?) (B=hp,, = 053: 10s cm. 


0) eyed 


: eel it). iG aren eee = Si 
(E+0-66D) =Rpy= —0-71 x 10-% em’, — = 5 >= = 1-47 x 10% em’, 


egies (6.47) 
where p,,=—0°31, pj5=0-09, p,,=—0:12, »=2-61. ‘This yields 
ag x 10-¥, e= — 30-8 ee 210 malip2 ie pel Goes 
a/m 4/m Vm 
26360 ip ge 21333 2 ip eee eos 
BT ers oatl ay cll lee Bei pa v= LS a : 
Krric (6.48) 


Substitute these values in (6.46) and multiply by 


Q 1 —26 Q 1 eb 
c(¥)- = gong, N75 x10 and C(g°..) = gaye #75 x10 


hen the intensities at room temperature (300° K.) are 


N N 
Ip, =4576 =— Be x 10°, Tip try) = 14050 Dae 5 eer oes (6.49) 
The intensity ratio between transverse and longitudinal Brillouin components 
is 
Law /lepHs2lee So Oe eee (6.50) 


If the phase term O90 On in (6.41) were neglected one would obtain the result 
Lg, t)/Le,1 = 64 which is in disagreement with experiment, as will be shown. 

The only available experimental data are communicated ina paper by Krishnan 
(1947). He gives two microphotometer records, which are reproduced in Figure 2. 


(a) (b) 


Figure 2. Microphotometer record of Brillouin scattering in diamond for two different orientations. 
L.B.C. and T.B.C.=longitudinal and transverse components. 


Because of the considerable width of the lines a quantitative interpretation of 
the photometer record is rather difficult. From Figure 2(a) it is possible to 
conclude that the ratio Ig 4,)/I(g,1 lies between the limits 


Phe ieee Pen (Si 
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_ Figure 2(b) does not allow a quantitative measurement of the intensity ratio. 
| Krishnan only points out that the transverse components are much brighter than 


the longitudinal. In any case, the calculated ratio 3-27 is in good agreement with 


f the experimental result, as given in (6.51). 


The intensity of the first order Raman line (Av =1332 cm!) depends only on 
the polarizability derivatives «,,.,(p 404m) which cannot be derived from the 
; photo-elastic constants. This different behaviour comes from the different 

_ properties of the terms Q.r ( oa and A, (x8 = ) which, according to (4.17) and 


io) 


_ (4.18), are characteristic for the intensities of the Brillouin components and the 


optical Raman lines. The sign of A, (x | and of the polarizability derivative 


| D of Table 2 is independent of /, whereas the sign of Q.r depends on / in 


i 
- kk’ 
| the same way as the polarizability derivatives A, B, E. Hence the products 


l nee, f q l 
oe te ) Au (x ®) cancel in avs fore 55%; (ee) =A, B, E; the products 
Reo! i Q.r cancel for «,, ,, G3 =i 
To derive the reer of the Raman line one can use the fact that in the case 
_ treated here (O,=Q,=Q.) the eigen-vectors are the same for the elastic and 
_ optical lattice vibrations. The same holds, of course, for the averages S,, ,,,» 
| Hence, the expression for the intensity must have the form given by (6.43) and 


(6.46). The only difference is that only the polarizability derivatives«,, ,, ( BE) = D 


_withpocremain. Consequently A and B, whichare of the form «,, ,, ( co must 


be omitted in (6.43) and instead of the term 2 Ma(E + 0-66D) one has to introduce 
the term 16ND/./m which is characteristic for the optical vibrations (see equation 
(6.42)). This gives with 


Q as 1 5 —30 
c( = gaye, 197 x10- from (4.13), 


ies Ny 2 — 26 
Iipy=9 a x 10 e 


ee .« (6.52) 
| 
J 

From Krishnan’s observation that the Raman line at room temperature is at 
least four times stronger than the total intensity of all Brillouin components, it 
follows that D=a,, ,(p AoA) is at least 0-64 x 10-1 cm?. ‘This is of the same 
order of magnitude as the other polarizability derivatives A, B, EF and gives a 
further support to the present theory. 
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ABSTRACT. The Bloch theory of electrical conductivity of metals assumes that the 
phonon distribution does not deviate from equilibrium. Peierls has criticized this assump- 
tion, showing that in the presence of an electric field a stationary state is achieved only through 
processes which do not conserve the total wave vector. These can be Umklapprozesse 
involving electrons and a phonon, or phonon—phonon interaction of the type responsible for 
thermal resistance in dielectric solids. He estimates that if only the latter were to act, 
strong deviations from Bloch’s formula would occur below 50° k. He therefore concludes 
that the electrical resistance at low temperatures is determined by Umklapprozesse. 

A more detailed treatment of the statistical equilibrium of phonons and electrons shows 
that the relaxation times of the important phonons is smaller than estimated by Peierls by a 
factor of about 100; consequently one would not, as a result of phonon-phonon interactions, 
expect important deviations from the Bloch theory above 10° x. It is thus unnecessary to 
consider Umklapprozesse. Since Umklapprozesse require a special shape of the Fermi 
surface in wave-vector space, Peierls has concluded from the apparent agreement with 
Bloch’s theory that the Fermi surface has this special shape. Hence it is now shown that 
this shape is not necessary. 


Sil JON MPRO ID OCA TON) 
VER since Peierls (1930, 1932, see also Sommerfeld and Bethe 1933) had 
pointed out that in the presence of an electric field interactions between 
phonons and electrons conserving the sum of the wave vectors cannot 
establish equilibrium, the agreement between Bloch’s theory and the observed 
electrical conductivities at low temperatures has remained unexplained. 

In order that equilibrium be obtained, there must be processes which do 
not conserve the wave vector: these could be Umklapprozesse or phonon—phonon 
processes. ‘The last-mentioned processes are the same as are responsible for the 
thermal resistance in non-metallic solids. Peierls deduced their effect from the 
known thermal conductivity of a typical non-metallic solid, NaCl, and came to 
the conclusion that these processes were sufficient to establish equilibrium at 
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room temperatures. At low temperatures, however, they would be so rare that 


| the phonon distribution would depart from equilibrium sufficiently to cause 


_ amarked increase in the electrical conductivity over that given by Bloch’s theory. 
| The Bloch theory would thus hold only to about 50°x. 

Since the electrical conductivity seems to agree with Bloch’s theory even at 
lower temperatures, Peierls considered the Umklapprozesse. ‘These are 
interactions between electrons and phonons in which the total wave vector is 
_ changed by an inverse lattice vector. The resistance due to these Umklapprozesse 


_ was estimated, and has the same temperature dependence as the resistance 
_ calculated by the Bloch theory. There are, however, two difficulties. Firstly 


the absolute value of the low-temperature conductivities are related to the 
high-temperature conductivities in the way predicted by Bloch’s theory, while 


no such estimation is possible for Umklapprozesse. Secondly, in order that 


Umklapprozesse can be effective in producing resistance, the energy surface in 
wave-vector space corresponding to the Fermi energy must come up to the 
surface of the first Brillouin zone, whereas a number of metals conform to the 
Bloch theory in spite of having Fermi energy surfaces well inside the first 
Brillouin zone. 

It thus seems probable that Bloch’s theory has a wider validity than is 
conceded in Peierls’ argument. It will be shown here that his estimate of the 
sufficiency of phonon-phonon interactions to maintain the phonon distribution 
near equilibrium was too pessimistic, and that in fact the Bloch theory is valid 
at all low temperatures. The two difficulties mentioned above are thus resolved. 


§2. RESUME OF BLOCH’S AND PEIERLS’ THEORIES* 

The thermal motion of the lattice destroys its periodicity, causing transitions 
in the state of the electron gas. Resolving the thermal motion into quantized 
lattice waves of wave number q and energy fw, the change of state due to 
interaction takes place as discontinuous processes in which an electron of wave 
number k’ is emitted and an electron k and a phonon q are absorbed or vice versa, 
so that the interference condition 


k’=k+q ...... (ia) @ ormikisks- qe lag, ‘Veo. ee% (1 5) 
is fulfilled. Here g is a vector of the reciprocal lattice, and has to be included 
in (1 5) since we specify an electron by its ‘reduced’ wave vector only. In addition 
we require for non-virtual transitions that energy be conserved, i.e. 

E(k) =£(K)+he ©  —aeaees (2) 


If N(q) is the number of q-phonons and /(k) the average occupation of the 
electron states about the wave number k, then the rate of change of the electron 
distribution due to interactions with phonons is given on a simplified model as 


af] 206-8 ag | 


x[{1—f(K)} (K+ gt N(q@) + 13 — fk) {1 fk + a)}.N(Q)] 
+ QUE(k + q) — E(k) +heo}[ {1 —f(k)} fk + q)N(—q) 
fk) fk+QMHM-M+U see (3) 


* See Sommerfeld and Bethe (1933). 
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where C is a coupling constant, G? is the number of ions per unit volume and. 
M is the mass of each ion.’ The resonance factor Q is defined by 


_ Q(x) = (2h? /x?)(1 — cos xt/h) 


and ensures conservation of energy. 

Correspondingly the rate of change in the number of phonons through 
interactions with electrons is 
aN ZEEG G0 
ale = “OM c 5p {E(k + 4) — E(k) —hw} | 

xf fk) AK + g)iN(q)+ AK) f+ Qi N(Q@)). (4) 

Both (3) and (4) disregard processes (1), i.e. Umklapprozesse, which only occur 
near the top of an energy band and need special treatment (Peierls 1932). 

It is easily verified that Of/dt], and 0N/dt], both vanish if we substitute for 
f and N the equilibrium distributions 


awe tl)” AH 1), eee (5) 
where «=(E—)/RT, x =hw/RT and Cis the Fermi energy, and if we take account ~ 
of the conservation of energy. ‘Thus the assembly of fermions and bosons. 
mutually interacting is in equilibrium when each component itself is in 
equilibrium. 

If the total wave vector is conserved in addition to the energy, the following 
distributions are also stationary if |A| is small 


jo =[exp {((2-—C€+4.k)/RT}$+1} ] 
=fot one fy if |A| is small 
kT Oc | (5’) 
AN” =lexpioreAaa)/Ri\= pa Gen tiles ae eee 
Be EON 
HN + Fa a> if [A] is small 


where A is a constant vector, and must be the same for fy’ and VW’. 

Hence, while it is true that an assembly of electrons and phonons in mutual 
interaction will generally tend to the equilibrium state (5) in the absence of 
external forces, when all interactions conserve the total wave number of the 
system J=>,/f(k)k+2%,N(q)q the assembly will tend to a quasi-equilibrium 
given by (5’). Here A is related to J by 


ae 


) 
J- 3kT 


a{2 ks Teak +3 |" Fe aap. hoe (6) 


obtained by substituting (5’) into the expression for J and summing. The fact 
that true equilibrium cannot be obtained solely by processes conserving the 
wave vector was first pointed out by Peierls (1930). 

Bloch’s theory of electrical conduction is based on the assumption that the 
phonon distribution is in true equilibrium, i.e. N=, and that all processes 
conserve the total wave vector. Then the stationary distribution f in the presence 
of a field F is given by 


oa ci ICE en | a ie zeae (7) 


at 1% dh ok So 
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To obtain Of/dt],, treat q as continuous, replace N by -/ in (3) and integrate 
_ over the angular variation of q. Taking 


Waele ay meet ie eat c Goeae ae (8) 


one obtains the following equation to first order in y: 

i qo SAO hw q 

. 2 0 ote) CS 
| tw) ee (ate ayn f * kdE|dk sal 


py l aliers hog? | 
a: ee (hw/RT) E ~ kdEydk ~ aA x(E—hw) 


[o(E +ho) fo(E —he) 9nMveF 1 df, (dE\? 
My om CalkT) + | anh = Fae 5 a8 


where v is the phonon phase velocity, Q) the volume of a unit cell and g, the Debye 
wave number. The solution of (9) gives us y(£), hence the electron distribution, 
and the current 


i ~§¢ | Gx) a a ree (10) 


whence the conductivity is found. 
If T>0, the Debye temperature, all terms on the left-hand side of (9) cancel 
except the terms in g?/2k?, and the solution of (9) is 


BG PELE Li ee (11) 
M k (dE 0 
where = he C2 (=) Qokbel 7, susie otets (12) 


whence the conductivity becomes 


_nyM (RdB)* kb 6 
Cdk} hka,T 
Here n, is the number of conduction electrons per atom, a) the Bohr radius and 


the barred quantities k etc. refer to the Fermi wave number. 
On the other hand if T<0 


| 

Mee is Ee Va cial ah Vance (14) 

where it can be shown that f(E) <x(E£) a 
Qn MveFk (dE\? ( 0 \> 
ya See TBO an) 

pO yt 

where J, = ie emp and hence ge 

ag hE M k dE “a (z) 133) 
EN TI ELS) 5 SAY eal ia ae ( 


According to Bloch’s theory if T,«<@<« T, from (13) and (13’): 


ao = 497: “(3 en sks (14’) 
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and this relation between the conductivity at high and at low temperatures seems 
well verified by observations. 

The temperature dependence can be qualitatively understood as follows: 
An electron moves in k-space in discontinuous steps due to its interaction 
with phonons. Consider the action of phonons in a wave number interval q, dq. 
The rate of change of the distribution function is characterized by a relaxation 
time 7,, which is proportional to the time required for an electron to move a 
given distance over the Fermi surface. The movement can be regarded as 
a random walk over the Fermi surface and the number of collisions required to 
move this distance is inversely proportional to the square distance of each step q?. 
We have 

of 1 __ number of steps per unit time 
Of lan te * ~humber of steps required 


and since the number of steps per unit time is proportional to (q) q?dq, the 
number of phonons in the interval g, dg, we have 


Of[Otlage NV (q)qtdq na ees (16) 
so that Z| x { AW(g)g-d¢d; | «=. =) = 8 nee (17) 


Hence if T7<0, the resistance is proportional to T°, while if T>0O, it is 
proportional to 7. 

We have seen that there are two assumptions in the Bloch theory, i.e. that 
N=WN and that all processes conserve the total wave vector. These are 
incompatible, as was already pointed out by Peierls, for the second assumption 
implies a continuous increase of J due to the electric field 


Ol Ot lpia = CF 1 ed Whee ee (18) — 


which, in the absence of Umklapprozesse (1b), requires that the phonon distribution 
deviates from equilibrium. Also if the phonon distribution were at equilibrium, 
the phonon-phonon processes, which are the only other processes capable of 
producing a stationary state, could not be effective, since 0 /dt],=0 by 
definition of .1. 

Hence if we disregard Umklapprozesse, we cannot assume that N=./, but 
must calculate the deviations of the phonon distribution from equilibrium, and 
then determine the effect of this on the electrical conductivity. This has been 
done by Peierls in the following way: 

The rate of increase of J due to the field (18) is balanced by the rate of decrease 
due to phonon processes, so that 


21S naseatgeeX NG hs tcc ie be 

ate ae al ei ot He reas oy, 
where Jg=2,q/(q). Peierls now assumes an overall relaxation time for the 
decrease of Jg, i.e. al le (20) 


and he further assumes that the relaxation time is the same as occurs in the 
expression for the heat conductivity of a lattice which does not interact with 


free electrons 
Kg Sere” | ee ee ee ee (21) 
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| since the processes not conserving the wave vector are the same in both cases. 
_ Here S is the specific heat per unit volume of the lattice. 


Now if 
F=fot i (22) and N=NV + 2A eee (22’) 
then, since 0f/dt],,=0 for A=, we can put 


al re (“=*) 2 M5 eee (23) 


where df/dt},, is the rate of change of f due to phonons if N=./, and dy is the 
corresponding deviation of the electron distribution. Since the rate of change 
due to the field is independent of the phonon distribution 


-F 1 nal Gates (23”) 
GE Ol, Ne Olay Lobe eae 
| so that A=Az+p. But the electrical conductivity is proportional to A, whence 
Och lis) gy bw 8 Pt Cake (24) 


where og is the conductivity according to the Bloch theory. 
Now vp is related to Jg by 


ree Ny aaa ae 
Je=2N@a=n (7) do at(z) win 5010-0 0.6 (25) 
obtained by substituting (22’) for N, while 
J ge OLNT EW = SENKe EF TORT, ww wie es (26) 


which is obtained from (19), using (18), (20) and (21). Thus » can be obtained 
from the lattice (thermal) conductivity. Similarly Ap can be obtained from the 
electrical conductivity using (12) and (13) or (12’) and (13’). In this way Peierls 
finally obtains for the ratio «/Ag, which is indicative of the deviation from 
Bloch’s theory : 


6kgen k® (0\4 1 = 
pis = = (7) ie So sodn (27) 
When T>6, p/Ap “kg, more generally 
Kk (0\41 2 J; 
wihnocx (7) p(s) Te” REsleiopete (28) 
and when T<6, B/Ay cK(O/T)?. 


Evaluating y/A, numerically from (27), there is an uncertainty about K,: 
the quantity which is required is the thermal conductivity of the lattice, if the 
lattice did not interact with the free electrons. To get the order of magnitude, 
Peierls uses the thermal conductivity of a typical dielectric solid, i.e. NaCl. 
He obtains p/Ap =9-21 at 290°K. But since xg increases as the temperature is 
lowered, Peierls estimates that 1/Ap =1 at about 50°k., and that phonon—phonon 
interactions are therefore not able to prevent appreciable deviations from the 
Bloch theory at such low temperatures. 
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§3. EQUILIBRIUM BETWEEN PHONONS AND ELECTRONS 
It is possible to criticize Peierls’ argument, since the relaxation time tg 
defined by (20) is not necessarily the same as 7, defined by (21). Neither equation 
takes account of the frequency dependence of the relaxation time. Equation (21) 
should be replaced by 3 


K=42 BG Eg. l) 07 7(q)d¢ 2 a) eee (29) 
J 


when E(q, T); dq is the energy content of modes of wave. number q, dg and 
polarization 7; and 7,(q) is the effective relaxation time of the (q,7) modes, 


defined by 
oN ey) LAGEN ENG) Tae Bee (30) 
dt Ip 7;(4) 
Since in the presence of a temperature gradient 
ON aN dN 
aeeeh acta —(v. | alee oe 
ot [ Ot _laritt Neate, aT BD 


the relaxation time is simply related to the deviation from equilibrium by 

hw exp (hw/RT) 

RT {exp (hw/RT)—1?”” 

Cer ee (et — 1) ee ee (32) 


This effective relaxation time depends not only on w and T, but also on the 
effective relaxation times (i.e. deviations from equilibrium) of all the other 
modes, as has been pointed out elsewhere (Klemens 1950, 1951). It is thus 
clear that to a certain degree + depends upon the frequency dependence of those 
processes which cause the deviation from equilibrium. This is an additional 
reason why one cannot identify 7, of equations (20) and (21), though, in fact, 
this will be shown to be less important than the frequency dependence of r. 

In the case of a temperature gradient in a dielectric crystal, when 
three-phonon processes conserving the wave vector have a negligibly small 
relaxation time, one can show, as a result of the frequency dependence of N— WV 
expressed in (32), that for all polarizations 


7(q)=Tu (RT/iv) for q<kT/hv 
=71 (q) for, g2RE/o; eee (33) 


when 7;;'(qg) is the relaxation time of transverse phonons of wave number g due 
to processes not conserving the total wave vector. 

In the case of electrical conductivity, the processes causing the deviation 
from equilibrium have a different frequency dependence. We obtain from (9), 
after integrating over all k in a manner outlined by Sommerfeld and Bethe (1933), 
that the rate of change of the phonon distribution due to interactions with 
electrons is 


N-—WN =(v.grad T) 


ONG een G A Bok ge" 
ot |. 9%nMkT (dE/dk): (e* — 1)? 


so that, since ON/ot| FON /Otl,=05) 2) see eee (35) 
we obtain from (30) 


so (Ens (34) 
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Consequently in the case of electrical conductivity, equation (33) has to be 


| replaced by the following, as is shown in the Appendix: 


Rie? /RT ; 
7;(q) = hog © Fe) for q<kT/hv 
=T1 (9) for g>kRT/hv. Reed Cs 


) Using (22’) for N—.W, and substituting into (34), we find, using (36) and 
| dividing both sides by xe*/(e"—1)?, that 


: Wee AMUATIN-= 38) 9 Bie ie ky ten es (38) 

or pe = ANVATHN I +heAtTe), -— ha wee (38’) 
2G-3 he 

fe where ee eee, Me. see (39) 


~ On Miv? (dE/dk)> 


Thus » depends upon frequency, since + depends upon frequency and 
x=hw/RT, while i, the deviation of the electron distribution, is the same for all 
phonons. To obtain the deviation from the Bloch theory, we must find 2. 
As before o/og =A/Ag. As in (23’) we consider the rate of change of the electron 


| distribution in the case when »=0 and A=A,, and in the actual case. But 


since . depends upon the phonon frequency, we must consider the rate of 
| change of f due to phonons of the wave number interval g, dq, and then integrate. 


DE sss EE A pes Sot I pe Cte 

131.27 a a= | 5 ae ieee (40) 
Now 0f/0t},,, the rate of change due to phonons gq, dq when p =0, is given by (16). 
Hence oy 

| Af SW (aa! d= idan sik mii? (41) 


and expressing yu in terms of 4 according to (38’) 
ian ay, 


9 e*—1° 


Op Ap 


(oo) uy 
“BB | re as 
J0 = 


This expresses in general the deviation of the electrical conductivity from the 
value given by Bloch’s theory in terms of 7; this in turn is related to 7,;’, the 
relaxation time of phonons due to processes not conserving the wave vector, 
by (37). 

In general o/oy as calculated by (42) will be smaller than the same quantity 
given by Peierls’ argument, 1.e. by (24) and (27). This is because the frequencies 
most important for the conduction of heat by phonons are not the same as those 
important for establishing electrical resistance. In general, higher frequencies 
are important in the latter case, and they have a shorter relaxation time. 


§4. THE RANGE OF VALIDITY OF BLOCH’S THEORY 
We shall now apply these results to the case treated by Peierls, i.e. we shall 
find the temperature for which o/c, =2 for the case of gold, assuming the lattice 
properties of a typical alkali halide, namely KCl. 
The thermal resistance of alkali halides in the liquid hydrogen region is due 
to lattice defects, which scatter phonons according to a fourth-power law. 
De Haas and Biermasz (1937, 1938) have found that in this region 


Rea Pan WAtt CIN OCB. oy kw ee (43) 
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whence from (29) with (33) one obtains, as shown by Klemens (1950, 1951), 
Ti OTs LON NSEC, fe eee (44) — 

Put ty =79(xT)-* then from (42), 


g,/o'= | Papeete oe (45) 


0 e*—1 


2CG3 BR kr 


e erie beet 
where a =hvAr)T4 = nMe Gmdee TT? (46) 
(es) xt 
Now op =$o0 when { ——= [1 +axn*]-*dx=13, 
9 e*—1 


i.e. when « =105, as has been found by numerical evaluation. 
We thus get for the critical temperature, below which the deviation from 
Bloch’s theory is too large, 
2C?G* R? 


6 ye = SS 
T= ae poa Te (47) 


; 3 2m k®R?0G-* ] 
Using (13), we get (i= Taal os Pantene T a) : 
Now for gold-k=1-21 x 10®cm-4, -G?=5:95 x 10” cm. 059) =5°76 x 10” Esra 
a =5:28 x 10-°cm. (Bohr radius), v=2 x 105cm/sec., 0=175°K. and we take 
my=1. One then obtains T,=10-8°K. We have neglected the difference between 
7,(q) and 7y'(q) for g<kT/hv; but this difference has very little effect in (42), 
since the interval «<1 contributes very little to the integral (45). 

We have thus seen that the deviations from Bloch’s theory due to the finite 
relaxation time of the phonons are less than estimated by Peierls and become 
large only at about 11°x. At lower temperatures, however, the size of the 
specimen will limit the phonon relaxation time (Klemens 1950, 1951), and this 
will again limit the deviation of the electrical conductivity from Bloch’s theory. 
Thus the deviation of the conductivity from Bloch’s theory will only be 
appreciable, if at all, in a small temperature interval just below 11°k. 

It is of course not justified to use the properties of the KCl lattice for a metal, 
and both Peierls’ and the present calculations can be criticized for this reason. 
Nevertheless the present calculations have shown that in general Peierls’ estimate 
is too pessimistic, and it becomes plausible that the observed agreement between 
Bloch’s theory and observation at low temperatures should be ascribed to 
interactions amongst the phonons rather than to Umklapprozesse. This would 
resolve the difficulty referred to in the Introduction. 
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APPENDIX 


The proof of (37) is very similar to the proof of (33) in the case of thermal 
conduction (Klemens 1950, 1951): it can be shown quite similarly in the case of 


WO1=W, GtG=Go eaves (A 1) 


@ have no appreciable effect on N(w) unless |w,|>]w| and therefore |w,|> |w. 
In that case 


= | =D Aye NN (Ne +1)-(N+1(N+1)NQ] «0... (A2) 


} and putting N=./ +7 etc. the factor in the bracket can be written 


(Ny —N 9) +N (ny — Mg) — (1yVg +123); 


neglecting 7,—n, this is approximately 


no 3 
sid Pin aa becte PPA my Ua helices aes 78k A3 
28i5, =o 1 1 ( ) 
The largest contribution to the rate of change comes from frequencies 
w,=kT/h; thus if w>kRT/h, the three-phonon processes do not contribute 
appreciably to the rate of change since interactions with w,;<w are not important. 


~ Only when w<kT/h will such interactions be important. In this case the ratio 


of the two terms in (45) becomes 


ON, wo __ (hw)? exp (w/RT) _{exp(ta/RT)—1} ro 
a ae ar fexp (hw/RT)—1? 7,(ho,/RTP 40,” 
Raat (A4) 


having substituted (36) for m and m,. Here 7 and 7, are the effective relaxation 
times for frequencies w and w,;. We now remember that w,=RT/h for maximum 
effect and obtain the effective relaxation time in terms of +’ and o,, where 


ON 
zi = Gece a 
i.e. o, is the relaxation time of the three- acl processes if all modes are at 
equilibrium except for w mode. We then get 
7'(w)[o,(w) + (RT /hw)r’(RT/hw,)] 
7(a) = ee + 
T'(w) + 04(w) 


and in the case of crystals, when o, <7’, we obtain (37). 
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A Combinatorial Problem in Counting Cosmic Rays 


By ERWIN SCHRODINGER 
Dublin Institute for Advanced Studies 


MS. received 1st August 1951. 


ABSTRACT. The probability is computed that precisely m out of ry counters be hit by 
at least one out of m cosmic particles, each of which has the probability 1/r of hitting 
any particular counter. The estimation of m from given 7 and m is discussed. 


cosmic-ray shower and hit by n(>1) rays. A particular ray hits one, and 

only one, counter and the probability of its hitting any particular counter 

is 1/r. What is the probability of precisely m counters going off? Thus m is 

at least one and at most 7 or n, whichever is the smaller. The ultimate aim is, 

of course, to estimate 7 from the observed m, but this is a second question, to be 
touched upon at the end. 

Let p(n, m) be the probability we want to find. ‘The clue is the recurrence 

formula 


Tes problem was put to me: A tray of 7 equal counters is exposed to a 


p(n+1,m)=™ p(n, m)+ romt* p(n,m-1), iden (1) 


Indeed, the additional (n+ 1)th ray can produce an m-case either from an m-case 
by hitting one of the m counters that are already hit, or from an (m—1)-case by 
hitting one of those that are not yet hit; and in no other way. By reflecting and 
trying one is led to the assumption (which will have to be proved) 


p(n, m) = 6b", “ae mS =, (Sy ey een (2) 
where the b’s are constants and ] 
fay (X) =ale—1)(4—2)<. ns. (Se ee (3) 
If we carry (2) into (1) we get 
beet =m bt rr (4) 


Now from common sense p(1, 1) =1, and p(1,m) =0 for m1. Thus (2) is correct 
for n=1, provided that we take 5';=1 and 5!,,=0 for m1; it will therefore 
remain correct for any n, if the b’s are successively determined by (4). This can 
conveniently be done by building up (or rather ‘ building down’) a number 
triangle, reminiscent of the schoolboy’s method for finding the binomial 
coefficients : 


1 

1 1 

1 3 1 

1 7 6 1 

it L5 25 10 Penk steep ts aa be 
1 31 90 65 1s 1 


< 
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uf (The mth number in any line is obtained on multiplying by m the number above: 
/ it and adding the latter’s left neighbour; empty places are considered to be zero.) 
This solves the practical problem. The 5-coefficients are familiar to the 


{| statistician. The sum of the probabilities (2) for fixed m and m=1, 2,3,....n must 
be unity, and that for any integer 7. It follows that the relation 
= Eb pty (2) 


m=1 

| holds for any integral value of x, and is therefore an identity. In statistics it is. 
} mainly used for expressing the ‘moments’ by the ‘factorial moments’. The 
» explicit expression for 5”,, 


Brg E (Den = ay E (— (may 


MMe aN Taya ae Ns Sita, aN (6) 
igi, Lam ant 
One easily confirms that it obeys (4) and gives the correct initial 64,=1. I[ hit 


~ on (6) ina roundabout way, not worth being reported. It is bound to be contained 
| as a special case in a very complicated formula quoted by Kendall (1948, p. 58, 
equation (3.20)). 

Let us now briefly review the question of estzmating the number of incident 
i rays n from the observed number m of counters that have gone off (while r—m did 
not). I believe the only reasonable method of estimation is here as in many other 
cases that of the ‘median’. Let the observed number m be m*. We would 
regard as the most likely value of m the one that gives equal probabilities to 
m<m* and tom>m*. One must have prepared, from (2), ‘for every n’ not only 
the series p(n, 1), p(n,2), . .. . p(n,r) (of which the sum is 1), but also the sums 


Se) (Gases 
m= 1 


up to that s for which the sum first exceeds 3. By scanning these tables one will 
_ find which » most nearly fulfils the condition of the median. Border cases, and 
_ particularly cases with small numbers, are baffling. For example, let m*=1. 
meshave Al t)=—1, p(2,1)=1)r; p(2,2)=@— Dir- ~Here n=1 fulfils the 
condition of the median exactly. If 7, the total number of counters, is not too 
_ small, this is a safe estimate. But if for instance 7 were only 2, then n=2 gives 
still a 50:50 chance to m=1, and the estimate n =1, though the best, is poor. 

The most embarrassing case is m*=r. ‘The median is then reached only in 
the limit noo. Obviously this case yields no reasonable estimate. 


My thanks are due to Professor C. B. A. McCusker of the School of Cosmic 
Physics who raised this problem. 
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Observations of the Penetrating Non-Ionizing Component of 
the Cosmic Radiation 


By J. C. BARTON 
Birkbeck College, University of London 


Communicated by E. P. George; MS. received 1st June 1951 


ABSTRACT. An anticoincidence arrangement of Geiger—Miiller counters was used to 
detect that part of the non-ionizing component which is capable of penetrating 10 cm. of 
lead. Measurements were made at sea level and at an altitude of 3,560 metres. The mean 
absorption length of the non-ionizing particles was found te be 170+5 gm. cm~ in air, 
190+35 gm. cm~? in carbon and 400+ 25 gm.cm~ in lead. The secondary particles 
were strongly absorbed in a few centimetres of lead. Studies were made of the angular 
distribution of the secondaries and of the production of two or more secondary particles in a 
single process. Less than 0-5% of the non-ionizing particles were found to be associated with 
extensive air showers, but about 10% appeared to be produced by ionizing particles inter-_ 
acting in the lead above the apparatus. 

The angular distribution of the secondary particles produced in carbon was approximately 
proportional to cos 0. 

The nature of the non-ionizing particles is discussed and it is concluded that the only 
known particles which could explain the results are high energy neutrons. On this hypo- 
thesis, their energy spectrum and intensity at the two altitudes are calculated. 


S17 INDRODUCTLION 

T the present time several types of non-ionizing particle are known to be 
present in the cosmic radiation. These include (i) photons, (ii) thermal 
energy neutrons, (ili) neutral mesons, (iv) the penetrating non-ionizing 
component discovered by Janossy and Rochester (1943), and (v) the particles 
responsible for the majority of the nuclear disintegrations observed in 
photographic plates. It is now generally believed that the same particles, high 
energy neutrons, are responsible for both (iv) and (v), but there has remained a 
discrepancy between the intensity of neutrons required to account for the two 
sets of results. Whereas Janossy and Rochester concluded that the intensity of 
neutrons at sea level was 0-14°% of the total intensity of charged particles, the 
equivalent figure given by George and Jason (1949) as a result of photographic 
plate work was 6%. Since, in addition, an alternative explanation for the stars, 
in terms of unstable neutral particles, had been put forward (Perkins 1947), it 
seemed important to repeat and extend the measurements of Janossy and 
Rochester. The present work has been carried out for this purpose and, more 

generally, to obtain further information about the non-ionizing component. 


§2. EXPERIMENTAL METHOD 
It was shown in a series of experiments by Janossy and Rossi (1940) and 
Rossi et al. (1940 b) that the only satisfactory method of measuring the non- 
ionizing component with Geiger—Miiller counters would be to use an efficient 
anticoincidence arrangement. This was first achieved by Janossy and Rochester 
(1943). The general arrangement of the present apparatus, Figure 1, is rather 
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similar to that used by the latter workers, although the geometry is slightly 


| different. 


The apparatus consisted of a counter telescope B, C, D, surrounded on all 
| sides except underneath by a screen of anticoincidence counters A. When 
the telescope recorded a threefold coincidence without any of the anti- 
coincidence counters being discharged, the event was interpreted as being caused 
by a neutral particle entering the apparatus and producing a secondary ionizing 
particle in the material at b. This type of event is usually referred to as an 


_— anticoincidence. The question of spurious events will be discussed later. It 


is useful to note here that any effects due to photons were greatly reduced by the 
lead screening a, which had a thickness of 10 cm. on top and 5 cm. at the sides. 
The Geiger—Miiller counters used in this work were of the normal copper foil 
cathode-glass envelope type. The counters in the anticoincidence screen were 
used with univibrator circuits which were designed to reduce the dead time of the 
counters (Simpson 1944). This enabled the ‘ blocking’ time of the apparatus 


Longitudinal Section “Tom, Transverse Section 


Figure 1. Arrangement of counters. 


(cf. Janossy and Rochester 1943), 1.e. the time during which the apparatus was 
made insensitive to the recording of true coincidences, to be less than 1%. 

Most of the circuits used were of a conventional type and will not be described. 
The actual anticoincidence circuit was derived from that used by Hodson (1948). 
An additional stage was added to avoid troubles arising from the fact that a 
univibrator circuit is insensitive to a pulse for a short period after returning to 
its stable state. 

Three separate sets of observations were made. ‘The first and third were 
carried out, during the periods January to June 1949 and December 1949 to 
February 1951, on the roof of the tower of the University of London Senate 
House, at a height of 90 metres above sea level. The second set of experiments 
was carried out from July to September 1949 in the Sphinx observatory of the 
Scientific Station, Jungfraujoch, at a height of 3,560 metres above sea level. For 
both the sea-level and high-altitude measurements the apparatus was set up in 
a hut which had a roof of light material. In neither situation was there any large 
amount of material screening the apparatus at angles greater than 20° to the 
horizontal. 
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§3. EXPERIMENTAL RESULTS 
3.1. General Remarks 


Most of this experimental work has consisted of measuring the anti- 
coincidence rates —A+B+C+D as a function of the thickness of material at 
a,bandc. The value of the telescope coincidence rate B + C+ D will not normally 
be quoted; it is sufficient to note that, with 10 cm. of lead at a and no absorber 
at 6 or c, this counting rate was 1,050 per hour at sea level and 2,100 per hour 
at the Jungfraujoch. 

Two corrections have been made to the results for the anticoincidence rates. 
Firstly, variations in the barometric pressure are assumed to have changed the 
rates by —0-65°, per millibar. This figure was obtained from the measured 
absorption of the radiation in air. Secondly, allowance has been made for those 
counts which were due to the leakage of mesons through the anticoincidence 
screen. The method of determining this residual background rate will be 
discussed later. It was about 7°% of the maximum anticoincidence rate at sea 
level and 3% at the Jungfraujoch. . 


3.2. The Absorption of the Non-Ionizing Particles 


Measurements of the absorption in lead at sea level and in lead and carbon ~ 
at the Jungfraujoch were made by varying the thickness of the absorber a. In 
addition, by comparing the rates at the two altitudes the absorption in air can be 
calculated. 

Table 1 
a (cm.) Rate per hour Absorption length 
Jungfraujoch 


b=15 cm. 10 Pb 18-3+0-4 
paraffin wax 


400+ 40 gm/cm? Pb 


17°5 Pb 14:-4+0-4 
c—O0 25 Pb 12:7+0°4 
10 Pb+ 
36 C TSO e5 190+ 35 gm/cm? C 
b= 7-SicmaC 10 Pb 15-3+0-4 
0 25 Pb 10-3£0-5 } 4305 0 or eee 
b—A5 cra. 10 Pb Slee 0)98} 
paraffin wax 380 + 100 gm/cm? Pb 
c=4:1 cm. Pb DNS; \218) 3°3+0:3 
Sea Level 
b=5 cm. Pb 10 Pb 2:05+0-10 7 
p= 20 Pb 1-75 40-09 700 + 200 gm/cm? Pb 
30 Pb 1-53+0-:09 | 
b=0-9 cm. Pb 10 Pb 2:13+0:15 
P20 20 Pb 1-32+0-14 300-+ 50 gm/cm? Pb 
30 Pb 1-214+0:12 


Corrected results for the absorption in lead and carbon are given in Table 1. 
The values for the absorption length have been calculated on the assumption 
of exponential absorption. It is seen from the Table that the absorption in lead 
could be the same at the two altitudes and is not markedly dependent on the nature 
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of the materials at 6 and c. The mean values for the absorption lengths of 


_ the non-ioniziny particles are found to be 400+25 gm/cm? in lead, and 


190 + 35 gm/cm? in carbon. 

In these experiments the area of the absorber was always sufficient to cover 
the solid angle defined by the counter telescope BCD. This means that the 
geometric conditions were satisfactory provided that the non-ionizing particle 
always produced an ionizing particle which maintained roughly the same direction. 

The corrected results which are relevant to the absorption in air are presented 
in Table 2. 


Table 2. (4—10 cm. Pb) 
c(cm.) Rate/hr. at Rate/hr. at 


oxen) sea level Jungfraujoch Rae 

15 paraffin wax 0 2:-47+0-10 18-3+0-4 7440°3 
TOE 0 2:14+0-15 15-3+0-4 PDs 
5 Pb 0 2-00 + 0-10 15:0+0-8 D2 0 
0 0 1-00 + 0-07 7:9+0°5 LI40-9 
15 paraffin wax 0-9 Pb 1-05+0-09 10-0+0°5 9-5+1-0 
15 paraffin wex Dep 136 0-76+ 0-06 PQERV SS 10-4+0°9 
Xe DDN 0:59+0-08 6:140°3 10:34+1-4 


The last column gives the ratio of the rates at the two altitudes for various values 
of b and c. The particular material at b appears to be unimportant but there 
is a slight increase in the ratio when there is an absorber at c. Again assuming 
exponential absorption and noting that the thickness of air between the two 
stations was 340 gm/cm?, the mean absorption length in air becomes 
170+2 gm/cm? with c=0, and 147 +4 gm/cm? with c40. 

The classification of the results into the two classes c=0 and c#0, is the 
simplest possible one. It is more likely that there is a systematic dependence of 
the absorption length on the thickness of material at c, but the data are much too 
scanty to establish this. Clearly the case c40 means that the apparatus was only 
detecting the more energetic particles. 

The errors given above are the statistical standard deviations. An estimate of 
the experimental error is less certain. For although the same apparatus, including 
mostly the same counters, was used, there is the possibility that the geometrical 
arrangements were not exactly the same. From the observed ratio of meson 
rates, as compared with the results of Rossi et al. (1940 a), this error cannot have 
been greater than a few percent. ‘Therefore, allowing for the experimental error, 
the final value for the absorption length for all the non-ionizing particles is 
170 +5 gm/cm?. 


3.3. Production and Absorption of the Secondary Particles 


Table 3 gives the anticoincidence rates for various materials at b and different 
thicknesses of lead at c. The masses of the lead, carbon and carbon content of 
the paraffin wax at b were identical. ‘The carbon and lead were used in the 
form of layers spread out to occupy the same volume as the paraffin wax, thus 
minimizing any geometric effect due to change in position of the producer 
material. In each case the thickness was sufficiently small for the absorption of 
the primaries to be neglected. It is seen from Table 3 that equal masses of 
carbon and lead are about equally effective in producing secondaries and that 
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the hydrogen present in the paraffin wax makes a significant contribution to the 
anticoincidence rates. 
The comparatively high counting rates when no material was placed at b need 
to be discussed. It seems possible to explain them in terms of the matter present 
in the walls of the counters, supports, etc. It is estimated that the average mass of 
matter between the active volume of the anticoincidence counters and the top bank 
of the telescope was several grammes per square centimetre. Also the geometric 
conditions were favourable for recording any event which occurs, for example, 


Table 3 
a (cm. Pb) b (cm.) ¢(cm.) 
Sea Level 0 0:9cm. Pb. 2:2cm. Pb. 4:1 cm. Pb. 
10 15 par. wax. 2:47+0:10 1:054+0:09 0°76+0:05 0-35+40-07 
10 756 2:13+0:14 — 0:53+0-08 — 
10 0:9 Pb 2-13+0:15 — 0:42+ 0:08 — 
10 0 1:02+0:07 0:38+0:07 0:234+0:05 0:14+0-06 
Jungfraujoch 
10 15 par. wax. 18 3+0:-4 10:3+0:5 7:9+0°3 5-1+0°3 
10 TASE 15:3+0-4 — 6:1+0°3 3-740:3 
10 0 7:9+0°5 — 2:3+0°3 1:6+0-2 
25 15 par. wax. 12°7+0°5 — — 3°3+0°3 
aS a=l0cm.Lead | ab) a-l0cm Lead | 
wet (5 


so 


“ 
16 
2 
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Figure 2. Absorption of secondaries. 


in the glass of the top bank of telescope counters. Further measurements were 
made at thicknesses of c greater than those shown in Table 3. From the 
appreciably constant rates which were recorded it was concluded that this residual 
background was due to a leakage of mesons through the anticoincidence screen; 
in this way it was possible to estimate the value of the necessary correction, as 
mentioned in §3.1. 

It is obvious from ‘Table 3 that the secondary particles are strongly absorbed ina 
few centimetres of lead. The results are shown graphically in Figure 2, where the 
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thickness of absorber has been converted into the equivalent mass of lead. In 


_ doing this, allowance has been made for the finite amount of material in the tele- 


scope even when c =0 and for self-absorption in the producer layerb. ‘The absorp- 
tion is seen to be similar in all cases but is rather greater at sea level. ‘This 
corresponds to the greater absorption in the air of the more energetic primaries, an 
effect which was noted in the previous section. 


3.4. Multiplicity of the Secondary Particles 


During part of the sea-level experiments, an additional circuit indicated whether 
the centre counter of the trays B or C was discharged simultaneously with either of 
the outer two counters of the same tray. Anticoincidences in which this occurred 
will be referred to as multiple events; the results, given as a fraction of the total 
number of anticoincidences, are shown in Table 4. 


Table 4. (a=10cm. Pb) 


Multiples Multiples 


c (cm. Pb) b (cm.) Total (%) 


in B in C 
0 15 boca wax Hiss Fash 3840-7 
open ee ae see Nil 1940-7 
2 Ue a emo ie SOROS 
ee cath) 08 


The proportion of multiple events does not appear to depend on which tray is 
‘subdivided nor on whether b=0 or 15cm. paraffin wax. It is possible that there 
are rather less multiple events for large values of c. ‘The average ratio from all 
the experiments is 2-7 +0-4%. 

The possibility that this result can be explained in terms of knock-on electrons 
has to be considered before the multiplicity can be regarded as a genuine effect. 
Indeed, when the anticoincidence counters were switched off and the telescope 
recorded fast mesons, the measured ‘ multiplicity’ was 2:8°%. ‘The similarity 
between these figures must, however, be accidental. For it can be shown that any 
particle which is absorbed as easily as the secondaries in these experiments cannot 
knock on an electron with sufficient energy to penetrate the walls of the counters. 

It therefore seems that the multiplicity is a genuine effect. From the fact that 
the effect did not depend on the amount of material at b it can be concluded that two 
or more secondaries are produced simultaneously and not in successive collisions of 
the non-ionizing particles. 

The measured multiplicity of 2:7°/, can be only a fraction of the true multi- 
plicity, as the geometrical conditions for recording two particles, without dis- 
charging any of the anticoincidence counters, are very stringent. It seems likely 
that more than one secondary particle is produced in at least 10°% of the inter- 
actions of the non-ionizing particles, and that the proportion of multiple events 
might well be still higher. 
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3.5. Association of the Non-Ionizing Particles with Extensive Air Showers 


During some of the experiments at the Jungfraujoch an attempt was made to 
determine whether any of the anticoincidence events were accompanied by exten- 
sive showers. ‘T'wo trays of counters, each having an area of 700 cm?, were placed 
in the same horizontal plane as the anticoincidence set and distant 2-5 metres from 
it and each other. 

It was found that out of 1,220 anticoincidences, with various values of 
a, b and c, none were accompanied by the simultaneous discharge of both the 
extensive trays. In the same period the expected number of accidental 
coincidences was ():2. 

It can therefore be concluded that less than 0-5 °%/ of the non-ionizing particles. 
are part of extensive air showers with a density greater than one particle per 
1,000 cm*. This conclusion would not apply to dense penetrating showers. 
which, by discharging the anticoincidence screen, would prevent the detection 
of any non-ionizing particles accompanying them. 


3.6. Local Production of the Non-Ionizing Particles 


An experiment was carried out at the Jungfraujoch to determine whether 
any of the non-ionizing particles were produced by ionizing particles interacting — 
in the absorber a. For this purpose an additional tray of counters was placed. 
above a. It was found that 10+2-5°% of the anticoincidences occurred 
simultaneously with the discharge of this tray of counters. In these experiments. 
a was varied from 10 to 30 cm. of lead, but there was no marked dependence on 
this thickness. A similar effect, of roughly the same relative magnitude, was. 
subsequently found to occur at sea level. 

It seems necessary to interpret these results as evidence for the local production. 
of a fraction of the non-ionizing particles. A possible alternative explanation, 
in terms of the simultaneous arrival of ionizing and non-ionizing particles, is- 
rendered unlikely by the negative result of the extensive shower experiment. 


3.7. Measurement of the Angular Distribution 


The counter assembly was modified for this experiment and is shown in 
Figure 3. The entire counter assembly was so mounted that it could be rotated 
inside the fixed lead screen a. ‘The axis of rotation was approximately the same 
as the axis of the counter B and pointed E-W. 


En eae 


Longitudinal Section 10cm. Transverse Section 


Figure 3. Angular distribution apparatus. 


At sea level, the telescope counting rate B+ C+D was about 400 per hour; 
the angular distribution for all the ionizing particles was found to agree, except 
at large angles, with the expected cos?@ law. The results for the anticoincidence 
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rates are presented in Table 5. Neither the coincidence nor the anticoincidence 
rates showed any marked N-S asymmetry. 

It may at first sight seem surprising that this much smaller apparatus should 
give anticoincidence counting rates comparable with those obtained previously. 
There are two reasons for this: firstly, the minimum amount of material in the 
telescope is about 25% less than before, because there is no longer a supporting 
framework of wood and tin-plate; secondly, the geometric conditions are more 
_ favourable for the detection of the secondary particles. 


‘Table 5. (a=5 cm. Pb, b=5 cm. C) 
6 (deg.) 0 30 60 90 


c=0 2:22-0:08 1:893:0:07 ~§1:29+0-:06 0:86 0-05 
c=1:9cm.Pb 0:-40+0-03 0:3240-025 0-:214+0-02 0-09+0-01 


If it is assumed that the angular distribution of the secondaries is given by 
cos” §, then the required value of n, for values of 6 between 0 and 60°, is 
1:0+0-2. The distribution does not seem to depend greatly on the thickness 
of absorber in the telescope at c, except that there is a larger proportion of low 
energy particles at directions close to the horizontal. It is unlikely that the 
secondaries maintain the direction of the non-ionizing primaries, so that the 
value of x for the latter must be greater than 1-0. 


§4. DISCUSSION 
4.1. Nature of the Non-Ionizing Particles 


Despite the comparatively high intensity of photons present in the cosmic 
radiation, it is easily seen that they cannot be responsible for the present results. 
In particular, the measured ratio of the absorption in lead and carbon (Table 1) 
was about one to two, whereas the absorption of photons is at least ten times 
greater in lead than in carbon: Many of the other results also serve to show that 
the photon hypothesis is untenable. 

Any contribution from high energy neutrinos, which may have been produced 
during the decay of mesons, is extremely unlikely, as the interaction cross section 
is thought to be exceedingly small. 

Recently the existence of neutral mesons in the cosmic radiation has been 
established (Carlson et al. 1950), but their lifetime is extremely short. The 
apparatus used in the present experiments is incapable of detecting any particle 
whose lifetime is insufficient for it to pass through the anticoincidence counters. 
‘Thus it would be necessary to postulate a second and more stable type of neutral 
meson. 

The neutral V-particles of high mass which have recently been found in 
penetrating showers are believed to have a lifetime just sufficient for them to be 
detected by this apparatus (Armenteros et al. 1951), but their occurrence in showers 
means that the anticoincidence arrangement will invariably discriminate against 
them. 

A comparison of the results with those to be expected if the particles were 
neutrons can be made with the help of results obtained from accelerator work 
on their behaviour. In particular it is possible to compare the absorption with 
the cross sections given by Fox et al. (1950) for 280 Mev. neutrons. The work of 
DeJuren and Moyer (1951) has shown that the absorption is constant for neutron 
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energies between 180 and 280 mev. In making this comparison it is important 
to note that in the cosmic-ray experiments the conditions are those which nuclear 
physicists specify as ‘ bad geometry’. This means that the comparison must be 
made with the znelastic cross sections, as the elastic collisions will have little effect 
on the measured absorption ranges. 

The comparison is set out in Table 6, which shows that the absorption of 
the cosmic-ray non-ionizing particles is slightly less than that expected for neutrons. 
of a few hundred mev., and that the dependence on atomic number is similar. 

In view of this relatively close agreement it seems worth while to consider 
whether the other experimental results are consistent with the neutron hypothesis. 
Firstly it can be concluded from the results on the absorption of the secondaries, 
which would now be identified as protons, that the necessary neutron energies. 
are greater than 100 Mev. ‘That these neutrons should sometimes produce more 
than one secondary ina single interaction is notimprobable. The local production 
of the non-ionizing particles would be explained as the nuclear interaction of high. 
energy protons leading to the ejection of neutrons. 


Table 6. Comparison of Cross Sections (expressed in barns) 


Carbon Air Lead 
180 to 280 Mev. 0:14+0-01 0:16+0-01 1-43+0-08 
neutrons 
Cosmic-ray 0-11+0-02 0:140+0:-005 0:87+0:06 
particles 


It is now usual to assume that the primary protons of the cosmic radiation are not 
absorbed catastrophically but give rise to a number of high energy nucleons, which 
in turn can make further collisions. ‘This picture has been developed theoretic- 
ally in some detail, particularly by Messel (1951). The particles observed in 
the present experiments can reasonably be interpreted as the lower energy neutrons 
and protons which are produced in these nucleon cascades. ‘This view is supported 
by the experimental results on the local production of the non-ionizing particles. 
Similarly the results for the angular distribution show that most of the neutrons 
must have been produced within a few interaction lengths of the apparatus, since 
even an isotropic production at high altitudes would give rise to a sharply colli- 
mated distribution at sea level. ‘The complicated nature of the events occurring 
in the nucleon cascades may explain the slight difference between the absorption 
of the non-ionizing particles and neutrons which is shown in Table 6. 

As the neutron hypothesis appears capable of explaining most of the observed 
effects, it will be developed more quantitatively in the next sub-section. 


4.2. The Energy Spectrum and Intensity of High Energy Neutrons 


On the assumption that the secondary particles are protons it is possible to 
convert the range spectra of Figure 2 to integral energy spectra. In the case of the 
paraffin wax—carbon difference the process occurring is the simple (n, p) one, where 
it is known that, owing to the influence of exchange forces, the proton usually 
acquires nearly all or else very little of the neutron’s energy. In the former case 
there is a strong correlation between the direction of the neutron before collision 
and of the proton afterwards. ‘Thus the proton spectrum should correspond 
directly to the neutron spectrum. If it is assumed that the differential energy 
spectrum is of the form N(E) dE=E~’1dE, then the value of y obtained is 1:0+0-5- 
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Similarly, the spectrum for all the secondary particles from paraffin wax is given by 
py=1-6+0-1.  Itis possible that these two spectra are really the same, although it 
would not be surprising if the energy spectrum of the particles ejected from carbon 


| nuclei was rather steeper. 


In estimating the neutron intensity it is more satisfactory to use the results for 
the paraffin wax—carbon difference. ‘The measured value for this difference at the 
Jungfraujoch was 3-0+0-6 per hour. After correcting for the absorption in the 
10 cm. of lead at a above the apparatus, this becomes 4-0 + 0-8. 

The important constants of the apparatus are : area of producer layer b=150 cm?, 
- thickness of b=12gm/cm? of paraffin wax, composition of wax (CH,),, mean 
solid angle subtended by telescope at b, Q. =0-12 steradians. 

It will be assumed that the angular distribution at the Jungfraujoch is not very 
different from that measured at sea level. The partial (n, p) cross-section, o (8), 
can be found from the results of Hadley et al. (1949) and Kelly et al. (1950). Values 
were interpolated to correspond to a mean energy of about 200 Mev. for the cosmic- 


ray neutrons. n/2 
Then the total counting rate becomes | I(@)27sin9@NQa(0)d0. This 


0 
can be integrated numerically and equated to the measured anticoincidence 
rate. In this way it is found that the vertical intensity of neutrons must be 14cm? 
sterad-! hour-!. The ratio of penetrating non-ionizing to the total ionizing 
_ radiation at an altitude of 3,560 metres is therefore 11°. The equivalent figure at 
sea level is found to be 4:5 %. 

The absolute accuracy of these estimates is not very good, chiefly because the 
standard deviation is 20%. There is a further uncertainty of about 10% in 
assessing the geometrical constants of the apparatus. For these reasons it has not 
seemed necessary to take account of the energy spectrum and the variation of the 
cross section with energy. 


4.3. Comparison with other Measurements of the Intensity 


The only comparable measurements which have been made with an anticoin- 
cidence arrangement are those of Janossy and Rochester (1943). ‘Their results 
are compared with those of the present experiment in Figure 4; the rates given by 
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Figure 4. Comparison with previous results. 


Janossy and Rochester have been multiplied by a normalizing factor of 1:5, which 
was obtained from the ratio of the telescope counting rates in the two experiments. 
The agreement is seen to be satisfactory, despite somewhat different geometrical 
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arrangements. However, Janossy and Rochester estimated that the intensity of 
non-ionizing particles at sea level was not greater than 0-14°%, which is only one 
thirtieth of the value calculated above. The difference can be accounted for as 
follows: (i) Janossy and Rochester always had 2-2 cm. of lead at c which, as shown 
in Tabie 3, reduces the rates to about one quarter ; (ii) they used 5 cm. of lead at b; 
again because of the strong absorption of the secondaries, not more than one third 
of this would be effective ; (iii) they assumed on theoretical grounds that the collision 
length in lead for the non-ionizing radiation was 10cm.; the absorption length is 
now found to be 35cm. and it is unlikely that there will be more than one high 
energy secondary produced in this distance. ‘These corrections bring the two 
sets of results into good agreement. 

Estimates of the intensity of particles producing stars, which are also assumed 
to be neutrons, have been made by George and Jason (1949). Working in terms of 
the integrated intensity, their results give values of 11° and 3-5 % at the same two 
altitudes. "The agreement here is rather better than would be expected as it is 
now known that about one-sixth of the stars are due to ionizing particles (Brown 
et al. 1949). On the other hand the mean energy for the smallest size star 
detected is larger, about 150 Mev. (Brown et al. 1949). ‘These two corrections will 
apply in opposite directions, which may be responsible for the final close agreement. 
It is satisfactory that the ratio for the two altitudes agrees with that found here in the 
case c<0, i.e. the higher energy neutrons. 

From similar observations Page (1950) concluded that the intensity of the 
star producing radiation was only about one half of that given by George and Jason. 
The difference is a result of assuming that the collision length for star production 
is half the absorption length, on the grounds that the latter corresponds to a 
cross section which is only one-half of the geometric area of the nucleus. The 
Berkeley results have now shown that although the total absorption cross section 
for neutrons of 270 Mev. is roughly geometric, about 50°% of the collisions cause 
only elastic scattering (DeJuren 1950). Hence the cross section for star pro- 
duction is half the geometric one, and the assumption made by Page is incorrect 
for the moderate energy neutrons which produce most of the cosmic-ray stars. 

Another estimate of the intensity has been made by Tongiorgi (1949 b) from 
measurements on low energy neutrons. It was concluded that the flux of high 
energy neutrons was 20% and 5% at altitudes similar to those used in the present 
experiments. In other experiments Tongiorgi (1949 a) has found that about 
2°% of the particles in extensive air showers are capable of producing slow neutrons. 
Compared with the total number of neutron producing particles, however, the 
number occurring in extensive air showers is only 0-:03°%. Hence this result is 
not in disagreement with the conclusion reached here: namely, that less than 
0-5 °% of the non-ionizing particles are accompanied by extensive air showers. 

Rogozinski (1948) reported the existence of a horizontal neutral component. 
No absolute value was given for its intensity, but it seems possible that the 
observed effects were due to neutrons, as it has been found in the present 
experiments (§ 3.7) that many of these have directions making large angles to the 
vertical. 


4.4. Production of High Energy Secondaries 


The results of §3.5 enable the relative number of secondaries from hydrogen 
and carbon to be compared directly. For the high altitude results the corrected 
tates were carbon 9-4+0-7 per hour, paraffin wax—carbon 3-0 + 0-6 per hour. 
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Hence it follows from the composition of the paraffin wax that the cross 

| section for the ejection of a fast ionizing particle from carbon is about six times 
7 the(n, p) cross section. Substituting the known value of the latter, it is found that 
7 the partial cross section for the forward ejection of a high energy secondary 
# from carbon is 120+40 millibarns per sterad. This value is comparable to the 
total absorption cross section for neutrons by carbon, and may therefore seem 
surprisingly high. It is, however, confirmed by the work of Hadley and York 
(1950) who found that the partial cross section for the forward knock-on of 
} ‘secondaries from carbon by 90 Mev. neutrons was about 100 millibarns per sterad. 
| They also showed that, for light nuclei, the angular distribution of the 
secondaries was strongly peaked in the forward direction, so that the total cross 
section remains reasonably small. 

The emission from cosmic-ray stars of protons with an energy too great to 
be explained on the evaporation theory is well known. Detailed results referring 
to the stars observed in nuclear plates exposed at mountain altitudes have been 
given by Brown e¢ al. (1949) and Bernardini et al. (1950). Their results are 
| expressed in rather different forms, but indicate that the mean number of high 
| energy secondaries per star is between 0-5 and 1:0. In about 20% of the stars 
_ more than one fast secondary is produced, in agreement with the results obtained 
here from the multiplicity experiments. 


e §5. CONCLUSIONS 

1. The mean absorption lengths of that part of the non-ionizing component 
capable of penetrating more than 10 cm. of lead are 170+5 gm/cm? in air, 
190 +35 gm/cm? in carbon and 400 +25 gm/cm? in lead. 

i 2. The average range of the secondary particles produced by the non-ionizing 
~ component is about 30 gm/cm? of lead. 

3. The number of secondary particles produced in equal masses of different 
materials is roughly constant. 

4. 'T'wo or more secondary particles are produced simultaneously in more than 
10% of the interactions of the non-ionizing particles. 

5. The angular distribution of the secondary particles produced in carbon is 
approximately proportional to cos 0. 

6. Less than 0:5% of the non-ionizing particles at 3,560 metres are part of 
extensive air showers with a density greater than one particle per 1,000 cm?. 

| 7. About 10% of the non-ionizing particles are produced by ionizing particles 
in the lead above the apparatus. 

8. All the experimental results are consistent with the hypothesis that the 
majority of the non-ionizing particles are neutrons, and are inconsistent with the 
behaviour of any other known particles. 

9. The differential energy spectrum of the neutrons is given by E~’! where 
y=1-0+0-5 in the energy range 100<E<250 mev. 

10. The vertical intensity of neutrons with energy greater than 100 Mev. is 
11% of that of the total ionizing radiation at a height of 3,560 metres, and 4-5 % at 
sea level. 
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PETTERS AO. ThE Ebi 
The Long-Wave Limit of Infra-Red Photoconductivity in PbSe 


Published measurements on the spectral distribution of photoconductivity in PbSe 
(Roth 1938, Blackwell et al. 1947, Moss and Chasmar 1948, Starkiewicz 1948, Milner 
and Watts 1949) indicate that the long wavelength ‘limit’ is not greater than, and is 
possibly identical with, that of PbTe. The ‘limit’, as measured, indicates a 50% drop 
in sensitivity at about 5-0 when the cell is at a temperature of 90°K. We have now 
established that the ‘limit’ of photoconductivity in PbSe is in fact at a significantly longer 
wavelength than that of PbTe provided the layers are prepared in a suitable manner. 

Recent work by R. V. Jones and his colleagues (Paul, Jones and Jones 1951) has shown 
that certain single crystals of PbS are relatively transparent beyond 3y, a sharp absorption 
edge occurring at this wavelength. This observation has been confirmed by us and in 
addition the absorption spectra of PbSe and PbTe crystals grown in this laboratory (Lawson 
1951) have also been measured over the temperature range 20° to 600° kK. This work will 


be reported in another publication. It is found that the position of the absorption edges — 


in PbS and PbTe is in good agreement with the observed limits of photoconductivity. 
For PbSe, however, the absorption edge occurs at about 4:7 at room temperature and 
6:7 at liquid oxygen temperature. These figures indicate that the long-wave photo- 
conductive limit of PbSe has been set in the past at much too low a value. This has now 
been confirmed by direct measurement. The photoconductivity of a PbSe photo-diode 
consisting of a PbSe p-type crystal and a tungsten cat’s whisker, has been studied. The 
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long-wave limits, defined as above, are found to be at 4:74 and 6:8 at room temperature 
and 90° k. respectively. These values are in good agreement with the absorption spectra. 
In addition it has recently proved possible to make photoconductive layers of PbSe with 
| long wavelength limits (50% down from maximum) at 4-8 at room temperature and 
7:14 at 90°. The spectral sensitivity curve for such a layer is shown in the Figure. The 
4 100 
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Spectral response of lead selenide cell. 


technique of layer sensitization has not, as yet, been studied in any detail and PbSe cells 
produced up to the present are of relatively low sensitivity. There is no obvious reason 
| why this deficiency should not be overcome and, if so, a cooled PbSe cell should be very 
valuable for infra-red spectroscopy in the 5—74u region. 

We should like to acknowledge our indebtedness to Messrs. A. S. Young and D. A. H. 
Brown, of this Establishment, for providing the sapphire windows and electronic equipment 
used in the above measurements. 

This letter is published by permission of the Chief Scientist, Ministry of Supply. 


Note added in proof. A later cell with a periclase window gave ‘limit’ at 8-1 at 20°x. 


Telecommunications Research Establishment, A. F. GIBSON, 
Great Malvern, Worcs. W. D. Lawson. 
20th August 1951. T. S. Moss. 
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Angular Distributions of the Neutron Groups from the 
Reaction ’C(d,n)'*N 


Using the same technique as described in our previous communication (El-Bedewi, 
Middleton and Tai 1951), we have observed, at a bombarding energy of 8 Mev., three neutron 
groups from the reaction C(d, n)!*N which correspond to the transitions to the ground 
state and two excited states of *N at 2:38+0-05 Mev. and 3:53+0-05 Mev. respectively. 
The angular distributions of these groups, plotted in the centre-of-mass system, together 
with the theoretical curves calculated by Newns (private communication) on a deuteron 
stripping theory similar to that proposed by Butler (1950) are shown in Figures 1, 2 and 3. 

- For the first and second excited states the best agreement between the theoretical and 
_ the observed angular distributions was achieved by assuming the angular momentum of the 
absorbed proton 1, equal to 0 and 2 respectively. It follows that the first excited state of 
18N must have spin } and even parity and the second excited state has spin 5/2 or 3/2 and 
even parity. Apart from the parity of the second excited state, which is in contradiction to 
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that given by Hornyak, Lauritsen, Morrison and Fowler (1950), the results agree with | 

previous predictions (Hornyak et al. 1950, Goldhaber and Williamson 1951). 
For the ground state the general appearance of the angular distribution suggests | 

1)=2, but this interpretation would involve an excessively large radius of interaction to be 
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used in the theoretical calculation. Fair agreement was obtained, as shown in Figure 1, 
by assuming /,=1, which would indicate this state to be Ps, or Pi, state. This agrees with 
the prediction of P4/. state from Mayer’s model (Mayer 1950) and also with the results 
obtained for the mirror nucleus °C, recently reported by Rotblat (1951). 


Nuclear Physics Research Laboratories, F. A. EL-BEDEwI. 
University of Liverpool. R. MIDDLETON. 
13th August 1951. CAaveIATs 
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The Reaction ’C(y, 3a) and a 16°9Mev. Level of *Be 
Emitting Alpha-Particles 


‘The stars formed in a nuclear emulsion by the 1C(y, 3x) reaction have already been 
studied for gamma-ray energies up to 25 mev. (Goward, Telegdi and Wilkins 1950, 
Wilkins and Goward 1951). The purpose of this letter is to report a study of the reaction 
using Ilford E1 and C2 emulsions (200 microns thick) exposed to 70 Mev. bremsstrahlung. 
The “C(y, 3a) stars were first selected from the considerable variety of stars observed 
(many of which are as yet unidentified) on the basis of appearance combined with a 
momentum balance criterion (Goward and Wilkins 1950). 180 stars selected in this way 
proved to be sufficient to determine the variation of the ##C(y, 3a) cross section as a function 
of gamma-ray energy, E,, up to E,~40 Mey. A second maximum in the cross section 
appears at E, ~27 mev., of approximately three-quarters the amplitude of the maximum 
already known to exist at 18:5 Mev. At 23:5 Mev. and 33 mev. the cross section falls to 
half the value at 27 Mey. 

Further analysis of the selected stars also reveals a (16:9+ 0:3) Mev. excited level of *Be, 
which decays by emitting two alpha-particles and whose existence has not been reported 
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previously. This level may be shown to play a predominant part in the formation of the 
stars at L,,>26 Mev. For each star (with E,, >26 mev.), three values were calculated for 
the energy E* released when an assumed intermediate *Be nucleus disintegrates into two 
alpha-particles. For eight of the stars an E* value~0-1 Mev. was obtained, indicating 
probable transition through the ground state of *Be; the remaining 52 yielded the E* 
histogram shown in the Figure. The curve fitted to the histogram was calculated by 
assuming that all 52 stars involve a 16-9 mev. excited level of ®Be, of apparent width 1-0 Mev. 
at half maximum, and the excellent fit obtained indicates the predominance of transitions 
through this new level. The estimate of 16:9 Mev. excitation energy is unlikely to be more 


Number of Values 


than 0:3 Mey. in error, and allowance for limitations of energy resolution (arising 
particularly from straggling) enables the real width at half maximum to be estimated as 
less than 0-8 Mev. ‘The nearest levels of *Be already known lie at 17-6 Mev. and 18-1 mev. 
excitation (Hornyak et al. 1950), but neither of these levels emits alpha-particles. 

The 60 stars with E, >26 mev. hardly provide a sufficient basis for detailed figures of 
the proportions in which the various ®Be levels occur. Eight stars probably involve the 
ground state, and all the remaining 52 are consistent with either the 3 Mev. or 16:9 Mev. 
excited levels (roughly in the proportion 5 : 47), but other levels are not definitely excluded. 

We are indebted to Dr. H. C. Pollock for irradiating our emulsions at the General 
Electric Co., Schenectady, U.S.A., where Dr. A. M. Gurewitsch is carrying out experiments 
similar to that reported here. 


Atomic Energy Research Establishment, J. J. WILkKINs. 
Harwell, Didcot, Berks. F. K. Gowarp. 
27th August 1951. 
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The Thermal Conductivity of Cadmium in a Magnetic 
Field at Low Temperatures 


In the course of measurements of the thermal conductivity of a number of metals 
(Au, Mo, Rh, U, Cd) at low temperatures the effect of a magnetic field was investigated. 
While in the other metals this produced only a very small decrease in the thermal 
conductivity, in the case of cadmium the effect was very marked even in the relatively 
small fields used. The specimen was a rod cast from H.S. cadmium (Lab. No. 10,547) 
whose impurity content was less than 1 in 10°. 

The graph shows the relation between the ratio Ry/Ry and H, where Ry, is the thermal 
resistance at a given temperature in a field H and R, is the thermal resistance at this 
temperature in zero field. For H=3,800 gauss, when the field was transverse to the 
temperature gradient, the thermal resistance increased five tinies at 4-4° kK. and over seven 
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times at 2:2°k. As can be seen from the graph, the effect of a longitudinal field was very 
much smaller. Except for the lowest fields the relation is approximately linear and there 
is no sign of any saturation in the highest fields used. 


BE 


Field (kilogauss) 


These results are of interest when compared with the experiments of Milner (1937) 
on the electrical resistance of cadmium, in which he found a very large change of resistivity — 
with field at low temperatures. A more detailed description of the work together with 
further results will be given at a later date. 


Clarendon Laboratory, K. MENDELSSOHN. 


Oxford. H. M. RosENBERG. 
31st August 1951. 
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ABSTRACTS FOR SECTION B 


The Diffraction Pattern of Cold Worked Metals : I—The Nature of Extinction, by 
W. H. Hatt and G. K. WILLIAMson. 


ABSTRACT. ‘The influence of crystal perfection on extinction is briefly discussed and 
the possible corrections given. Experimental results for chemically produced copper and 
annealed aluminium filings have been obtained with a Geiger counter spectrometer using 
monochromatic radiation. ‘The intensity measurements, and for copper the experimentally 
determined particle size, were inconsistent with the views that primary extinction occurred, 
but could be interpreted whol:y in terms of secondary extinction, the corrected f values 
agreeing within +2 with accepted values. Evidence to show that annealed metals 
contain residual lattice strain is advanced, in accord with the dislocation theory, and it is 
suggested that secondary extinction will predominate over primary extinction in all crystals 
containing such dislocations. 


The Diffraction Pattern of Cold Worked Metals: II—Changes in Integrated 
Intensity, by W. H. Hatt and G. K. WILLIAMson. 


ABSTRACT. Measurements have been made, using a Geiger counter spectrometer, of 
the changes produced by cold work on the integrated intensities in the Debye—Scherrer 
spectrum of aluminium. It is shown that cold work reduces the integrated intensity of all 
reflections, this decrease being exactly compensated by a rise in the diffuse background 
scatter, so that the total reflected intensity is unchanged. ‘These results are at variance with 
those of previous workers, but the discrepancy can be readily explained by the absence of a 
correction for extinction in previous work. Brief reference is made to measurements of 
line width and it is shown that the results support the dislocation theory of plastic 
deformation. 


Outer Brillouin Zones for Face-Centred and Body-Centred Cubic Lattices, by 
J. F. NicHo.as. 


ABSTRACT. In order to resolve some inconsistencies in the literature on the second 
Brillouin zones of cubic lattices, the shapes of these zones have been recalculated. The 
results, together with descriptions of the third and fourth zones, are given. 


Imaging Properties of a Series of Magnetic Electron Lenses, by G. LIEBMANN and 
E. M. Grab. 


ABSTRACT. ‘The paraxial lens data and the more important first order lens aberrations 
were determined for a series of symmetrical magnetic lenses of equal lens diameter D but 
variable air gap width S for a wide range of lens excitations. ‘The range of gap widths 
covered was S/D=0:2 to S/D=2. The field distribution within the lenses was measured to 
a high degree of accuracy with a resistance network analogue, and the lens data were then 
obtained through numerical trajectory tracing. ‘The results are represented as a series of 
graphs. ‘These results are representative of the most common magnetic electron lenses, and 
should have useful applications in lens design. An approximate formula for the focal length 
is f =~50(DS)1/?V,/(ND?. A few design examples are discussed. 


The Symmetrical Magnetic Electron Microscope Objective Lens with Lowest 
Spherical Aberration, by G. LIEBMANN. 


ABSTRACT. Using results recently published it is shown that each of the investigated 
symmetrical magnetic lenses has an optimum value of the lens excitation parameter, for which 
the spherical aberration constant is smallest. Taking account of the dependence of the 
maximum obtainable axial field strength Hy on the maximum field strength Hy in the pole 
piece gap, one finds that there is an optimum design of electron microscope objective lens 
of conventional type which gives a lowest value of spherical aberration and highest resolving 
power. A modification of this optimum lens as adapted for practical use is described. 
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The Photographic Action of X-Rays, by J. R. GREENING. 


ABSTRACT. ‘The variation of photographic blackening with intensity, exposure, the 
thickness of material round the film, the angle of incidence and the wavelength of x-rays, are 
examined theoretically and, with one exception, experimentally. A detailed theoretical 
investigation is made of the energy absorption in the silver bromide grains, consideration 
being given to the energy lost by scatter and fluorescence and to the electron interchange 
between the silver bromide grains and the surrounding gelatine and other material. It is 
concluded that for the majority of emulsions there will be no variation of blackening with 
X-ray intensity and that the photographic density will be proportional to the exposure at 
low densities. At very short wavelengths the blackening is critically dependent upon the 
thickness of material round the emulsion. Photographic density decreases with increase 
of angle of incidence of x-rays, the maximum decrease observed being about 17%. Within 
the wavelength range examined experimentally (0-020 to 0-354 a.) the photographic density 
was approximately proportional to the energy absorbed by the silver bromide. The 
density per réntgen varied by a factor of 42 in this range. The variation of density per 
roéntgen with wavelength and angle of incidence limits the applications of photographic 
methods of x-ray measurement. Methods by which emulsion manufacturers can reduce 
these variations are suggested. 


A Study of Rectification Effects at Surfaces of Germanium and Lead Sulphide, by — 
C. A. HocartH and J. W. GRANVILLE. 


ABSTRACT. Previously reported results of polishing and chemically etching the surfaces 
of germanium crystals are confirmed. The effects of vacuum heat treatment are investigated 
and it is found that treatment of this kind at temperatures between 500 and 900°c. can 
improve the rectification characteristics considerably. The treatment does not spoil the 
optical quality of a polished surface and in this connection may provide an advantage over 
chemical etching. Similar effects are obtained with deficit-conducting lead sulphide. 
Modifications to the surface structures of the crystals by various treatments are investigated 
by electron diffraction. 


A Vibrating Cantilever Method for the Investigation of the Dynamic Elasticity of 
High Polymers, by K. W. HILtier. 


ABSTRACT. A method is described for obtaining dynamic values of Young’s modulus 
at low frequencies (20 c/s. to 100 c/s.) by using the forced resonance of a vibrating cantilever. 
Provision has been made for keeping the apparatus at a constant temperature and a few 
measurements of the dynamic modulus are given for polythene at 30° c. and 40° c. The 
theory of the vibrations of a cantilever including viscous internal damping is discussed. 
The results continue previous work carried out at higher frequencies by a different method 
and indicate the continuity of the dynamic modulus—frequency relation between 
20 ke/s. and 20 c/s. for polythene. ‘ Static ’ measurements have also been carried out and 
agree with those published earlier. Even at the lowest frequency the dynamic modulus 
is still three times the “‘ static’ value, and it is concluded that lower frequencies should 
be investigated. 


A Simplified Shearing Interferometer, by R. L. Drew. 


ABSTRACT. A simple new design of interferometer based on the wave-front shearing 
principle is described and illustrations given of its application to the solution of a variety of 
‘through-type’ and autostigmatic testing problems. 
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